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I. Abstract for the project ('reprinted from the proposal 1 ) 

Mariner 9 UV spectrometer data would be reinverted for the ozone 
abundance, cloud abundance, dust abundance, and polar-cap albedo. The 
original reduction of the spectra ignored the presence of atmospheric 
dust and clouds, even though their abundance is substantial and can mask 
appreciable amounts of ozone if not accounted for (Lindner, 1988). The 
Mariner 9 ozone data has been used as a benchmark in all theoretical 
models of atmospheric composition, escape, and photochemistry. A sec- 
ond objective is to examine the data for the interrelationship of the 
ozone cycle, dust cycle, and cloud cycle, on an annual, inter-annual, and 
climatic basis, testing predictions by Lindner (1988). This also has 
implications for many terrestrial ozone studies, such as the ozone hole, 
acid rain, and ozone-smog. A third objective is to evaluate the efficacy 
of the reflectance spectroscopy technique at retrieving the ozone abun- 
dance on Mars. This would be useful for planning ozone observations on 
future Mars missions or the terrestrial troposphere. 

II. Summary of Research to Date 

11.1 Model development. 

The photochemical-radiative transfer model of Lindner (1985; 

1988) has been updated. The code has been modified to run on AER's 
computer system, has been updated to use new cloud and dust scattering 
parameters as determined by Clancy and Lee (1991), and has been 
updated with new improvements to the Discrete Ordinate Method by 
Stamnes et al. (1988). Briefly, the code runs in 3 parts. First a code has 
been written which sets up the base atmosphere, including photochem- 
istry, and then computes the wavelength-dependent atmospheric opacity, 


4 


single-scattering albedo and phase function. Figure 1 shows the atmo- 
spheric opacity as a function of wavelength for one atmosphere scenario. 
This opacity includes Rayleigh scattering (strongest at 1500 Angstroms, 
and decreasing rapidly with wavelength), CO 2 absorption (the prominent 
component of the opacity shortward of 2000 Angstroms), and ozone 
absorption (the "bump" in opacity between 2000 and 3000 Angstroms). In 
addition, there is a virtually wavelength-independent opacity of cloud 
and dust. It is this "bump" that leaves a characteristic signature of the 
ozone abundance in the measured UV spectra. The radiative transfer 
model then computes intensities as a function of wavelength based on 
this atmospheric opacity and scattering. Finally, a code has been writ- 
ten which converts these intensities into a ratio of UV spectra which 
can be compared to the data. 

ILZ Acquire Mariner 9 UV Spectrometer data . 

While at the MSATT workshop held in Boulder, Colorado Sept. 23- 
25, 1991, the principal investigator discussed the procedure for obtain- 
ing the Mariner 9 UV Spectrometer data with Steve Lee, the director for 
the node of the Planetary Database System (PDS) at the University of 
Colorado. This and subsequent discussions have indicated that the data 
are not on the PDS system as promised years ago, and are not expected to 
be on the PDS for several years, due mostly to budget cuts. Steve has 
offered to help us retrieve the original data, however, we have been 
denied access to this data by Charles Barth of the University of Colorado, 
despite our requests to him to release this 20 year old data. 

Nonetheless, there have been several good spectra published in Lane et 
al. (1973) and Wehrbein (1977), which are quite adequate for us to com- 
plete all of our objectives for this project. A sample data spectra is 
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Figure 1. Vertical opacity of the atmosphere for 30 pm-atm of 
ozone, no dust and no cloud. A winter polar temperature profile is 
adopted (see Lindner, 1988). 
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shown in Figure 2. 

11.3. Simulate UV Spectra. 

The theoretical model was run with a variety ot inputs to generate 
synthetic UV spectra. Parameters within the base atmosphere were 
varied as shown in Table 1 and synthetic spectra were calculated for 
each case. This is done to examine how changes in atmosphere affect 
ozone retrieval efficacy, affect errors in retrieved ozone abundance, and 
affect scatter in data. For ease of comparison to data (see Fig. 2), we 
have computed ratio of spectra to that at 20N latitude, where little 
ozone is present. Computed spectra for 20N latitude show little wave- 
length dependence. Figure 3 shows how spectra can vary for a single 
base atmosphere depending upon the viewing geometry of the spacecraft. 
Hence, comparison of synthetic spectra to observed spectra must prop- 
erly account for the correct geometry. Synthetic spectra for various 
scenario for base atmospheres are shown in Figures 4 through 7. 


TABLE 1 

Synthetic spectra calculated for the following ranges in parameters 


Ozone 

Ozone distribution 
Solar Zenith Angle 
Dust opacity 
Cloud opacity 
Ice albedo 
Look Angle 

Azimuth of Look Angle 


0 to 100 pm-atm. 

well-mixed, and concentrated down low 

50 to 85 degrees 

0 to 1 

0 to 1 

0.3 to 0.8 

0 to 90 degrees 

0 to 90 degrees 
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Figure 3L Ultraviolet spectrum measured by Mariner 9 at 57*N latitude on orbit 

144 (upper plot). To enhance the ozone absorption feature, this spec- 
trum was divided by one obtained at 20°N latitude on orbit 144, where 
ozone abundances are minimal. The ozone absorption feature for orbit 
202 is compared to laboratory data in the lower plot (laboratory data 
is offset to facilitate comparison) . Both figures are taken from Lane 
et al. (1973). 
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Wavelength (in Angstroms) Wavelength (in Angstroms) Wavelength (in Angstroms) 

Figure 3. Synthetic ultraviolet spectra, ratioed to synthetic spec- 
tra for 20N latitude, shown for various viewing geometries (viewing an- 
gle theta relative to the zenith and azimuthal angle phi relative to the 
sun). As with Lane et al. (1973; see Fig. 2), intensity is unitless and 
relative. 30 pm-atm ozone, no dust or cloud, a solar zenith angle of 75 
degrees, ice albedo 0.5, and uniform mixing ratio of ozone are assumed. 
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Figure 4. Synthetic ultraviolet spectra, ratioed to synthetic spectra for 
20N latitude. Viewing geometry in all cases is viewing angle of 60 
degrees, and azimuthal angle of 0 degrees. Atmospheric and surface 
parameters for the plots are as follows (showing ozone abundance, solar 
zenith angle (SZA), ice albedo, and cloud and dust opacity): 

4(a) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, No dust, No cloud 
4(b) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, dust 0.3, cloud 1.0 
4(c) Ozone 100 pm-atm, SZA 75 degrees, Albedo 0.5, dust 0.3, cloud 1.0 

1 0 
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Figure 5. Synthetic ultraviolet spectra, ratioed to synthetic spec- 
tra for 20N latitude. Viewing geometry in all cases is viewing angle of 
60 degrees, and azimuthal angle of 0 degrees. Atmospheric and surface 
parameters for the plots are as follows (showing ozone abundance, solar 
zenith angle (SZA), ice albedo, and cloud and dust opacity): 

5(a) Ozone 30 pm-atm concentrated low in the atmosphere, SZA 75 
degrees, Albedo 0.5, No dust, No cloud 

5(b) Ozone 30 pm-atm concentrated low in the atmosphere, SZA 75 
degrees, Albedo 0.5, dust 0.3, cloud 1.0 

5(c) Ozone 100 pm-atm concentrated low in the atmosphere, SZA 75 
degrees, Albedo 0.5, dust 0.3, cloud 1.0 

1 1 
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Figure 6. Synthetic ultraviolet spectra, ratioed to synthetic spectra for 
20 N latitude. Viewing geometry in all cases is viewing angle of 60 
degrees, and azimuthal angle of 0 degrees. Atmospheric and surface 
parameters for the plots are as follows (showing ozone abundance, solar 
zenith angle (SZA), ice albedo, and cloud and dust opacity): 

6(a) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, No dust, No cloud 
6(b) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, dust 0.3, No cloud 
6(c) Ozone 100 pm-atm, SZA 75 degrees, Albedo 0.5, dust 0.3, No cloud 
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2100 2400 2700 3000 3300 


Wavelength (in Angstroms) 

Figure 7. Synthetic ultraviolet spectra, ratioed to synthetic spectra for 
20N latitude. Viewing geometry in all cases is viewing angle of 60 
degrees, and azimuthal angle of 0 degrees. Atmospheric and surface 
parameters for the plots are as follows (showing ozone abundance, solar 
zenith angle (SZA), ice albedo, and cloud and dust opacity): 

7(a) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, No dust, No cloud 
7(b) Ozone 30 pm-atm, SZA 75 degrees, Albedo 0.5, dust 1.0, cloud 1.0 
7(c) Ozone 100 pm-atm, SZA 75 degrees, Albedo 0.5, dust 1.0, cloud 1.0 
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11.4. Compare synthetic and observed spectra. 

The synthetic and observed spectra compare quite favorably. In 
fact, the inclusion of dust and cloud absorption and scattering improve 
the synthetic fit to the slope in the data from 2100 to 2300 Angstroms. 
Comparison of Figures 4 and 5 shows how difficult it is to try to infer 
the ozone distribution from UV spectra. Two radically different ozone 
profiles were used, and yet the spectra exhibit similar behavior. 

Wehrbein (1977) examined the data in depth to retrieve the ozone distri- 
bution, and concluded that ozone is fairly well mixed in the atmosphere. 
Lindner (1988) also showed that ozone should be fairly well mixed based 
on theoretical simulations of atmospheric chemistry. 

11.5. Error in inferring On without considering dust/cloud. 

Ozone has been classically inferred from the spectra by fitting the 

depth of the absorption at 2500 Angstroms. We have summarized the 
depth of absorption in Table 2 for many cases we have tried. What can be 
most clearly seen is that for all the cases we have tried, the amount of 
ozone inferred from the UV spectra is underestimated by about a factor 
of 3 when the inversion uses simply ozone absorption (i.e. ignoring cloud 
and dust). This is seen in Table 2 in that approximately 3 times as much 
O3 is needed to produce the same depth in absorption for scenario in 
which we fully include the masking effects of cloud and dust to the case 
for which cloud and dust are ignored. This is particularly noteworthy in 
that even if the effect of cloud and dust is included in the retrieval, the 
uncertainty in cloud and dust opacity and scattering properties will still 
leave an uncertainty in retrieved ozone abundance of a factor of 2! 
Considering that the published ozone abundances retrieved from the 
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TABLE 2 

2550/3200 ANGSTROM RATIO 



CLOUD ALBEDO 1 03 DIST. 


mixed 

mixed 

mixed 

mixed 


mixed 

mixed 


mixed 

mixed 


low 

low 

low 


mixed 

mixed 

mixed 


mixed 

mixed 

mixed 


mixed 

mixed 

mixed 


mixed 

mixed 

mixed 


RATIO 


0.15 

0.43 

0.30 

0.14 


0.27 

0.07 


0.52 

0.20 


0.18 

0.56 

0.30 


0.12 

0.41 

0.13 


0.21 

0.45 

0.14 


0.14 

0.40 

0.10 


0.17 

0.35 

0.10 





















































Mariner 9 data are used as a benchmark by a!! photochemical models of 
the Mars atmosphere, this has important implications for much of the 
modeling work done over the last 20 years. 

Figure 6 shows that even simply ignoring a minimal dust opacity of 
0.3 in a cloud-free atmosphere will underestimate the ozone abundance 
by a factor of 2. This is noteworthy in that dust was completely ignored 
in the earlier inferences of ozone from the spectra by Lane et al. (1973) 
and others. Dust abundances are always of this order or higher (Pollack 
et al., 1979). 

Figure 7 shows that for high dust/cloud loading the inferred ozone 
abundance is a factor of 4 or 5 below what is truly there. Pollack et al. 
(1979) note that perhaps 20% of the martian year during the Viking 
observation period had this extent of dust loading or greater. 

Table 2 does show that errors in retrieval are the same for all 
solar zenith angles we looked at, and for all ice albedo below 0.5. For 
high ice albedo, ozone becomes even more severely underestimated. 

11.6. Reflectance Spectroscopy efficacy. 

Figures 3 through 7 and Table 2 show that to properly retrieve the 
ozone abundance from the Mariner 9 UV spectra with the reflectance 
spectroscopy technique, the effects of cloud and dust must be consid- 
ered. Therefore, the efficacy of the reflectance spectroscopy technique 
is only as good as the accuracy with which we know the opacity and 
scattering properties of the dust and cloud itself. Clancy and Lee (1991) 
have done the definitive study to date of retrieving these cloud and dust 
parameters. However, even with their study, there remain large uncer- 
tainties in the spatial and temporal variability in opacity and scattering 
properties, and also in the wavelength dependence of these properties. 



The amount of other spacecraft sensors and observations needed to 
obtain these properties to the degree needed to obtain ozone abundance 
to within a 10 to 20 % uncertainty is enormous. This raises serious 
doubts about the efficacy of this technique for retrieving ozone. Also, 
the reflectance spectroscopy technique does not convey any information 
about night-time ozone, and since most of the winter polar region where 
ozone is abundant is in night, most of the region where ozone is abundant 
is unexplored. 1 seriously recommend that other techniques be examined 
before designing future Mars Aeronomy sensors; perhaps detecting ozone 
at other wavelengths such as in the infrared, or by detecting ozone in- 
situ, as I have proposed to NASA. 

11-7. Presentations made at conferences. 

Several conferences were attended which were partly or fully 
funded by this contract. Reprints of papers and abstracts published at 
these meetings are included in the Appendix. 

IIJGG Assembly. I attended the International Union of Geodesy and 
Geophysics Assembly in Vienna, Austria, in August 11-24 1991, and pre- 
sented a paper entitled "Mars seasonal C02-ice lifetimes and the angular 
dependence of albedo" in the special Mars climate session. The abstract 
appeared in the conference proceedings (see Appendix). Half of expenses 
were paid by NASW-4614 and half were paid by another contract. 

PCI Symposium. I attended the International Symposium on the 
chemistry and physics of ice, held in Sapporo, Japan Sept. 1-7, 1991, and 
presented a paper entitled "Why is the north polar cap on Mars different 
than the south polar cap?" in the extraterrestrial ice session. The 
abstract appeared in the conference proceedings (see Appendix). 

Contract NASW-4614 paid for registration and publication costs only. 
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MSATT Workshop. I attended the Mars Surface and Atmosphere 
Through Time Workshop, held in Boulder, Colorado Sept. 23-25, 1991, and 
presented a paper entitled "Simulations of the seasonal polar caps on 
Mars". The abstract appeared in the workshop proceedings put out by the 
Lunar and Planetary Institute in Houston (see Appendix). NASW-4614 
paid for expenses, but these were limited mostly to airfare and regis- 
tration due to gratis accommodations, food, and car (paid by LPI). 

MSATT Workshop. I attended the Mars Surface and Atmosphere 
Through Time Workshop on innovative instrumentation for the in situ 
study of atmosphere-surface interactions on Mars held in Mainz, Germany 
October 8-9, 1992. Results presented in this report were presented in a 
paper entitled "Does UV instrumentation effectively measure ozone 
abundance?" (see Appendix). Grant NASW-4614 paid partial support 
(45%) (with remaining support coming from another contract). 

Lab. Research Planetary Atmospheres Conference . I attended the 
Fourth International Conference on Laboratory Research for Planetary 
Atmospheres held in Munich, Germany October 10-11, 1992. A paper 
entitled "How well is martian ozone inferred with reflectance spec- 
troscopy?" presented results from this contract (see Appendix). Grant 
NASW-4614 provided for partial support (note that this conference was 
concurrent with the MSATT Workshop, limiting costs). 

Planetary Science Conference . I attended the American 
Astronomical Society Division of Planetary Science Conference held in 
Munich, Germany October 12-16, 1992. A paper entitled "Martian polar 
cap seasonal regression simulations" was presented (see Appendix). 

Grant NASW-4614 provided for partial support (note that this conference 
was concurrent with the MSATT Workshop, limiting costs). 
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MSATT Workshop. I attended the Mars Surface and Atmosphere 
Through Time Workshop on The Polar Regions of Mars: Geology, 

Glaciology, and Climate History held in Houston, Texas November 13-15, 
1992. A paper entitled "Is carbon dioxide ice permanent?" was presented 
and published by LPI (see appendix). Grant NASW-4614 provided support. 

AMS Annual Meeting Chemistry Session . I attended the American 
Meteorological Society Annual Meeting Session on Atmospheric 
Chemistry held in Anaheim, California January 1993. A paper entitled 
"Atmospheric chemistry on Mars" presented results from this contract 
(see Appendix). The bulk of support came from another contract, but 
grant NASW-4614 paid for registration and publication costs for this 
session. 

AAAS Annual Meeting. I attended the American Association for the 
Advancement of Science annual meeting held in Boston, Mass., Feb. 1993. 
A paper presenting results from this contract entitled "The abundance of 
ozone on Mars" was presented and published by the AAAS (see appendix). 
Grant NASW-4614 provided support (limited to registration since I live 
in Boston). 

Optical Remote Sensing of the Atmosphere. I attended the Optical 
Remote Sensing of the Atmosphere conference held by the Optical 
Society of America in Salt Lake City, Utah, March, 1993 and presented a 
paper entitled "Probing the martian atmosphere in the ultraviolet" which 
presented results from this contract (see Appendix). Grant NASW-4614 
provided partial support, with the bulk of support supplied by another 
contract. 

MSATT Workshop. I attended the Final Mars Surface and 

Atmosphere Through Time Workshop on Mars: Past, Present, and Future - 
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Results from the MSATT Program held in Houston, Texas November, 1993 
and presented two papers entitled "The effect of polar caps on 
obliquity" and "How well was total ozone abundance inferred with 
Mariner 9?", which presented results from this contract. Both papers 
were published by the Lunar and Planetary Institute (see appendix). 

Grant NASW-4614 provided support. 

Low-Cost Planetary Missions Conference . An abstract entitled 
"In-Situ Mars Ozone Detector" has been accepted for presentation at the 
International Academy of Astronautics International Conference on Low- 
Cost Planetary Missions to be held in Baltimore, Maryland in April 1994. 
The abstract (copy enclosed in the appendix) will be published in a book 
of abstracts. 

M.8. Publications 

An invited contribution for the Questions and Answers section was 
published in The Planetary Report. Volume 14, p.21, 1994. This is an 
overview of ozone on Mars for the general public. A reprint is included 
in the Appendix. 

A paper entitled "Mars Ozone: Mariner 9 Revisited" has been 
submitted to Geophysical Research Letters . This paper summarizes the 
principal findings within this project. A preprint is included in the 
Appendix. 

A paper entitled "In-Situ Mars Ozone Detector" will soon be 
submitted to Acta Astronautica. This paper will describe a possible 
new ozone measurement technique which would avoid most of the error 
and uncertainty involved with the reflectance spectroscopy technique. 

Dr. Elliot Weinstock of the Harvard University Applied Physics Lab. has 
been part of a team which has measured ozone in the Earth’s 
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stratosphere via balloon, and we have been involved in examining the 
feasibility of sending this instrument to Mars. An abstract by the same 
name is in the appendix describing more details. A proposal by Dr. 
Weinstock and myself to undertake these studies is currently under 
consideration at the NASA PIDDP Program and the NASA Advanced 
Studies Program. 

A paper entitled "Ozone heating in the martian atmosphere" was 
published in Icarus. Volume 93, pp. 354-361, 1991. This paper 
describes the contribution of ozone to the energy budget of the martian 
atmosphere. This work was completed earlier, but would not have been 
published had I not received assistance under this contract. A reprint is 
included in the Appendix. 

A paper entitled "Sunlight penetration through the martian polar 
caps: Effects on the thermal and frost budgets" was published in 
Geophysical Research Letters. Vol. 19, pp. 1675-1678, 1992. This paper 
describes modeling studies of the polar cap frost and thermal budgets 
and the effect of including the penetration of sunlight into the caps, 
which was previously always assumed to be absorbed only at the surface. 
This work was completed earlier, but would not have been published had 
I not received assistance under this contract. A reprint is included in 
the Appendix. 

A manuscript entitled "The hemispherical asymmetry in the mar- 
tian polar caps" was published in J. Geophvs. Res. , Vol. 98, pp.3339- 
3344, 1993, in their special MSATT issue commemorating papers 
presented at the Workshop on the Martian Surface and Atmosphere 
Through Time. This paper describes how the radiative effects of clouds 
and dust result in a hemispherical asymmetry in the residual polar caps, 
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as has been observed. This work was completed earlier, but would not 
have been published had I not received assistance under this contract. A 
reprint is included in the Appendix. 

A chapter in the book Physics and Chemistry of Ice entitled "CO2- 
ice on Mars: Theoretical Simulations" has been published by Hokkaido 
University Press, Sapporo, N. Maeno and T. Hondoh, ed.s, 1992. This 
chapter describes how C02 ice stability can be modeled and what has 
been learned. This work was completed earlier, but would not have been 
published had I not received assistance under this contract. A reprint is 
included in the Appendix. 

A manuscript entitled "Cooling the martian atmosphere: The 
spectral overlap of the CO2 15 pm band and dust" has been accepted for 
publication in Planet. Space Sci. and is in Press. This paper describes 
the proper procedure for calculating the infrared radiative transfer in 
the martian atmosphere. This work was completed earlier, but would 
not have been published had I not received assistance under this 
contract. A preprint is included in the Appendix. 

A manuscript entitled "Martian atmospheric radiation budget" has 
been submitted to Publ. Astron. Soc. Japan . This paper describes the 
contribution of C02 and dust to the solar and IR energy budgets of the 
martian atmosphere. This work was completed earlier, but would not 
have been published had I not received assistance under this contract. A 
preprint is included in the Appendix. 

A paper entitled 'Comment on "Mars secular obliquity change due to 
the seasonal polar caps" by David Parry Rubincam' has been submitted to 
J. Geophvs. Res. . This paper describes the interaction between the 
orbital parameters of Mars and the polar caps of Mars. This work was 
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completed under this contract, but is really nothing more than a 
comment and hence does not merit any discussion here other than a 
reference to the preprint which is included in the Appendix. 

11.9. Honors Awarded. During the course of this research, Dr. 
Lindner received 2 honors for his Mars research. Dr. Lindner was elected 
as a new member of the International Association of Meteorology and 
Atmospheric Physics Commission on Planetary Atmospherics and their 
Evolution (a part of the International Union of Geodesy and Geophysics). 

Dr. Lindner was also recently nominated for the European 
Geophysical Society Young Scientists' Publication Award for his 
Planetary and Space Science Manuscript entitled "Cooling the martian 
atmosphere: The spectral overlap of the CO 2 15 pm band and dust". The 
winner of this award will be announced soon. 

III. Becommendations-and Conc isio ns 

The primary objectives for the work have been completed. Mariner 9 UV 
Spectrometer data have been reinverted for the ozone abundance. The spectra 
were fit by models which covered the full range in observed solar zenith angle, 
cloud, dust and ozone amount, ice albedo and look angles. This work has shown 
that significant underestimation of ozone occurred in earlier analysis of Mariner 9 
data, and that much of the observed variability in Mars ozone is due to masking of 
ozone by clouds and dust. As these previously published abundances have been 
used as a benchmark for all theoretical photochemical models of Mars, some of this 
work will have to be reinvestigated. The efficacy of the reflectance spectroscopy 
technique at retrieving the ozone abundance on Mars is questioned. The technique 
is incapable of retrieving ozone at night and during dust storms, and has as much as 
a factor of 3 error at other times. An in-situ measurement by balloon is 
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recommended as it is the only technique capable of accurately inferring the ozone 
abundance in all conditions. Studies of the effectiveness of such a mission are 
recommended. 


24 



Publications wholelv or partly under this contract 
(reprints are in the appendix) 

1. Refereed Papers 

Lindner, B. L., Ozone heating in the martian atmosphere, Icarus . 93 . 354- 
361, 1991a. 

Lindner, B. L., Sunlight penetration through the martian polar caps: 

Effects on the thermal and frost budgets, Geophvs. Res. Lett. . 1 9 . 
1675-1678, 1992a. 

Lindner, B. L., C02-ice on Mars: Theoretical simulations, in Physics and 
Chemistry of Ice . N. Maeno and T. Hondoh, ed.s, pp. 225-228, 
Hokkaido University Press, Sapporo, 1992c. 

Lindner, B. L., The hemispherical asymmetry in the martian polar caps, jL 
Geophvs. Res. . 98 . 3339-3344, 1 993a. 

Lindner, B. L., Questions and Answers, in The Planetary Report . 14 . 
21,1994a. 

Lindner, B. L., Cooling the martian atmosphere: The spectral overlap of 
the CO 2 15 pm band and dust, Planet. Space Sci. . in Press, 1994b. 
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Ozone is shown to be a minor, although nonnegligible, heat 
source in the martian atmosphere, in contradiction to earlier re- 
ports which had suggested a larger contribution to the atmospheric 
heat budget by ozone. It is further shown that the observed hemi- 
spherical asymmetry in ozone abundance is unlikely to be a signifi- 
cant factor in the observed hemispherical asymmetry in the perma- 
nent polar caps. Ozone heating may still be an important 
consideration in the formation and decay of clouds and snowfall, 
although not as important as previously suggested, o 1991 Academic 

Pre», Inc. 


INTRODUCTION 

Prabhakara and Hogan (1965) constructed the first and 
only photochemical-radiative transfer model of the Mar- 
tian atmosphere. Their model indicated that ozone had a 
significant influence on the temperature structure. How- 
ever, they ignored the odd hydrogen catalytic cycle and, 
thus, overestimated the ozone abundance by several or- 
ders of magnitude. Also, since the abundance and compo- 
sition of the atmosphere were poorly known prior to 
spacecraft exploration of Mars, many inaccurate assump- 
tions regarding the atmospheric properties were made. 

Ozone was first measured on Mars by the Mariner 7 
UV spectrometer (Barth and Hord 1971) and was mapped 
across the planet for half a Mars year by Mariner 9 (Barth 
et al. 1973). Based on these measurements, ozone was 
suggested to be an important contributor to the atmo- 
spheric temperature in the polar regions by Kuhn et 
al.( 1979). Radiative equilibrium temperatures computed 
by Kuhn et al. were up to 10 K higher when ozone heating 
was included. Kuhn et al. also suggested that ozone may 
be important as a heat source in determining precipitation 
and cloud formation, as temperatures at polar winter lati- 
tudes are very near the condensation temperature of C0 2 
and HjO. They further speculated that the observed hemi- 
spherical asymmetry in 0 3 (Barth et al. 1973) may explain 
the observed hemispherical asymmetry in the composi- 
tion of the residual polar caps (Kieffer et al. 1976, Kieffer 
1979). 


This study examines the absorption of solar radiation 
by ozone in the Martian atmosphere more closely. Both 
earlier studies neglected dust heating rates, which are 
included here. Furthermore, the possibility that ozone 
may have an effect on polar-cap ice abundance is ex- 
amined. 

MODELING APPROACH 

The radiative intensity as a function of altitude from the 
surface to 40 km is calculated for 57 solar wavelength 
intervals. The discrete ordinate method of Stamnes and 
Conklin (1984) is used to compute the scattering and ab- 
sorption of solar radiation by 0 3 , C0 2 , 0 2 , H 2 0 , H 2 , H0 2 , 
H 2 0 2 , and dust. The mixing ratios of C0 2 , 0 2 , and H 2 
are taken as 95.4, 0.13, and 0.002%, respectively (An- 
derson 1974, Liu and Donahue 1976, Owen et al. 1977, 
Krasnopolsky and Parshev 1979). Significant reduction in 
the mixing ratio of CO z over the winter polar region be- 
cause of the condensation of C0 2 is unlikely, based on 
meterological constraints (Hess 1979). The mixing ratio 
of 0 2 is based on observations (Trauger and Lunine 1983), 
but is also reproduced through photochemical modeling 
(e.g., Lindner and Jakosky 1985). The winter polar atmo- 
sphere appears to be saturated in H 2 0 at all altitudes, 
although evidence is sparse (Davies, 1979, Jakosky, 1985). 
Thus the vapor pressure there of H 2 0 as a function of 
temperature is given by the Clausius-CIapeyron equation 
for saturation pressure. HO z and H 2 0 2 abundances are 
taken from Lindner (1988). The existence of H 2 0 2 has not 
been confirmed by direct measurement, and thus remains 
speculative. However, it has been predicted by atmo- 
spheric modeling (e.g., Kong and McElroy 1977b). The 
altitude profiles of C0 2 , 0 2 , H 2 , and the pressure are 
determined from the hydrostatic equation. The surface 
pressure for 57°N latitude in late winter is taken to be 8 
mbar, based on Viking lander observations (Hess et al. 
1980) corrected for elevation (Lindal et al. 1979, Jakosky 
and Farmer 1982). 

Late winter (L $ = 343°) conditions for a latitude of 57°N 
were chosen for this study, as these were the conditions 
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when Mariner 9 observed the maximum 0 3 column abun- 
dance of 57 jxm-atm (Barth et al. 1973). (The aerocentric 
longitude of the Sun, L s , is the seasonal index; L s = 0°, 
90°, 1 80°, and 270° correspond to northern spring equinox , 
summer solstice, autumnal equinox, and winter solstice, 
respectively). The 0 3 mixing ratio used here is virtually 
constant between 10 and 30 km altitude. The surface loss 
of ozone (Kong and McElroy 1977b) causes a reduction 
in the mixing ratio at the surface of 30% (Lindner 1988). 
Ozone surface loss is not as active in the winter polar 
regions because of the cold temperatures and frost-cov- 
ered surface. Ozone decreases above 30 km because the 
primary production mechanism is a three-body reaction 
which becomes ineffective at high altitudes. The altitude 
dependence of 0 3 abundance has not been measured in the 
polar regions, but is predicted by photochemical models 
(Lindner 1988), and is therefore subject to the usual uncer- 
tainties of these models. 

The properties of dust are extrapolated from observa- 
tions at mid-latitudes, given the lack of observations at 
winter polar latitudes. Background (i.e., non-dust-storm) 
dust opacities over the Viking landers varied seasonally 
from 0.2 to 1.0 (Pollack et al. 1979), although year-to- 
year differences seem likely. The altitude profile of the 
background dust mixing ratio is approximated here with 
a Gaussian profile of the form exp( - z 2 /650), where z is 
altitude in kilometers. This results in background dust 
opacity confined mostly below 20 km altitude. The wave- 
length dependence of the dust opacity is taken from Toon 
et a/.( 1977). The single-scattering albedo of airborne dust 
as a function of wavelength is taken from Zurek (1978, 
1982), using a solar average of 0.86 (Pollack et al. 1979). 
The Haze-L phase function (Diermendjian 1969) is used 
to describe the scattering of radiation by martian dust. The 
Haze-L phase function provided the best fit to Mariner 9 
IR spectra (Toon et al. 1977). 

Small temperature inversions are common for winter 
polar latitudes. The temperature profile at 57°N latitude 
that is used in this study rises linearly with altitude from 
150 K at the surface to 170 K at 6 km altitude and then 
falls linearly with increasing altitude to 130 K at 40 km. 
The temperature profile at 70°N latitude is the same as 
that used at 57°N latitude except that it has a peak temper- 
ature of 155 K at 6 km. These temperature profiles are 
extrapolated from radio occultation observations made 
near winter polar latitudes during background dust condi- 
tions (Lindal et al. 1979) and agree with temperatures 
obtained by the Viking IRTM experiment (Kieffer 1979). 
Temperatures are markedly higher at the time of the global 
dust storms, but global dust storms are not considered 
here. The surface at 57°N and 70°N latitude in late winter 
is covered with old ice (Iwasaki et al. 1982), with a typical 
albedo of 0.5 (Kieffer et al. 1976). The surface is assumed 
to reflect according to Lambert’s law. 


TABLE I 

Absorption Cross-Section References 

C0 2 Shemansky 1972, Lewis and Carver 1983 
Oj Daumont et at. 1983, Freeman et at. 1984, Paur and Bass 
1985; Griggs 1968 
O, Ogawa 1971 

H,0 Thompson et al. 1963, Hudson 1971 

H ; Oj Lin et al. 1978, Molina and Molina 1981, Demore et at. 1985 
HO, Hochandel et at. 1972, Demore cl at. 1985 

H, Hudson 1971 


Absorption cross-sections (or) for C0 2 , Oj, 0 2 , H 2 0, 
H 2 0 2 , H0 2 , and H 2 are taken from the references shown 
in Table I, taking into account the temperature dependen- 
cies of or(C0 2 ) and <r(0 3 ). The temperature dependence 
causes o-(C0 2 ) to change as much as a factor of 3 at 
some wavelengths. The importance of the temperature 
dependence of cr(C0 2 ) was stressed in a study of martian 
photodissociation and heating rates by Parisot and Zuc- 
coni (1984). The U V heating rate of C0 2 and other constit- 
uents differs by as much as a factor of 2 or more depending 
on whether or not the temperature correction is included 
(Parisot and Zucconi 1984). However, Parisot and Zuc- 
coni (1984) did not use the improved measurements of the 
temperature dependence of o-(C0 2 ) by Lewis and Carver 
(1983). Using the Lewis and Carver measurements, 1 find 
an even greater change in UV heating rates due to the 
temperature dependence. Furthermore, the temperature 
dependence of o-(C0 2 ) is even more important to consider 
in the winter polar region where temperatures are the 
lowest on the planet (Parisot and Zucconi studied equato- 
rial latitudes). 

The only appreciable (>10%) temperature-dependent 
changes in 0 3 cross-sections are longward of 2700 A. Total 
0 3 heating rates are lowered by up to 30% when the 
temperature dependence is included. Temperature de- 
pendencies for the cross-sections of 0 2 , H z O, HO,, and 
H 2 0 2 are not available, although this is expected to be a 
small effect (Frederick and Hudson 1980). All absorption 
of solar radiation is assumed to go completely into instan- 
taneous and local heating. 

The cross-sections and solar intensities were treated as 
constant over 50-A intervals. These intervals were then 
summed from 1500 to 3200 A, with an added term for the 
Chappuis bands of ozone from 4000 to 7000 A, in 100-A 
intervals. At wavelengths less than 1500 A, cross-sections 
were not needed since C0 2 is optically thick and no contri- 
bution to the heating rate is made, even at 40 km altitude. 
Wavelength intervals of 1 gm are used out to 5 /xm to 
compute dust heating rates. The solar intensity is taken 
from Rottman (1981), corrected for the orbital eccentricity 
and orbital radius of Mars. 

Atmospheric properties are zonally averaged and as- 
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sumed azimuthally independent. To account for vertical 
inhomogeneity, the atmosphere is divided into twenty 
layers, each of which is characterized by one pressure 
and one temperature. Rodgers and Walshaw (1966) have 
shown that an altitude spacing of 0.2 scale height results 
in errors in heating and cooling rates of less than 0.03 
degrees per day, and so a grid spacing of 2 km was chosen. 

Heating rates are diurnally averaged as in Cogley and 
Borucki (1976), except that their derivation is repeated 
using the Chapman function instead of the secant func- 
tion, integrating over time numerically. Because diurnal 
variations are characteristically double-valued (one value 
for daytime, switching quickly to another nighttime 
value), a two-level parameterization is adopted here. The 
day/night ratio in the number of densities of 0 3 , H0 2 , 
and H 2 0 2 is taken from Shimazaki (1981), which shows 
agreement with day/night ratios from Krasitskii (1978) 
despite the large difference in latitude (equatorial versus 
polar). 

The value of the specific heat was taken to be 840 J kg ~ 1 
K" '. Values ranging from 750 to 900 have been used over 
the years by the various researchers studying the Mars 
atmosphere (see Lindner 1985). Further experimental re- 
search as to the value of the specific heat of the martian 
atmosphere is recommended in view of its importance 
in atmospheric heating rates. An exponentially averaged 
density over each layer was used in the calculation of 
heating rates to prevent an artificially sloping profile of 
the heating rate with altitude. Aerosols were not included 
in the density as the column abundance of airborne aero- 
sols for normal dust loadings is 5 x 10~ 4 g cm -2 (Pollack 
et al. 1979), five orders of magnitude below the column 
abundance of gas . It is also important to remember that the 
atmosphere is only 95% C0 2 and that the flux calculations 
based on absorption by C0 2 use 95% of the atmospheric 
number density, while the heating rate applies to the den- 
sity of all species combined. Ignoring this effect can 
change heating rates by 5%, and resultant temperatures 
by 6 K (Kondrat’ev et al. 1973). 

RESULTS 

Figure 1 compares the vertical optical depths of dust, 
Rayleigh scattering, C0 2 , and 0 3 as functions of wave- 
length. Wavelengths shortward of 1700 A are optically 
thick even at 40 km altitude, due to absorption by C0 2 in 
the upper atmosphere. Appreciable contribution to atmo- 
spheric heating rates does not occur until 1800 A, at which 
point the opacities of Rayleigh scattering and dust are 
both appreciable. Gaseous absorption shortward of 1900 
A is due primarily to C0 2 , with some contribution by 0 2 . 
Between 1950 and 2150 A, Rayleigh scattering and dust 
each have optical depths greater than that of C0 2 and 0 3 . 
Ozone absorption is appreciable between 2000 and 3000 
A, having vertical optical depths of almost 2 at 2500 A. 


Because solar zenith angles are very large in the winter 
polar region, the optical depth the solar radiation actually 
traverses is many times the vertical optical depth shown 
in Fig. 1. This means that very little ultraviolet radiation 
shortward of 3000 A actually reaches the surface. At visi- 
ble wavelengths, the optical depth is dominated by dust 
absorption and scattering. The only gaseous absorption is 
that of the ozone Chappuis bands. As these bands are far 
weaker than the ultraviolet bands, the optical depth is 
much smaller. 

Rayleigh scattering is an important consideration in the 
computation of U V heating rates. Including Rayleigh scat- 
tering causes an increase in C0 2 UV heating rates above 
20 km altitude, and a decrease in C0 2 UV heating rates 
below 20 km altitude. A 30% increase in C0 2 UV heating 
rates occurs at 40 km altitude, and a 50% decrease occurs 
at the surface (Lindner 1988). The Rayleigh scattering 
optical depth as a function of wavelength is given by 
Hansen and Travis (1974). The refractive indices of C0 2 , 
Ar, N 2 , CO, and 0 2 (the major gases in the martian atmo- 
sphere) are 1.000450, 1.000281, 1.000300, 1.000300, and 
1 .000272, respectively (Handbook of Chemistry and Phys- 
ics 1983). Upon evaluation, it is found that CO, provides 
well over 99% of the Rayleigh scattering optical depth, 
due not only to comprising 95% of the atmosphere, but 
more to the larger index of refraction for C0 2 . Rayleigh 
scattering is significant in the ultraviolet; at 2000 A, r R 
(Rayleigh scattering opacity) = 0.224; at 2500 A, r R = 
0.084; at 3000 A, r R = 0.038. Note that the Rayleigh 
scattering optical depth at 3000 A reported by Kuhn and 
Atreya (1979) of less than 0.01 is in error, most likely a 
typographical error since their value of 0.07 at 2500 A is 
good. 

Ozone heating. Figure 2 shows the total atmospheric 
heating due to ozone at 57°N latitude in late winter (L s = 
343°), when the maximum ozone abundance of 57 /xm-atm 
was observed by Mariner 9 (Barth et al. 1973). Note that 
this heating rate is directly related to the assumed altitude 
dependence of the 0 3 mixing ratio, and is therefore depen- 
dent on the accuracy of the photochemical model. Also 
shown in Fig. 2 is the contribution to the overall heating 
by a few select wavelength bins, which provides an idea of 
the wavelength dependence. Extreme ultraviolet heating 
(EUV) at wavelengths less than 1700 A is unimportant in 
the lower atmosphere of Mars due to absorption in the 
upper atmosphere. Ozone heating at wavelengths less 
than 1900 A is only important at altitudes above the sur- 
face. Large optical depths of C0 2 prevents solar radiation 
from reaching lower altitudes, and explain the strong alti- 
tude variation in the heating rate at these wavelengths. At 
2000 A, the atmospheric optical depth has a minimum (see 
Fig. 1). Solar flux at 2000 A is only weakly attenuated at 
the surface and the heating rate is virtually constant in 
altitude. 

Beyond 2000 A, absorption cross-sections of ozone in- 
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FIG. 1. Contribution of Oj, C0 2 , dust, and Rayleigh scattering to the vertical optical depth of the atmosphere at the surface at 57°N latitude 
and L , of 343° at UV wavelengths (a) and visible wavelengths (b). 


crease with wavelength to a maximum at 2500 A (see Fig. 
1). Increasing cross-sections and increasing solar fluxes 
result in larger heating rates, as shown in Fig. 2. Large 
optical depths attenuate the solar flux at lower altitudes, 
which results in a strong altitude dependence in the ozone 
heating rate at 2500 A. While heating increases with wave- 
length from 2000 to 2500 A at altitudes greater than 10 km 
above the surface, heating at the surface actually de- 
creases with increasing wavelength due to the attenuation 
of solar flux. From 2500 to 2700 A the absorption cross- 
sections of ozone decrease with increasing wavelength. 
However, the solar flux increases strongly with increasing 


wavelength from 2500 to 2700 A, resulting in a maximum 
ozone heating above 10 km at 2700 A. Beyond 2700 A, 0 3 
cross-sections drop dramatically, which results in de- 
creasing heating rates, but also allows more solar flux to 
reach the surface. As a result, the maximum ozone heating 
at the surface occurs at 2900 A. 

Ozone cross-sections reach a secondary maximum at 
visible wavelengths in the Chappuis bands. While not 
the major source of ozone heating, the Chappuis bands 
provide over 10% of the total ozone heating near the 
surface. As the atmosphere is optically thin at visible 
wavelengths (see Fig. 1), the heating is virtually constant 



OZONE HEATING RATE (K/DAY) 

FIG. 2. Total atmospheric heating rate due to O, at 57°N latitude and the contribution by a few 50-A wavelength intervals in the UV and by 
a 500-A wavelength interval in the visible, centered at the wavelengths shown (in A). 
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with atltitude. The maximum ozone heating in the Chap- 
puis bands occurs at 5200 A. 

The downward solar flux as a function of wavelength 
and altitude is shown in Fig. 3. Flux shortward of 2000 A 
is absorbed by CO,. Note the strong decrease in the solar 
flux striking the polar cap between 2200 and 2800 A due 
to ozone absorption. In the center of the Hartley bands at 
2500 A, the solar flux impinging on the polar cap is reduced 
by two orders of magnitude due to ozone absorption. 
Above 20 km altitude, the solar flux is not significantly 
attenuated due to the low ozone abundances at these 
altitudes (Lindner 1988). 

Significantly lower ozone abundances exist at 70°N lati- 
tude (Barth et al. 1973) and result in lower ozone heating 
rates. To emphasize latitudinal variations, the ratio of the 
0 3 heating rate at 70°N latitude to 0 3 heating rates at 
57°N latitude is shown in Fig. 4, along with a few select 
wavelength intervals. While ozone abundances are lower 
at 70°N latitude, the slant path is larger. Larger slant paths 
result in greater solar flux reductions near the surface, 
particularly where large CO, absorption (e.g., 1825 A) or 
large 0 3 absorption (e.g., 2525 A) exists. This causes 
larger reductions in heating rates near the surface at those 
wavelengths. As Rayleigh scattering optical depths are 
greater than CO, and 0 3 absorption optical depths at 2000 
A (see Fig. 1), the reduction in the heating rate is less at 
70°N latitude at 2000 A. At 3000 and 5000 A the altitude 
behavior is the same at both latitudes because the vertical 
optical depth of the atmosphere is small. 

Ozone heating relative to dust heating. Figure 5 com- 
pares ozone and dust heating rates for two latitudes, all 
summed over the solar spectrum. Dust heating is larger 
than 0 3 heating at all altitudes, particularly above 20 km. 


0 3 heating is at least a factor of 3 less than dust heating. 
Figure 5 also shows that 0 3 heating is less important 
relative to dust at 70°N latitude than it is at 57°N latitude, 
due to the lower 0 3 abundances at 70°N latitude. There- 
fore, while dust heating is less at higher latitudes due to 
the decreased solar flux, the importance of dust heating 
relative to 0 3 heating has increased. The relative impor- 
tance of dust heating continues to increase right up to the 
edge of the polar night. 0 3 heating will also not be im- 
portant relative to dust at equatorial and mid-latitudes 
because negligible ozone abundances exist there (Barth 
et al. 1973). 



FIG. 4. Ratio of the total Oj heating rate at 70°N latitude to that at 
57°N latitude. Also shown are the ratios of the Oj heating rate in a few 
50-A wavelength intervals. 
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FIG. S. Total atmospheric heating due to O s and dust (r v = 0.2) at 57°N latitude (A) and 70°N latitude (B). 


Ultraviolet C0 2 , 0 2 , H 2 0, H0 2 , and H 2 0 2 heating. 
C0 2 , 0 2 , HjO, H0 2 , and H 2 0 2 also absorb in the UV. 
Integrating the heating from 1500 to 3200 A gives the total 
heating rates for each molecule, shown in Fig. 6 as the 
ratio of the heating to the total 0 3 heating. H 2 0, H0 2 , and 
H 2 0 2 number densities are too low for any appreciable 
UV heating. However, H 2 0 UV heating rates are much 
higher at polar summer, equatorial, or mid-latitudes where 
H 2 0 abundances are appreciable. C0 2 and 0 2 U V heating 
rates are appreciable in the upper atmosphere, being in 
fact the dominant source of atmospheric heating above 40 
km altitude, but are ineffective at low altitudes due to 
optically thick absorption bands. However, C0 2 and 0 2 
U V heating rates are also higher at polar summer, equato- 



FIG. 6. Ratio of the UV heating rates of C0 2 , 0 2 , H 2 0, HO,, and 
H 2 0 2 to those of 03 at 57°N latitude. 


rial, or mid-latitudes, due to lower solar zenith angles 
which result in less optically thick absorption paths. 

DISCUSSION 

Prabhakara and Hogan (1965) and Kuhn et a/. (1979) 
have shown that 0 3 UV heating rates are comparable to 
C0 2 near-IR heating rates. However, neither work consid- 
ered dust heating, which is at least a factor of 3 larger 
than 0 3 heating rates at all altitudes and latitudes. 0 3 
abundances are usually much smaller than the maximum 
observed, which will further reduce the importance of 0 3 
heating to the atmospheric thermal budget. Moreover, 
the UV measurements of Martian 0 3 cannot distinguish 
vapor-phase 0 3 from 0 3 adsorbed in the polar cap (Lane 
et al. 1973), and therefore the atmospheric burden may 
be less than the assumed 0 3 abundance. Furthermore, 0.2 
vertical optical depths of dust as used in Fig. 5 represent 
the minimum amount of dust observed (Pollack et al. 
1979). Larger dust loading will provide proportionally 
greater dust heating (Lindner 1985). Therefore, 0 3 is a 
minor, although nonnegligible, heat source in the Martian 
atmosphere. UV heating from CO,, 0 2 , H,0, H0 2 , and 
H 2 0 2 is also a minor heat source in the lower winter polar 
atmosphere. 

Mariner 7 and 9 observations indicate almost twice as 
much 0 3 over the northern polar cap as over the southern 
polar cap (Barth et al. 1973, Lane et al. 1973), due primar- 
ily to the hemispherical differences in airmass (Kong and 
McElroy 1977a, Lindner 1988). Kuhn et al. (1979) sug- 
gested that the observed hemispherical asymmetry in the 
composition of the permanent polar caps (Kieffer et al. 
1976, Kieffer 1979) could be related to this hemispherical 
asymmetry in 0 3 . 0 3 affects surface frost formation by 






360 


BERNHARD LEE LINDNER 


changing the atmospheric temperature, thus changing the 
heat received by the surface from atmospheric thermal 
radiation and from atmospheric heat transport. However, 
ozone has been shown here to be a minor solar heating 
source, compared to dust, and has been also shown to be 
a minor thermal radiator on Mars by Crisp (1990). Hence, 
the hemispherical asymmetry in 0 3 should not cause any 
significant hemispherical asymmetry in atmospheric tem- 
perature. 0-, also affects surface frost formation through 
absorption of solar radiation which would otherwise strike 
the polar cap. However, dust absorbs or scatters to space 
most UV light before it strikes the polar cap (Lindner 
1990). Furthermore, while 0 3 strongly absorbs at UV 
wavelengths as shown in Fig. 3, 0 3 absorbs less than 1% 
of the total solar flux which strikes the polar cap. Hence, 
the hemispherical asymmetry in 0 3 should not cause any 
significant hemispherical asymmetry in the total solar flux 
absorbed by surface frost. Since atmospheric tempera- 
tures are near the condensation temperature of C0 2 and 
H 2 0 at winter polar latitudes, 0 3 may still be important in 
the formation of clouds and snow, as Kuhn et al. (1979) 
suggested. The importance of 0 3 in cloud and snow forma- 
tion is less than Kuhn et al. predict since they ignored 
meridional transport (Pollack et al. 1990) and dust heating, 
but an 0 3 heating rate of 1 K/day will still significantly 
reduce atmospheric condensation. However, while the 
ratio of CO ; frost to C0 2 snow in the polar cap is difficult 
to determine. Pollack et al. (1990) have predicted that 
surface frost formation dominates over snowfall. There- 
fore, the hemispherical asymmetry in 0 3 is unlikely to 
play an important role in explaining the hemispherical 
asymmetry in the permanent polar caps. 
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Abstract. An energy balance model of the seasonal polar 
caps on Mars is modified to include penetration of solar radia- 
tion into and through the ice. Penetration of solar radiation has 
no effect on subsurface temperature or total frost sublimation if 
seasonal ice overlies a dust surface. An effect is noted for sea- 
sonal ice which overlies the residual polar caps. For the case 
of an exposed water-ice residual polar cap, the temperature at 
depth is calculated to be up to several degrees wanner and the 
calculated lifetime of seasonal C02 frost is slightly lower 
when penetration of sunlight is properly treated in the model. 
For the case of a residual polar cap which is perennially cov- 
ered by CX >2 frost, the calculated lifetime of seasonal CO 2 
frost is very slightly increased as a result of sunlight penetra- 
tion through the ice. Hence, penetration of sunlight into the 
ice helps to stabilize the observed dichotomy in the residual 
polar caps on Mars, although it is a small effect. 

Introduction 

Computer simulation of the condensation and sublimation 
of CO 2 frost in the martian polar caps has been fairly success- 
ful in reproducing the annual cycle in atmospheric pressure 
observed by the Viking Landers [Leighton and Murray, 1966; 
Cross, 1971; Briggs, 1974; Davies et al., 1977; James and 
North, 1982; Lindner, 1985, 1986]. However, these studies 
have not been able to uniquely explain why CO 2 frost survives 
summer in the southern hemisphere and not summer in the 
northern hemisphere. Several theories have been advanced to 
explain this discrepancy, including hemispherical asymmetries 
in albedo [Paige, 1985; Paige and Ingersoll, 1985; Lindner, 
1985; 1986], in snowfall [Pollack et al., 1990], in the radiative 
effects of clouds and airborne dust [Briggs, 1974; James and 
North, 1982; Lindner, 1990; 1992], in subsurface heat 
conduction [Jakosky and Haberle, 1990], and in ozone [Kuhn 
et al., 1979; Lindner, 1988; 1990; 1991]. It is also possible 
that a permanent reservoir of CO 2 frost from an earlier epoch 
is being uncovered in the southern hemisphere, although this 
seems unlikely [e.g., Jakosky and Barker, 1984], 

Prior models have used the assumption that all non-re- 
flected sunlight is absorbed at the surface, when in reality sun- 
light penetrates into the surface, sometimes to several meters 
depth [Clow, 1987]. This paper describes studies of how this 
phenomenon affects the thermal and frost budgets of the polar 
caps and the subsurface. 
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Modeling Approach 

To perform this study, the energy-balance model of 
Jakosky and Haberle [1990] is combined with the ice micro- 
physics model of Clow [1987]. The energy budget at the sur- 
face of the geographic poles involves balance between ab- 
sorbed sunlight, absorbed and emitted thermal-IR radiation 
(from the surface to space and from the atmosphere.to the sur- 
face), conduction to or from the subsurface (the subsurface is 
taken to begin at the base of the seasonal ice), and condensa- 
tion or sublimation of C02- The instantaneous energy balance 
at the surface can be written as 


So 6T(z,t)| 

p0-A)P(t)cos(i) - eoT^tet) - K ' 


dm(t) 

*=o* L -5r*° (1> 


where So is the solar constant at 1 AU corrected for atmo- 
spheric absorption and scattering [Lindner, 1985; 1990; 
Lindner et al., 1990]; R is the distance of Mars from the Sun 
in AU; A is the bolometric Bond albedo of the surface material 
at the pole; P(t) is the fraction of non-reflected sunlight which 
is absorbed by the seasonal polar cap at time t (the remainder 
being absorbed by the surface underneath the seasonal ice); 1 is 
the solar incidence angle; e is the effective emissivity of the 
surface; o is the Stefan-Boltzmann radiation constant; T(z,t) is 
the temperature at depth z (surface temperature is set to the 
condensation temperature of C02 if CO 2 frost is present); K is 
the thermal conductivity of the surface and subsurface materi- 
als; L is the latent heat of sublimation of CO 2 frost; dm(t)/dt is 
the time derivative of the CC>2 surface frost abundance (set to 
zero when no frost is present); and the temperature gradient 
8T/5z is evaluated at the surface (z=0). 

Sunlight penetration and heat conduction into the subsur- 
face are accounted for by 


5T(z,t) 1 5^T(z,t) So SF(z,t) 

— ‘ Tc' K + ^ “ S< ° "li - 


( 2 ) 


where p is the bulk density of the subsurface material, C is its 
specific heat, and 8F(z,t) is the fraction of the sunlight which 
passes through the seasonal cap and is absorbed within depth 
5z of the subsurface. The derivatives in (2) apply locally. The 
values used for p and C of the subsurface dust are 0.93 gcm‘3 
and 0.15 cal g"‘K'l, respectively, and K is allowed to vary 
depending on the thermal inertia. The albedo and thermal in- 
ertia for dust materia] is assumed to be comparable to the typi- 
cal regolith values (0.25 and 0.006 cal cm'^s'^K"^, respec- 
tively). The albedo and thermal inertia for a residual polar cap 
are taken to be like those of the north residual polar cap (0.45 
and 0.03 cal curbs' respectively [see Kieffer et al., 
1976; Paige, 1985]. When the model uses an albedo of 0.65 
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for seasonal CO 2 ice, all seasonal ice sublimes during the 
summer and the underlying surface becomes exposed. Hence, 
an albedo of 0.65 is used to simulate the case of an exposed 
residual polar cap. An albedo of 0.74 is needed in this model 
to retain seasonal CO 2 ice through the summer. While energy 
is deposited at depth by solar radiation, it is not effectively 
transported by IR radiation relative to conduction. Optical 
absorption coefficients are large in the IR, and temperatures 
are low, both of which inhibit transport by radiation. Hence, 
IR radiation is not included in (2). Terms in (1) and (2) other 
than P and F are evaluated as in Jakosky and Haberle [1990]. 

Clow [1987] calculated the multiple scattering of sunlight 
by the matrix of water-ice grains and dust particles on Mars by 
inserting the average single-scattering parameters for these 
particles into the 5-Eddington method of solving the radiative- 
transfer equation. Figure 1 shows the integral of net short- 
wave flux with respect to depth, which Clow calculated for 
pure snow and dirty snow (1000 ppmw dust) and for fine- 
grained snow (icc-grain radii of 50 Jim, bulk density of 50 kg 
m'3) and coarse-grained snow (ice-grain radii of 1000 jim, 
bulk density of 400 kg m'3). Values for P, the fraction of 
non-reflected solar radiation absorbed by the seasonal polar 
cap, are extracted from Figure 1 using the thickness of the sea- 
sonal polar cap as predicted by the energy balance model. 
Values for F(z), the fraction of non-reflected solar radiation 
absorbed by the surface under the seasonal ice, are then com- 
puted by taking the remainder (1-P) of the non-reflected solar 
radiation and apportioning it with depth as shown in Figure 1. 



Fig. 1. The integral of net shortwave flux with respect to 
depth for clean and dusty water snow on Mars, using the 
mean-annual incident solar flux at latitude -38°. Curves (A) 
correspond to a snow with ice grain radii of 50 pm and bulk 
density 50 kg m*3 while curves (B) are for a snow with 1000- 
jtm ice grains and a density of 400 kg m‘3. (Reproduced by 
permission from Clow [1987]). 


Kicffer [1990] and Moore [1988] have both deduced that 
ice-grain radii in the rcsiduai polar cap are over 100 pm and 
that dust concentrations are less than 1/1000. This means that 
ice in the residual polar cap is believed to be between fine- 
grained and coarse-grained and between pure and dirty, as I 
have defined them. There is no evidence for the grain size and 
dust content of seasonal ice, although it is likely to be less 
dirty and more fine grained than ice in the residual polar cap. I 
have assumed that the residual polar cap is at least several me- 
ters thick, which is believed to be the case [Toon et al, 1980], 

The Clow [1987] calculations presented in Figure 1 as- 
sume water ice, perfectly valid for the residual polar cap com- 
posed of water ice. Recent calculations by Warren et al. 
[1990] suggest that CO 2 ice on Mars is more absorbing than 
water ice on Mars. Greater absorption of sunlight by ice 
would decrease the penetration depth of sunlight shown in 
Figure 1 (G. Clow, personal communication, 1992). Hence, 
the penetration of sunlight calculated here for seasonal ice may 
be an upper limit case, which enhances the conclusions as 
shown below. 

To apply the model, (2) is numerically integrated at the ge- 
ographic pole through a martian orbit around the Sun, subject 
to the boundary condition of (1). The model is iterated for six 
Mars years, by which time convergence is achieved. The sub- 
surface is divided up into discrete layers, each of which is de- 
scribed by a center temperature. A 5 cm surface layer thick- 
ness is used. This is required since the layering needs to be on 
a finer scale than the sunlight penetration depth (see Figure 1) 
and the thermal skin depth. The thermal skin depth varies 
from 1 m to 10 m at the pole, depending on the choice of sub- 
surface material. Succeeding layers are thicker by a factor of 
1.2 each; this reaches to a depth of 37 m (over 4 annual skin 
depths). This subsurface thermal model produced excellent 
agreement to the depth profiles of temperature computed by 



Fig. 2. Temperature versus depth in a water-ice residual polar 
cap for the case where seasonal ice sublimes completely in the 
summer. The surface is still covered by CO 2 ice at this time 
(late summer, Ls = 260°). The dashed curve is simulated by 
the model including the effect of penetration of solar radiation, 
and the solid curve is simulated by the model with no penetra- 
tion of solar radiation. Surface values arc shown at 1 cm 
depth. Seasonal ice albedo and thermal inertia is 0.65 and 
0.03 cal cm‘2 g-lfe K*l, respectively. The flux profile for 
clean, fine-grained ice in Figure 1 is used for seasonal ice and 
that for clean, coarse ice is used for the residual polar cap. 
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Kieffer f 1990] for the same case he presented in his Figure 3. 

Results 

In the case where seasonal ice overlies a dust surface, our 
model calculates no changes in either the frost budget of the 
polar cap or the thermal structure of the subsurface due to sun- 
light penetration. The top grains of dust under the ice absorb 
the solar radiation which passes through the ice, and it is much 
easier for these top grains to conduct or radiate heat back to the 
overlying ice than deep into the subsurface. Effectively, all 
non-reflected sunlight is absorbed by the seasonal frost, as 
was assumed in previous polar-cap models. 

The behavior for seasonal C02 ice overlying a residual 
polar cap is different Solar radiation which passes through 
the seasonal ice may penetrate quite deeply into the residual 
polar cap (see Figure 1). For a water-ice residual polar cap, 
the sunlight which penetrates the seasonal polar cap heats the 
residual polar cap by up to 3°K, primarily in late summer when 
the seasonal ice is thinnest, with the greatest heating occurring 
at 1 m depth in the residual polar cap (see Figure 2). Most of 
the penetrating radiation actually gets absorbed near the surface 
of the residual polar cap, but that easily conducts or radiates 
away to the seasonal CO2 ice, accounting for the low heating 
rate at the surface of the residual polar cap. Dirty ice exhibits 
the same behavior seen in Figure 2, but with only half the 
magnitude of subsurface heating. The results are essentially 
the same no matter whether fine-grained or coarse-grained ice 
is used. However, the heating may be less than shown in 
Figure 2 since the seasonal ice may be more absorbing in the 
visible than assumed here. Different values for thermal inertia 
also change the degree of subsurface heating, with lower 
thermal inertia producing greater heating. For a residual polar 
cap which contains CO2 ice, any subsurface heating goes into 
subliming local CO2 ice, and does not change the local tem- 
perature until all local CO2 ice has sublimed. 

For the case of an exposed water-ice residual polar cap, the 
maximum surface temperature at the pole occurs just after the 
last of the seasonal ice sublimes. This is also when the maxi- 
mum change due to sunlight penetration occurs in the calcu- 
lated temperature at depth, as shown in Figure 3. I find that 
surface temperatures are cooler in early summer (L$ = 270° to 
Ls = 310°) for the model with light penetration than for the 
model which has no penetration. Solar radiation which pene- 
trates into the subsurface must conduct to the surface to heat 
the surface, which is not as effective at heating the surface as 
direct absorption of all the solar radiation by the surface. In 
late summer (Ls = 310° to 360°), the surface is actually warmer 
for the model with light penetration, due to the increased heat 
conduction by the warmer subsurface. After L s = 360°, the 
surface becomes covered with seasonal CO2 frost. Again, the 
maximum subsurface heating occurs at 1 m depth in a water- 
ice residual polar cap. 

The model which includes penetration of sunlight predicts 
less seasonal ice all year for the exposed water-ice residual cap 
case, as shown in Figure 4. By allowing solar radiation to be 
absorbed at depth, the surface does not become as warm in 
early summer when surface temperatures are their warmest. 
This decreases the annual-total energy lost as infrared radiation 
emitted by the surface and increases the amount of energy 
stored in the subsurface, which in turn decreases the amount 



Fig. 3. As in Figure 2, except after all seasonal ice has sub- 
limed in the model (Ls = 280°, 37 Mars days after Figure 2). 

of frost condensation needed in the early winter to maintain 
balance (see equation 1). Thus, there is an increasing differ- 
ence in predicted frost abundances between the models from 
Ls = 0° to 180°. After Ls = 180°, the difference in the (X>2 
frost amount between the models decreases because some 
sunlight penetrates the seasonal ice to heat the residual polar 
cap instead of subliming seasonal ice. However, this effect is 
not enough to compensate for the increased conduction of heat 
during early winter, and the net result is that the seasonal ice 
sublimes away earlier in the year for the model which includes 
penetration of sunlight Furthermore, seasonal ice may be 
more absorbing than assumed here, which will decrease heat- 
ing Of the residual polar cap while seasonal ice is present, 
resulting in a slightly larger difference between models in 
Figure 4. The effect of light penetration on the frost budget 
for dirty ice shows the same behavior as in Figure 4, but only 
half the magnitude. 

Discussion 

The inclusion of light penetration in an energy balance 
model slightly decreases the albedo needed in the model to 



Fig. 4. Annual variation in seasonal CO2 frost amount at the 
south pole for the case where seasonal ice sublimes completely 
in the summer. The difference between the two models is also 
plotted, using the scale on the right side. The same ice prop- 
erties as described for Figure 2 are used. 
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keep seasonal CO2 ice on a residual polar cap through the 
summer. Furthermore, the number of days that an energy-bal- 
ance model calculates for which an exposed water-ice residual 
polar cap is free of seasonal CO2 ice is increased by approxi- 
mately 5 days when light penetration is included in the model. 
Thus, since the effect of light penetration slightly decreases 
CO2 ice lifetimes when the ice does not exist year-round and 
slightly increases CO2 ice survivability when ice does exist 
year-round, it also enhances the Jakosky and Haberle (1990) 
conclusion that there are two stable states for the residual polar 
caps; perennially covered by CO2 ice or exposed every sum- 
mer. These results also enhance the conclusion of Jakosky 
and Haberle that conduction of heat to and from the subsurface 
plays an important role in the energy balance of the polar cap. 

Currently CO2 ice does not survive summer on the north- 
ern residual polar cap while CO2 ice does survive summer on 
the southern residual polar cap. The CO2 ice at the south pole 
could have originated from an earlier epoch, but that would 
make the large abundances of water vapor observed in south- 
ern summer in 1969 harder to explain [Jakosky and Barker, 
1984]. The inclusion of light penetration in an energy balance 
model makes it slightly easier for energy balance models to 
maintain seasonal CO2 ice on an existing CO2 residual polar 
cap, as now exists in the south, and to totally sublimate sea- 
sonal CO2 ice on an exposed water ice residual polar cap, as 
now exists in the north. Hence, the effect of sunlight penetra- 
tion makes it easier for energy balance models to explain the 
dichotomy in the residual polar caps. 

While of some importance directly at the poles, the pene- 
tration of light into the polar cap has only a small effect on the 
globally-integrated energy budget of the polar cap, specifically 
the globally-integrated CO2 sublimation and condensation as 
inferred by the Viking Lander measurements of atmospheric 
pressure and as predicted by theoretical models of the general 
atmospheric circulation. Given the uncertainties currently pre- 
sent in albedo and other parameters, the effect of light penetra- 
tion is second order, and can be currently neglected in models 
of the globally-integrated energy budget of the polar caps. 
Further work remains to be done on other processes involved 
in the frost budget of the polar caps, as discussed in the intro- 
duction. These processes would not affect the major results of 
our work, but may have a significant impact on the globally- 
integrated frost budgets. 
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ABSTRACT: A theoretical model of the energy budget of the polar caps of 
Mars has been created which is used to study the hemispherical asymmetry 
in CO 2 ice. The observations which show survival of seasonal CO 2 ice in 
the southern hemisphere in summer and not in the northern hemisphere in 
summer have been reproduced. 


1. INTRODUCTION 

One of the most puzzling mysteries 
about the planet Mars is the hemi- 
spherical asymmetry in the polar 
caps. Every spring the seasonal 
polar cap of CO 2 recedes until the 
end of summer, when only a small 
part, the residual polar cap, re- 
mains. During the year that Vik- 
ing observed Mars, the residual 
polar cap was composed of water 
ice in the northern hemisphere (1) 
but was primarily carbon dioxide 
ice in the southern hemisphere 
(2). Scientists have sought to 
explain this asymmetry by modeling 
observations of the latitudinal 
recession of the polar cap and 
seasonal variations in atmospheric 
pressure (since the seasonal polar 
caps are primarily frozen atmo- 
sphere, they are directly related 
to changes in atmospheric mass). 
These models reproduce most as- 
pects of the observed annual vari- 
ation in atmospheric pressure 
fairly accurately. Furthermore, 
the predicted latitudinal reces- 
sion of the northern polar cap in 
the spring agrees well with obser- 
vations, including the fact that 
the CO 2 ice is predicted to com- 
pletely sublime away. However, 
these models all predict that Che 
carbon dioxide ice will also sub- 
lime away during the summer in the 
southern hemisphere, unlike what 


is observed, as is shown in Figure 

1. This paper will show how the 
radiative effects of ozone, 
clouds, and airborne dust, light 
penetration into and through the 
polar cap, and the dependence of 
albedo on solar zenith angle af- 
fect CC >2 ice formation and subli- 
mation, and how they help explain 
the hemispherical asymmetry in the| 
residual polar caps. These ef- 
fects have not been studied with 
prior polar cap models. 

2. MODEL DESCRIPTION 

The energy budget of the surface 
of Mars has been studied with a 
model which includes all of the 
processes shown schematically in 
Figure 2. The sources and sinks 
of energy for a square centimeter 
of surface include solar insola- 
tion which strikes the surface, 
modified for the absorption due to 
clouds and aerosols; infrared 
emission by the clear atmosphere 
and by clouds and aerosols to the 
surface; infrared emission by the 
surface to space; penetration of 
solar radiation into the surface; 
atmospheric heat transport as re- 
presented by a thermal wind; heat 
conduction in the subsurface; and 
latent heat of condensation of 
CC >2 . The net gain or loss of 
energy integrated over one martian 
day is used to compute either a 



change in the surface temperature 
or a change in the amount of CO 2 
frost present. Details on the mo- 
del are presented elsewhere 
(9,10,11). 
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Since O 3 is more prevalent in the 
northern hemisphere than in the 
southern hemisphere, O 3 was sug- 
gested as a cause for the hemi- 
spherical asymmetry in the residu- 
al polar caps (12). However, 03 
has since been shown to have a mi- 
nor effect on the atmospheric tem- 
perature (13), and hence on the 
infrared radiation which strikes 
the polar cap, and it has been 
shown that O 3 absorbs less than 1 % 
of the total solar radiation ab- 
sorbed by the polar cap (9) . 
Thus, O 3 is not an important con- 
sideration in the polar cap energy 
budge t . 
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The solar and thermal flux strik- 
ing the polar cap of Mars has been 
computed for various dust and 
cloud abundances and for three so- 
lar zenith angles (9). These cal- 
culations have been inserted in 
earlier versions of polar-cap mo- 
dels (10,11). Vertical optical 


Fig. 1. The seasonal recession of 
the south polar cap as observed 
over the last 20 years (3) and as 
predicted by (4,5,6, 7, 8 ). (The 
aerocentric longitude of the sun, 


L , is the seasonal index; L s 
. 5 ~ . a a « t n n n o 


0°, 90° , 180° and 270° correspond 
to northern spring equinox, summer 
solstice, autumnal equinox and 
winter solstice, respectively.) 
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depths of dust and cloud ranging 
from zero to 1 increase the ab- 
sorbed flux significantly in polar 
night, where the pole spends half 
of the year, as shown in Table 1 
(9). Observed hemispherical asym- 
metries in dust abundance, cloud 
cover, and surface pressure com- 
bine to cause a significant hemi- 
spherical asymmetry in the total 
flux absorbed by the residual po- 
lar caps (9) , which helps to ex- 
plain the dichotomy in the residu- 
al polar caps on Mars. 

Penetration of solar radiation in- 
to the cap itself is included, 
based on theoretical work (14). 
The inclusion of light penetration 
slightly decreases the albedo 
needed in the model to keep CC^- 
ice year-round at the south pole 
by on the order of 1%. The re- 
quired albedo is decreased because 
some solar radiation is used to 
heat the subsurface, and not all 
of this heat is transported back 
to the surface. Overall, penetra- 
tion of light into the polar cap 
has only a small effect on the po- 
larcap energy budget. 

Calculations of the dependence of 
the albedo of the martian polar 
caps on solar zenith angle (15) 
have also been included in the mo- 
del. Since the albedo of ice in- 
creases and becomes more forward 
scattering at higher solar zenith 
angles, and since the solar zenith 
angle becomes higher as one ap- 
proaches the pole, the albedo is 
greatest at the pole. This de- 
creases absorption of sunlight, 
hence increasing survivability of 
CC >2 ice. In fact, this increases 
the survivability of ice enough to 
offset the decrease in survivabi- 
lity of ice due to the radiative 
effects of clouds and dust. 

3. DISCUSSION 

The combination of the effects of 
solar zenith angle on albedo and 
the radiative effects of clouds 
and dust act to extend the life- 
time of CO 2 ice on the south pole 
relatively more than on the north 
pole, explaining the hemispherical 
asymmetry in the residual polar 
caps without the need of a hemi- 


TABLE 1. Flux Striking the Surrace ol 
the Polar Cap in Polar Night (90°N 
Latitude, L g = 343°) 


Wavelength 
Interval, pm 


Dust Opacity 

0.0 

0.2 

0.5 

5.4-10.0 

0.0 

0.7 

1.3 

10.0-21.0 

14.7 

29.4 

43.4 

21.0-99.9 

0.0 

13.5 

29.9 

Absorbed flux 

14.7 

43.6 

74.7 

The vertical 

optical 

depth of 

dus t is 

given. -Flux 

values . 

are given 

in units 

of J cm* sol 

- 1 




spherical asymmetry in polar cap 
albedo. Another positive aspect 
of this solution is that neither 
the inclusion of solar zenith an- 
gle effects on ice albedo nor the 
radiative effects of clouds and 
dust should appreciably change mo- 
del predictions of the annual cy- 
cle of pressure or polar cap re- 
cession equatorward of 75° lati- 
tude, since approximately 90% of 
the seasonal CO 2 frost is equator- 
ward of 80° latitude. Hence, the 
good model agreement noted by pri- 
or researchers to the seasonal cy- 
cle in atmospheric pressure and to 
the recession of the polar cap 
equatorward of 80° latitude is re- 
tained . 
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An energy balance model is used to study the behavior of CX>2 ice on Mars. The effect of the solar 
zenith angle dependence of albedo is to lengthen CO 2 ice lifetimes at the poles. Hemispherical asymmetries 
in cloud and dust abundance could result in the survival of seasonal CO 2 ice through summer in the south and 
not in the north, in agreement with observations. CO 2 ice observed in the summertime polar cap in the south 
could be of recent origin, although a permanent CO 2 polar cap cannot be ruled out. 


Introduction 

During the year that Viking observed Mars, the summer- 
time polar cap consisted of water ice in the northern hemi- 
-sphere [Kieffer et al„ 1976], but was primarily carbon dioxide 
ice in the southern hemisphere [Kieffer, 1979]. CO 2 ice was 
also inferred to be present in the summertime south polar cap 
during the Mariner 9 encounter [Paige el al., 1990]. Open ques- 
tions remain as to why a hemispherical asymmetry exists, and 
whether the summertime CO 2 ice in the southern hemisphere is 
primarily seasonal in origin, or whether a reservoir of CO 2 
from an earlier epoch is being uncovered. A permanent 
exposed CO 2 ice cap would be one source of CO 2 for possible 
past climates on Mars with dense CO 2 atmospheres [e.g., 
Sagan et al., 1973; Toon et al., 1980; Hoffert et al., 1981; 
Fanale et al., 1982; Francois et al., 1990], However, Earth- 
based measurements made in 1969 of water vapor in the 
martian atmosphere suggest that all CO 2 ice sublimed from the 
southern polar cap and exposed underlying water ice [Jakosky 
and Barker, 1984], which would be inconsistent with the exis- 
tence of a permanent CO 2 ice reservoir. 

Theoretical models of the energy budget of the surface have 
been created which simulate the formation and dissipation of 
carbon dioxide ice and which have been compared to observa- 
tions of polar cap regression and the annual cycle in atmo- 
spheric pressure [Leighton and Murray, 1966; Cross, 1971; 
Briggs, 1974; Davies et al., 1977; Narumi, 1980; James and 
North, 1982; Lindner, 1985; Moreau et al., 1991; Wood and 
Paige, 1992]. As shown in Figure 1, these models have been 
unable to preserve seasonal CO 2 ice at the south pole and still 
obtain agreement with observations of the polar cap regres- 
sion and the annual cycle in atmospheric pressure. No combi- 
nation of best-fit frost albedos and emissivities was consistent 
with the stability of permanent C02 deposits at either pole 
[Wood and Paige, 1992], This implies that either these models 
improperly treated the energy budget or that CC >2 ice from an 
earlier time is exposed. 

Clouds have been observed to be more prevalent over the 
edge of the northern seasonal polar cap than over the edge of 
the southern seasonal polar cap, although observations of 
clouds are far from complete [James, 1983; Christensen and 
Zurek, 1983; James et al., 1987]. Atmospheric dust experi- 
ences an annual cycle, with order of magnitude variations in 
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dust opacity [Pollack et al., 1979]. Surface pressure is lower 
during southern winter over the south pole than during north- 
ern winter over the north pole, due to differences in elevation 
[Woiceshyn, 1974; Lindal et al., 1979] and due to the annua] 
cycle in atmospheric pressure [Hess et al., 1980]. The combi- 
nation of these phenomena appear to result in a hemispherical 
asymmetry in the annual-total amount of solar and infrared 
radiation which strikes the residual polar caps, favoring the 
survival of CC >2 ice at the south pole [Lindner, 1990]. 

Modeling Procedure 

Figure 2 shows schematically all of the processes included 
in the model. The model is similar to those previously devel- 
oped except as specified below. The sources and sinks of 
energy per unit area of surface include solar radiation, modified 
for the absorption and scattering due to cloud and dust; infrared 
emission by cloud, dust, and gas in the atmosphere; infrared 
emission by the surface to space; atmospheric heat transport; 
heat conduction into/from the subsurface; and latent heat of 
condensation of CO 2 . The energy balance at the surface is 
represented by 

F{x,t)A{x,i) - ecT'Ho.x./) + L bM(x.t) - KhT{z,x,t) + H(.x,l) = 0 (1) 

8 r 5 1 

where F(x,t) is the solar and infrared radiation which strikes the 
surface at latitude x and time 1; A(x,t) is the co-albedo (1- 
albedo); E is the effective emissivity of the surface; o is the 
Stefan-Boltzman radiation constant; T(zjc,l) is the temperature 
and 8775z is the gradient in temperature with depth z, evaluated 
at the surface; L is the latent heat of condensation of CO 2 : 
5M(x,t) is the change in mass of CO 2 frost; 8 / is the time 
interval (taken to be 1/50 of a sol); K is the thermal conductiv- 
ity of the surface and the subsurface; and // (x,t) is the transport 
of heat by the atmosphere. The net gain or loss of energy over 
the time interval is used to compute either a change in surface 
temperature or a change in the mass of CO 2 frost. The solar 
and infrared radiation which strikes the surface as a function of 
latitude, lime, and cloud and dust amount is taken from Lindner 
[1990). The radiative effects of O 3 were suggested to be impor- 
tant for the energy budget [Kuhn et al., 1979], but have been 
shown to be negligible [Lindner, 1991]. Penetration of radia- 
tion into and through the ice has been shown to have a minor 
effect on the energy budget [Lindner, 1992]. The solar zenith 
angle dependence of ice albedo is taken from Warren el al. 
[1990], using a solar zenith angle which is reduced by 10° as a 
rough correction for the effect of surface roughness. The 
albedo of ice in the infrared is taken to be zero. Heat transport 
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Fig. 1. The seasonal recession of the south polar cap as observed over 
the last 20 years [James and Lummc. 1982] and as predicted by Leighton 
and Murray [1966], Briggs [1974], Davies el al. [1977], Narumi [1980], 
and James and North [1982]. (The aerocentric longitude of the Sun, L s . 
is the seasonal index; L s = 0°, 90°, 1 80°, and 270° correspond to northern 
spring equinox, summer solstice, autumnal equinox, and winter solstice, 
respectively.) 



Fig. 2. Schematic of the model for the energy budget of the polar cap, 
showing the physical processes included. The sum of all sources and 
sinks of energy is balanced by either a change in temperature or a 
change in ice mass. 


by the atmosphere is parameterized as by James and North 
[1982], that is, it simply depends on the latitudinal gradient of 
surface temperature multiplied by a constant [see also North el 
al., 1981; Hoffert el al., 1981). Hence the heat transport term 
is most important near the edge of the polar cap in spring and 
summer when the surface temperature gradient is very large. 
This is obviously an oversimplification of the transport pro- 
cess, but is easily tractable and has produced satisfactory 
results in the past. Conduction of heal into the subsurface is 
represented by the heat conduction equation 

ZT(zs,t)/St = ( K/pC)(b 2 T(z,x,t)/&z 2 ) (2) 


where p is the bulk density of the subsurface material and C is 
its specific heat. Equations (1) and (2) are applied and evalu- 
ated as by Jakosky and Haberle ( 1 990], every 2° in latitude. 

Clouds are assumed to exist only where CO 2 frost exists. 
The existence and opacity of cloud are highly variable and 
uncertain, although observers have noted a “polar hood", a 
layer of clouds which cloaks the polar cap [c.g., Leovy et al., 
1972; James et al., 1987). Several values for cloud opacity are 
studied here. Atmospheric dust opacity is taken from Pollack el 
al. [1979], although the global dust storms are ignored because 
of the strong increase in atmospheric heat transport. However, 
global dust storms do not appear to have significantly affected 
polar cap recession [ Briggs , 1974; Davies, 1979; James et al., 
1979; Kieffer, 1979; Marlin and Kieffer, 1979; James and 
North, 1982; Paige, 1985). Furthermore, the annual cycle in 
atmospheric pressure has repeated remarkably well from year to 
year, even though dust storms do not, which also implies that 
any linkage between dust storms and the CO 2 cycle is not 
strong. 


Comparison of the Model Results 
to Observations 

Figures 3-6 show the regression of the polar caps as pre- 
dicted by the model and as observed. Agreement is good, par- 
ticularly in the northern hemisphere. An exact comparison is 
difficult, considering that the edge of the polar cap is usually 
patchy and ill-defined [Leovy el al., 1972; Christensen and 
Zt trek, 1983], in large part due to terrain. The edge of the polar 
cap is also diurnally variable since ice frequently forms at 
night and sublimes during the day. There is also some year-to- 
year variability in polar cap regression [Iwasaki et al., 1990; 
James et al., 1990], Observations hint that cloud opacity is of 
the order of 0.2 in the southern hemisphere and 0.5 in the 
northern hemisphere [e.g., Briggs and Leovy, 1974; James et 
al., 1987], perhaps due in part to airmass differences. Such an 
asymmetry in cloud is used in Figures 3-6. An asymmetry in 
cloud opacity of 0.5 and 1.0 in the southern and northern 
hemisphere, respectively, will result in similar recession 
curves, although a higher ice albedo is required to counteract 
the increased polar night heating which results from increased 
cloud opacity. If no asymmetry in cloud opacity is assumed, 
then the model predicts seasonal CO 2 ice is more likely to be 
retained at the north pole than the south pole, as noted in 
earlier models. 

These results are obtained with an albedo for old ice of 
0.72. This value is consistent with those obtained by earlier 
models [e.g., James and North, 1982; Jakosky and Haberle, 
1990] and, given model simplifications and uncertainties in 
albedo measurements, would appear to be consistent with 
observations [Jakosky and Haberle, 1990], Fresh ice is 
assumed to have an albedo of 0.95. Based on terrestrial 
analogs, ice albedo is assumed to shift linearly from new to old 
over a period of 20 sols after net sublimation has started. This 
assumption did not result in significantly different results than 
just assuming an instantaneous shift from new to old ice. The 
inferred albedos are not unique, as there are many combinations 
and permutations of free parameters that will yield similarly 
good fits (e.g., changes in emissiviiy, cloud cover, heat trans- 
port, etc.). In fact, this lack of uniqueness did not allow for 
any determination of a latitudinal or seasonal dependence of ice 
albedo, and therefore given the lack of data on the albedo, no 
latitudinal or seasonal dependence in ice albedo was assumed. 
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Fig. 3. The regression of the south polar cap, as observed for various years (taken from Iwasaki et al. [1990]) and as simulated 
by my model (thin line), as a function of the aerocentric longitude of the Sun ( L s ). The cap radius is that which would be mea- 
sured on a polar stereographic projection of the south polar region; the units of the radius are fractions of the planetary radius 
of Mars. 
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However, a hemispherical asymmetry in some parameter, such 
as in the albedo or in the opacity of cloud and dust, is needed to 
allow the models to retain some seasonal CO 2 ice year-round at 
the south pole and not at the north pole. 

Atmospheric temperatures used in this model ranged from 
130 K to 180 K [ Lindner , 1990], some of the coldest inferred 
in the polar regions by the Viking orbiter [Lindal el al., 1979]. 
Higher atmospheric temperatures than those used here would 
actually increase the thermal emission by dust and cloud even 



Fig. 4. The regression of the south polar cap, as observed in 1986 (solid 
circles), 1971 (crosses), and 1977 (pluses) (taken from James el al. 
[1990]) and as simulated by my model (thin line), as a function of the 
acrocentric longitude of the Sun ( L s ). Two-sigma error bars are indi- 
cated for the 1986 data; the errors are smaller for the denser 1971 data 
and for the 1977 Viking data. 


more, making their radiative effects yet more important to the 
overall energy budget However, the temperatures used are not 
as important as one might expect because most atmospheric 
dust and cloud lie near the surface (about 60% is within 8 km of 
the surface; one atmospheric scale height at 160 K), and atmo- 
spheric temperatures near the surface Temain very close to the 
CO 2 frost point temperature during winter, except during 
global dust storms [Lindal et al, 1979]. 

Discussion 

If a hemispherical asymmetry in polar clouds is included in 
an energy balance model together with the solar zenith angle 
dependence of ice albedo, then the observations of polar cap 
regression can be reproduced while retaining seasonal CO 2 ice 
at the south pole and not at the north pole. The addition to the 
model of the radiative effects of cloud and background dust and 
of the solar zenith angle dependence of ice albedo causes little 
change in the total flux absorbed by the polar cap near its edge 
but becomes more important farther poleward from the polar 
cap edge [Lindner, 1990]. This would explain why earlier 
models were successful at reproducing the regression of the po- 
lar caps equatorward of 80° latitude but not at the pole (see 
Figure 1). Moreover, the integral of seasonal CC>2 ice over the 
80° to 90° latitude band is an order of magnitude less than the 
integral of seasonal CO 2 ice over the planet. Thus any effect 
which is most important only for ice at the pole will not have a 
major effect on the atmospheric pressure, explaining the good 
fit to the atmospheric pressure data also obtained by earlier 
models. 

James et al. [1979] and Paige and Ingersoll [1985J present 
evidence for the existence of a hemispherical asymmetry in 
polar cap albedo, with the south pole being 25% brighter. 
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SOLAR LONGITUDE 


Fig. 5. The regression of the north polar cap. as observed for various years (taken from Iwasaki et al. (1982); symbols refer to 
Dollfus (1973], Baum (1974], Briggs ( 1974], James [1979], Capen and Parker [1981 ], and Iwasaki el al. [1982]) and as simu- 
lated by my model (thin line), as a function of the aerocenlric longitude of the Sun (L s ). 


although the specifics of how a hemispherical asymmetry in 
ice albedo maintains itself remain to be completely worked 
out. If a hemispherical asymmetry in albedo plays an impor- 
tant role in the asymmetry in the polar caps, then the asymme- 
try in albedo would have to be as extreme as that suggested by 
Paige and Ingersoll [1985]. A substantial part of the energy 
budget of the polar caps comes from IR emission by cloud and 
dust, atmospheric heat transport, and heat conduction from the 
subsurface; sources which arc unaffected by albedo variations at 
solar wavelengths. Furthermore, polar cap albedo does not 
have a linear effect on the energy budget, as cloud and dust will 
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Fig. 6. The regression of the north polar cap, as simulated by my model 
(thin line) and as observed, as a function of the acrocentric longitude of 
the Sun ( L s ). Data are taken from James el al. [1987]; also showing data 
of Iwasaki el al. (1979; 1982]. The symbols stand for 1977-1978 and 
1980 Viking data (crosses and pluses, respectively) and ground-based 
1975-1980 data (circles). 


reflect some of the solar flux reflected from the polar cap back 
onto the polar cap. For these two reasons, the presence of 
cloud and dust lessens the ability for albedo changes to affect 
CC>2 ice lifetimes. 

Many processes have yet to be examined as to their effect 
on the energy budget of the polar caps. The emissivity of CO 2 
ice could be much lower, based on the theoretical simulations 
of Warren et al. [1990]. Ice is strongly forward scattering at 
high solar zenith angles [Taylor and Stowe, 1984], which has 
implications for multiple reflections between polar cap and 
polar hood. Surface roughness, particularly if present in the 
form of penitentes [Svilek and Murray, 1988], would change 
the effective solar zenith angle and IR surface cooling. Ice 
microphysics also requires further study with regard to the ice 
energy budget and may offer some hemispherical asymmetry 
[J. Eluszkiewicz, personal communication, 1992]. Snowfall 
also needs to be accounted for, particularly during periods of 
high dust opacity [Pollack el al., 1990], An insulating residue 
on the ice [Saunders et al., 1986] and wind shifting of ice 
[Briggs, 1974; James et al., 1979; Kieffer, 1979] may also be 
significant for the ice energy budget. Improved modeling of 
atmospheric heat transport is also needed. Possible hemi- 
spheric asymmetries in weather systems may contribute to the 
asymmetry in the polar caps [Barnes, 1991). Moreover, more 
observations of the polar hood are needed to show the magni- 
tude of a hemispherical asymmetry in cloud opacity, particu- 
larly in polar night, and to answer questions with regard to 
interannual variability. 

Despite the work that remains, it appears possible to con- 
struct an energy balance model which maintains seasonal CO 2 
ice at the south pole year-round and not at the north pole year- 
round and still reasonably simulates the overall polar cap 
regression data. Hence earlier model difficulties in retaining 
seasonal ice at the south pole have been removed. CO 2 ice 
observed in the summertime south polar cap could be seasonal 
in origin. Seasonal CO 2 ice at the south pole is also suggested 
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by the water vapor observations of Jakosky and Barker (1984]. 
However, further research remains before it is certain whether 
the CO 2 ice observed in the summertime south pole is seasonal 
or is part of a permanent reservoir. 
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[ Jupiter's Celestial Necklace 



the motion of a very rapidly mov- 
ing object; for repetitive patterns 
such as the spokes of a wheel, the 
speed can be measured by adjust- 
ing the light’s flash rate to mimic 
the time it takes for one spoke to 
advance exactly to the position of 
another. This effect would not 
apply here because sunlight has 
no such regular (or even irregular) 
flickering. Neither are the Voyager 
cameras susceptible to this sort of 
effect. 

There is surely a lot that we 
don’t understand about Saturn’s 
spokes and about the rings in gen- 
eral. When Cassini arrives at Sat- 
urn in 2004, it will observe the flux 
of interplanetary material with its 
dust detector, observe the changes 
in the magnetic fields near the 
rings with charged-particle detec- 
tors, and observe the structure, for- 
mation, color and evolution of the 
spokes. It will discover whether 
the distribution of spokes varies 
with time and tilt of the rings, as 
might be expected based upon a 
meteoroid impact hypothesis. 

Cassini will also answer many 
other questions about the rings: 
how the composition varies from 
dark, grayish material in the C ring 
and Cassini Division to bright, red- 
dish material in the B and A rings; 
whether small moonlets reside in 
the many empty gaps in the rings 
(other than the one 10-kilometer 
object discovered recently in the 
A ring's Encke gap); whether new 
ring features have appeared or 
whether some of the features 
Voyager saw have since evolved 
into different forms. 

— JEFFCUZZI. MSA Ames 
Research Center 

Does Mars have an ozone layer 
similar to that of Earth? If so, 
are there fluctuations in the layer 
(presumably due to natural causes ) 
similar to those resulting in the 
ozone hole(s ) currently causing 
so much concern on Earth? 

— Mike Clean’, Jefferson City, 
Missouri 

Mars' ozone has been studied 
from Earth-based observatories 
and from Russian and American 
spacecraft, most notably Mariner 9. 
Mars does in fact have an “ozone 
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layer” at equatorial and midlati- 
tudes, much as Earth does. 

In those regions. Mars' ozone 
concentrations are highest at an 
altitude of around 30 to 40 kilo- 
meters (19 to 25 miles), with a 
maximum concentration of about 
0.5 parts per million (ppm). For 
comparison, ozone concentrations 
on Earth reach a maximum of 
about 5 ppm. However, the verti- 
cally mixed amount of ozone is 
many, many times less than it is 
on Earth. 

Hence, while ozone on Earth 
shields the surface from much of 
the Sun’s harmful ultraviolet radia- 
tion, the ozone on Mars does not. 

Furthermore, while the ozone in 
Earth’s stratosphere exerts a strong 
control on the temperatures there. 

Mars' ozone exerts only a negligi- 
ble control on the planet’s atmo- 
spheric temperatures. 

At Mars’ poles in winter, ozone 
does not form a “layer," since it 
also exists in high quantities near 
the surface. The vertically inte- 
grated amount of ozone is higher 
at the poles than at equatorial lati- 
tudes in winter. Thus, ozone on 
Mars fluctuates much more than 
what we’ve observed on Earth, 
but for different reasons. 

Earth's ozone hole is related in 
pan to human-made chemicals 
that are not present on Mars. There, 
concentrations are more strongly 
influenced by the amount of water 
vapor, and the wide variability of 
water vapor on the planet results in 
a large fluctuation in the amount 
of ozone. 

This explains why ozone on 
Mars is more abundant at the very 
cold wintertime poles, where water 
vapor is frozen out of the atmo- 
sphere. In fact, observations by 
Mariner 9 in 1972 showed that 
ozone rose and fell in conjunction 
with the passage of meteorological 
cold fronts, due in part to the re- 
duction in water vapor in colder air 
and in part to the obscuration of 
ozone by clouds and dust. 

Our group is proposing to 
NASA that a balloon-borne ozone 
detector be sent to Mars to study 
this variability more closely, as a 
possible Discovery-class mission. 

— B. LEE LINDNER. Atmospheric 
and Environmental Research. Inc. 
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Introduction 

C0 2 is the only gas which produces appreciable infrared 
cooling of the Martian troposphere (Goody and Belton, 
1967 ; Crisp, 1990 ; Lindner, 1991 ; Savijarvi, 1991 ; Hour- 
din, 1992); hence, the transmission on Mars of the C0 2 
15 pm band has been well studied (Gal’tsev and Osipov, 
1979; Crisp, 1990; Lindner et at., 1990a,b; Savijarvi, 
1991 ; Hourdin, 1992). Airborne dust also produces 
appreciable cooling, even in periods other than global dust 
storms (e.g. Lindner, 1985, 1993). However, virtually all 
previous studies of Martian radiative cooling have ignored 
the presence of dust when calculating C0 2 cooling rates. 
The effect of dust-C0 2 spectral overlap (hereafter abbrevi- 
ated as dust-C0 2 interaction) was studied by Kondrat’ev 
and Moscalenko (1975) and Egan et al. (1980) for the 
effect on line profiles, by Pallman (1983) for the near- 
surface boundary layer, by Pollack et al. (1979) for invert- 
ing data, and by Habcrle et al. (1982) for circulation 
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models. However, none of these studies published or dis- 
cussed dust-C0 2 interaction. The importance of both C0 2 
and dust in Martian radiative transfer requires a model 
capable of accurately handling the scattering by dust and 
the complex wavelength structure of C0 2 . To simplify 
the problem, most prior work has studied C0 2 and dust 
radiative transfer separately and avoided their over- 
lapping opacities. Here I examine the errors that can result 
from such a treatment. 


Model 

The discrete ordinate method of Stamnes et al. (1988) is 
used to treat the scattering, absorption and emission of 
monochromatic radiation through the Martian atmo- 
sphere. The exponential sum program of Evans et al. 
(1980) converts the banded wavelength structure of C0 2 to 
allow for monochromatic treatment (Freeman and Liou, 
1979 ; Lindner, 1985 ; Lindner et al., 1990a,b). The 15 pm 
band is broken into 16 subintervals to properly account 
for the multiple scattering (Lindner et al., 1990a,b ; Laris 
and Oinas, 1991). Transmission functions for the 15.0 pm 
band of C0 2 are taken from the Iine-by-Iinc model results 
of Gal’tsev and Osipov (1 979). The transmission function, 
Tr, as a function of temperature T, pressure P, and C0 2 
column abundance U is extrapolated from the Gal’tsev 
and Osipov results (subscript G) to temperatures below 
200K by: 

TriT, P,U) = 1 — {(1 — 7>7200K, P, U)) (77200K) 0 ‘}. 

(I) 

The exponential 0.8 is found when the temperature depen- 
dencies for the Pollack et al. (1981) transmission functions 
are recast in this form. Using a modified version of the 
FASCOD transmission model (Clough et at., 1986), the 
accuracy of the transmission functions in equation (1) 
arc confirmed for the range in temperature and pressure 
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present in the atmosphere to 40 km altitude (Lindner et 
al., 1990a,b; Pollack et al., 1990). Local thermodynamic 
equilibrium is assumed (Gierasch and Goody, 1967 ; 
Hourdin, 1992). 

Dust opacities vary from 0.2 to 1.0 for conditions other 
than global dust storms (Pollack et al., 1979). However, 
dust opacities over winter polar latitudes may be slightly 
less (e.g. Lindner, 1990). A Gaussian profile describes the 
vertical distribution of dust, being well mixed to 20 km 
altitude for conditions other than global dust storms 
(Anderson and Leovy, 1978; Zurek, 1982; Korablev et 
al., 1993). The wavelength dependence of the dust opacity 
and the dust single-scattering albedo are given by Toon et 
al. (1977). Scattering of radiation by dust is represented 
by the Haze-L phase function (Toon et al., 1977). Com- 
putational difficulties which accompany highly asym- 
metric phase functions are removed with the Delta-L 
method (Wiscombe, 1977). The emissivity of airborne dust 
is high and has been calculated as a function of wavelength 
from theory and observations (Toon et al., 1977 ; Simpson 
et al., 1981). C 

— ^ Atmospheric composition is taken as 95% gQ 2 (Owen 

et al., 1977). Atmospheric composition may have been 
quite different in past epochs, with C0 2 being perhaps a 
minor constituent (e.g. Lindner and Jakosky, 1985), but 
this study focuses on the present epoch. Season-dependent 
Ct\ abundances (Hess et al., 1980) are corrected for 
cir«. ilation-induced pressure gradient (Haberle et al., 
1979) and elevation (Lindal et al., 1979). The surface 
pressure is 8 mbar at 57°N latitude in late winter. Atmo- 
spheric properties are zonally averaged and are assumed 
azimuthally independent. The region from the surface to 
40 km altitude is broken into 20 layers, each 2 km thick, 
to account for vertical inhomogeneity. The temperature 
profile rises linearly with altitude from 150K at the surface 
to 160K at 10 km, and then falls linearly with altitude to 
130K and 40 km, typical for winter polar conditions (Lin- 
dal et al., 1979 ; Kieffer, 1979 ; Martin, 1984). The albedo 
of the polar cap is assumed to be zero (Kieffer, 1970; 
Smythe, 1975; Wiscombe and Warren, 1980). An average 
ice emissivity of 0.9 is adopted, with the wavelength depen- 
dence given in Ditteon and Kieffer (1979) and Hunt et al. 
(1980). 


Results and discussion 

Total atmospheric cooling rates owing to thermal emis- 
sion by a vertical optical depth of dust, x r , of 0.2, and the 
contribution by various wavelength intervals are pre- 
sented in Fig. 1 for winter polar conditions. The maximum 
cooling owing to dust below 30 km altitude occurs in the 
12- 1 8 pm wavelength interval. The strong dust cooling in 
the 12-18 pm range is owing to the cold temperature, to 
the high optical depth of dust at 12-18 /mi (Toon et al., 
1977), and to the secondary maximum in dust emissivity 
at 12-18 pm (Toon et al., 1977). Hence, COj impedes the 
ability of dust to cool the atmosphere because C0 2 absorbs 
strongly in the 12-18 pm region where dust cools most 
effectively. 

The presence of dust will affect C0 2 radiative transfer 



Dust Cooling Rate (K/Day) 


i Fig. 1. Total infrared cooling rate by dust at winter polar lati- 
: tudes (57°N latitude, L, = 343°) as a function of altitude for 
j dust vertical optical depth, of 0.2, and the contribution from 


each wavelength interval from 5 to 40 pm (wavelength range of 
each bin is labeled in pm) 


by increasing the optical depth, and absorbing light that 
is emitted by C0 2 . Dust is particularly important in the 
radiative cooling in the line wings of C0 2 . Line wings are 
important to cooling in that the strongly absorbing line 
center is often optically thick, which renders the line center 
ineffective at cooling (e.g. Lindner, 1993). Optically thin 
line wings cool effectively, in that photons emitted in the 
line wings can escape and cool the atmosphere. The pres- 
ence of dust affects the line wings by absorbing emitted 
photons that would otherwise escape. 

Hence, dust and C0 2 hinder each other’s attempts to 
cool. Thus, simply adding dust cooling to C0 2 cooling is 
not accurate in those wavelength bands where both are 
important. This competition affects C0 2 cooling more 
than dust cooling. Dust cooling occurs over the entire 
infrared spectrum, which is mostly free of C0 2 absorption. 
However, C0 2 absorption and emission must always com- 
pete with the absorption by dust. Figure 2 shows that the 
cooling is overestimated up to 30% if dust and C0 2 
cooling rates are computed separately and then summed. 
The effect is strongest near the surface where the optical 
depth is largest. C0 2 cooling is optically thick, making it 
very difficult for emitted photons to escape in the lowest 
10 km. Some of the photons emitted by C0 2 line wings 
that would normally escape and cool are being absorbed 
by dust, which results in a lower cooling rate. C0 2 also 
absorbs some of the photons emitted by dust that would 
normally escape and cool. 

The error induced by ignoring dust-C0 2 interaction is 
increased at higher dust loading (Fig. 3). Dust is now even 
more capable of choking off the escape of emitted photons, 
not only those of C0 2 , but also photons emitted by dust 
itself. Dust-C0 2 interaction is highest near the surface for 
all three cases. The relative importance of the line wings 
for cooling is greater at lower altitudes. Since dust has its 
strongest effect on the line wings, dust-C0 2 interaction is 
strongest near the surface. Note that the errors do not 
change as much when dust opacity is changed from 
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15 Micron Region Cooling Rate (K/Day) 

Fig. 2. Total 15 /on region (12-18 pm) cooling rates as a function 
of altitude for a pure COj atmosphere (labeled “gas only”) and 
for a pure dust (t, = 0.2) atmosphere (labeled “dust only") at 
winter polar latitudes. Also shown are the cooling rates for a 
simple sum of the C0 2 only and dust only cooling rates (labeled 
‘.‘gas only + dust only"), and the cooling rates calculated when 
dust and C0 2 are treated simultaneously (labeled “simultaneous 
gas & dust”) 

r. = 0.5 to 1 .0 as they do when dust opacity is changed 
from t v = 0.2 to 0.5 or from 0.0 to 0.2. This occurs because 
line wings are effective in optically, thin media. Once large 
dust opacity is reached, further increases in dust opacity 
do not affect line wing optical depths as much. 

Note that the error increases with dust optical depth, 
even when dust cooling exceeds C0 2 cooling. The expla- 
nation lies in the two-fold nature of dust-C0 2 interaction, 
and underlies the importance of its inclusion even for 
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Ratio of Dust + Gas / Dust & Gas 


Fig. 3. Ratio of the 15pm region (12-18 pm) thermal cooling 
rates computed ignoring dust-C0 2 interaction (labeled “gas 
only + dust only” in Fig. 2) to those computed including dust- 
COj interaction (labeled “simultaneous gas & dust" in Fig. 2) i 
for r,(dust) of 0.2, 0.5 and 1.0 at winter polar latitudes. Hence, ' 
these curves represent the overestimation error of computing 
cooling rates in the 12-18 pm wavelength interval without pro- 
per simultaneous treatment of dust and C0 2 , as a function of I 
altitude and dust opacity 


calculations for dust storm conditions. For low dust opac- 
ities, dust-COj interaction involves the effect of dust in 
choking off emission in the line wings of the C0 2 15 pm 
band. However, for large dust opacities, it is C0 2 that 
is the parasite. For large dust opacities, dust cooling is 
dominant and as much as half of all dust cooling occurs 
in the 1 5 pm wavelength interval (Fig. 1). However, dust is 
unabte to cool effectively in the 1 5 pm wavelength interval 
because emitted photons are quickly reabsorbed by C0 2 . 
In fact, the parasitic effect of C0 2 grows as dust opacity 
increases, as shown in Fig. 3. The error reaches an asymp- 
tote at large dust opacities. 


Summary and conclusions 

As much as half of all atmospheric cooling owing to ther- 
mal emission by dust occurs in the 12-18 pm wavelength 
range. Since this is also the wavelength region where most 
cooling owing to thermal emission by C0 2 occurs, it is 
important to consider the spectral overlap between dust 
and C0 2 . C0 2 infrared cooling and dust infrared cooling 
are both overestimated if C0 2 and dust cooling are com- 
puted separately and then summed. Indeed, ifdust-C0 2 
interaction is ignored, then the total 12-18 pm cooling is 
overestimated by almost 30% for minimal dust loading 
at winter polar latitudes. Furthermore, for one vertical 
optical depth of dust at winter polar latitudes, a simple 
sum of dust cooling and C0 2 cooling from 12-18 pm will 
overestimate the actual total 12-18 pm cooling by over 
40% near the surface. Hence, dust-C0 2 interaction is 
important not only for calculating C0 2 cooling rates, since 
dust opacities overlap the C0 2 1 5 pm band, but also for 
calculating dust cooling rates, since the 15 pm band of 
C0 2 overlaps a significant portion of the infrared region 
where dust cooling is most effective. In other words, even 
when the Martian atmosphere is very dusty, when CO : 
cooling and heating rates are minor compared with dust 
(e.g. Lindner, 1993), the interaction of C0 2 will sig- 
nificantly affect the ability of dust to cool the atmosphere. 
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ABSTRACT 


A computer model is used to study the radiative transfer of the martian winter-polar atmo- 
sphere. Solar heating at winter-polar latitudes is provided predominately by dust. For normal, 
low-dust conditions, CO2 provides almost as much heating as dust. Most heating by CO2 in the 
winter polar atmosphere is provided by the 2.7 Jim band between 10 km and 30 km altitude, and 
by the 2.0 |im band below 10 km. The weak 1.3 }im band provides some significant heating near 
the surface. The minor CO2 bands at 1.4, 1.6, 4.8 and 5.2 |im are all optically thin, and produce 
negligible heating. O3 provides less than 10% of the total heating. Atmospheric cooling is 
predominantly thermal emission by dust, although CO2 15 (im band emission is important above 
20 km altitude. 

Key Words: Mars, Radiative Transfer, Atmosphere, Planets 
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INTRODUCTION 


While the past 20 years of spacecraft exploration have expanded our knowledge of the 
planet Mars, it seems more interesting problems exist today than ever before. At the heart of 
most of these mysteries is the winter polar atmosphere. Observations of the winter polar atmo- 
sphere have been limited by orbital constraints and darkness, and few models have been con- 
structed to date. CO 2 and H 2 O condensation commonly occur in the winter polar atmosphere of 
Mars, resulting in extensive cloudiness and a massive ice sheet. The structure and composition 
of the polar clouds as functions of latitude, altitude and season are poorly understood. Although 
observations of the southern polar hood are limited, it appears strikingly different from its north- 
ern counterpart, with much less coverage and different altitude structure (James, 1983; 
Christensen and Zurek, 1983; Martin and James, 1984; Akabane et al., 1990). The seasonal 
recession of the south polar cap (e.g., Iwasaki et al., 1990) cannot be accurately reproduced with 
energy balance models without consideration of atmospheric radiative effects (Narumi, 1980; 
James and North, 1982; Lindner, 1990; 1991a; 1992a; 1992b; 1993a), and the relative fraction of 
snow and frost in the cap is unknown (Pollack et al., 1990). Atmospheric dynamics at winter 
polar latitudes is also dependent on atmospheric heating and cooling (Haberle et al., 1979; 1982). 
Ozone abundances on Mars are also dependent on the radiative effects of dust (Lindner, 1988), 
and the observance of ozone on Mars is impaired by the radiative effects of dust (Lindner, 
1992c). 

Understanding the dominant radiative heating and cooling mechanisms in the winter polar 
atmosphere is crucial to solving many of these mysteries. While dynamical and latent heat 
mechanisms do provide as much as half of the atmospheric energy budget (Pollack et al., 1990), 
it is the radiative transfer mechanisms which drive the dynamical and latent heat mechanisms. In 
fact, due to the strong CO 2 heating and cooling, radiative processes are relatively more important 
in determining the temperature structure of the martian atmosphere than of the terrestrial 
atmosphere (e.g., Pollack et al., 1990). This work intends to establish the relative importance of 
O 3 , CO 2 and dust in the radiative heating and cooling of the winter polar atmosphere of Mars, 
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studying the importance of all wavelength bands. CO 2 was shown to be the only gas which 
produced appreciable infrared cooling at winter polar latitudes by Goody and Belton (1967), 

Crisp (1990), and Savijarvi(1991). Also, significant solar heating occurs in all the near-infrared 
(NIR) bands of CO 2 (Pollack et al., 1981; Lindner, 1985; Savijarvi,1991). Ozone was suggested 
to be an important contributor to the atmospheric temperature in the polar regions by Kuhn et al. 
(1979). However, the contribution of ozone to the atmospheric heat budget was later shown to 
be minor when compared to the contribution of dust (Lindner, 1991b). However, atmospheric 
heating due to dust absorption of solar radiation was shown to be important during dust storms 
by Gierasch and Goody (1972), Moriyama (1975), Zurek (1978), and Davies (1979). The 
importance of both CO 2 and dust in infrared radiative transfer, as well as the interaction between 
gas and dust, was demonstrated for certain cases by Kondrat'ev and Moscalenko (1975). The 
importance of both gas and dust in martian radiative transfer requires an advanced model capable 
of accurately handling the scattering by dust and the complex wavelength structure of carbon 
dioxide. To simplify the problem, most prior work has studied gas and dust radiative transfer 
separately and avoided their overlapping opacities, a technique which is inaccurate (Lindner, 
1993b). To more properly assess the relative importance of dust and CO 2 in the thermal budget, 
this study treats them simultaneously. 

MODEL 

The discrete ordinate method of Stamnes et al. (1988) is used to treat the scattering, 
absorption and emission of monochromatic radiation through the martian atmosphere. The 
exponential sum program of Wiscombe and Evans (1977) and Evans et al. (1980) converts the 
banded wavelength structure of CO 2 to allow for monochromatic treatment (Lindner, 1985; 
1993b; Lindner et al., 1990a, b). Lindner (1993b) clearly shows that errors as large as 50% occur 
when dust and CO 2 cooling and heating rates are computed separately and then summed. Hence, 
an approach such as the exponential sum technique which allows for simultaneous treatment of 
CO 2 and dust in the solution of the radiative transfer equation is necessary in order to properly 
assess their relative importance (Lindner, 1993b). Transmission functions for the 2.759, 4.301 
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and 14.93 (lm bands of CO 2 (hereafter abbreviated as 2.7, 4.3, and 15.0, respectively) are taken 
from the line-by-line model results of Gal’tsev and Osipov (1979). The transmission function, 

Tr, as a function of temperature T, pressure P, and CO 2 column abundance U is extrapolated 

from the Gal'tsev and Osipov results (subscript G) to temperatures below 200 K by 

Tr(T,P,U) = 1 - [(1 - Tr G (200 K, P, U)) (T/200 K)Q] (1) 

The exponential Q is found to be 0.45, 0.3, and 0.8 for the 2.7, 4.3, and 15 |±m bands, respec- 
tively, when the temperature dependencies for the Pollack et al. (1981) transmission functions 
are recast in this form. Using a modified version of the FAS COD transmission model (Clough et 
al., 1986), the accuracy of these transmission functions is confirmed, and transmission functions 
are obtained for the 1.316, 1.455, 1.600, 2.020, 4.840 and 5.200 Jim bands of CO 2 (hereafter 
abbreviated as 1.3, 1.4, 1.6, 2.0, 4.8, and 5.2, respectively) covering the range in temperature and 
pressure present in the atmosphere to 40 km altitude (Lindner et al., 1990a,b). Additionally, 
ozone absorption cross-sections from 1500 A to 3200 A (Daumont et al., 1983; Freeman et al., 
1984) and from 4000 A to 8000 A (Griggs, 1968) and ultraviolet (UV) cross-sections for CO 2 
(Shemartsky, 1972; Lewis and Carver, 1983), O 2 (Demore et al., 1988), H 2 O (Thompson et al., 
1963; Hudson, 1971), HO 2 (Demore et al., 1988) and H 2 O 2 (Demore et al., 1988) are included 
[e.g., Lindner, 1988; 1991b]. 

Dust opacities vary from 0.2 to 1.0 for conditions other than global dust storms (Pollack 
et al., 1979; Lumme and James, 1984). However, dust opacities over winter polar latitudes may 
be slightly less [e.g., Lindner, 1990]. A gaussian profile describes the vertical distribution of 
dust, being well-mixed to 20 km altitude for conditions other than global dust storms (Anderson 
and Leovy, 1978; Zurek, 1982; Korablev et al., 1993). Dust storm conditions are not considered 
here because of the dramatic increase in dynamical processes during dust storms. The 
wavelength dependence of the dust opacity is given by Toon et al. (1977). The single scattering 
albedo of airborne dust as a function of wavelength is given by Zurek (1978; 1982) and Toon et 
al. (1977) for solar and infrared wavelengths, respectively, using a solar average of 0.9 (Clancy 
and Lee, 1991). Scattering of radiation by dust is represented by the Henyey-Greenstein phase 
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function (Toon et al., 1977; Clancy and Lee, 1991). Computational difficulties which 
accompany highly asymmetric phase functions are removed with the Delta-M method 
(Wiscombe, 1977). The emissivity of airborne dust is high and has been calculated as a function 
of wavelength from theory and observations (Toon et al., 1977; Simpson et al., 1981). Dust 
optical properties in the near-IR (1-5 |im) are highly uncertain, hence making the calculated 
heating and cooling rates uncertain. However, as will be shown later, dust heating and cooling 
rates in the near-IR are minor, making the uncertainty in dust optical properties in the near-IR 
unimportant Dust optical properties in the 15 Jim region are also highly uncertain, and factor of 
2 uncertainty in the computed cooling rates is quite possible. Clouds will also affect atmospheric 
radiative transfer. However, since the cloud opacity is highly variable (i.e. Briggs and Leovy, 
1974), the cloud particle scattering properties are very uncertain, and even the composition of the 
clouds is unclear, the effect of clouds is highly speculative and variable. But clouds should affect 
atmospheric radiative transfer similarly to how dust does, since dust single scattering albedos are 
very high (Lindner, 1990; 1993a). 

The Rayleigh scattering optical depth is computed as in Hansen and Travis (1974), using 
parameters appropriate for Mars. Solar fluxes are taken from Smith and Gottlieb (1974) and 
Rottman (1981), after adjusting for the eccentricity and orbital radius of Mars. Solar heating 
rates are diumally averaged (e.g., Cogley and Borucki, 1976). Atmospheric properties are 
zonally averaged and assumed azimuthally-independenL The region from the surface to 40 km 
altitude is broken into 20 2-km-thick layers to account for vertical inhomogeneity. The improved 
Curtis-Godson approximation (Yamamoto et al., 1972; Ramanathan and Coakley, 1978) is used 
to treat vertical inhomogeneity at thermal wavelengths. The Chapman function is used to 
approximate the slant path in place of the secant function [e.g., Smith and Smith, 1972], because 
the winter polar atmosphere always has large solar zenith angles, and the secant function is in 
error for large angles. 

Atmospheric composition is taken as 95% CO 2 (Owen et al., 1977; see also Kondratiev et 
al., 1973). Atmospheric composition may have been quite different in past epochs, with CO 2 
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being perhaps a minor constituent (e.g., Lindner and Jakosky, 1985; Lindner, 1993c), but this 
study focuses on the present epoch. Season-dependent CO2 abundances (Hess et al., 1980) are 
corrected for circulation-induced pressure gradients (Haberle et al., 1979) and elevation (Jakosky 
and Farmer, 1982; Lindal et al., 1979). The surface pressure is 8 mbar at 57°N latitude in late 
winter, which is when the maximum O3 column abundance of 57 [im-atm was observed (Barth et 
al., 1973). The altitude dependence of O3 is based on model results (Lindner, 1988). 

As this study uses late northern winter conditions (Ls = 343°), the surface is covered by 
somewhat dirty ice with an albedo of 0.5 (Kieffer, 1979; James and Lumme, 1982). [Ls, the 
solar longitude, is a seasonal index; Ls of 0°, 90°, 180°, 270°, correspond to northern spring 
equinox, summer solstice, autumnal equinox, and winter solstice, respectively]. The wavelength 
dependence of the ice albedo is taken from Hapke et al. (1981) and Warren and Wiscombe 
(1980). The infrared (5.4 - 100 (lm) albedo of the polar cap is assumed to be zero (Kieffer, 1970; 
Smythe, 1975; Wiscombe and Warren, 1980). An average ice emissivity of 0.9 is adopted, with 
the wavelength dependence given in Ditteon and Kieffer (1979) and Hunt et al. (1980). The 
temperature profile rises linearly with altitude from 150K at the surface to 160K at 10 km, and 
then falls linearly with altitude to 130 K and 40 km, typical for winter polar conditions (Lindal et 
al., 1979; Kieffer, 1979; Martin, 1984). Atmospheric temperatures are poorly known above 30 
km altitude, and therefore results above that altitude are speculative and are not presented here. 
Local thermodynamic equilibrium is assumed (Gierasch and Goody, 1967; Uplinger et al., 1984; 
Hourdin, 1992). 

RESULTS AND DISCUSSION 

Figure 1 presents the ozone, carbon dioxide, and dust heating rates, and the carbon dioxide 
and dust cooling rates, for late winter (Ls = 343°) conditions at 57°N latitude with 0.2 vertical 
optical depths, Tv, of dust (averaged over the solar spectrum). The net heating and net cooling 
are virtually identical at all altitudes. Because thermal emission is a strong function of the tem- 
perature and heating is virtually independent of temperature, less than a 10K adjustment in the 
assumed temperature profile will yield a perfect balance of heating and cooling. These modifi- 
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cations would still be consistent with observations of temperature (Lindal et al., 1979). Because 
the variation in dust with altitude is not well understood, the discrepancy could also be due to the 
assumed altitude profile of dust. Indeed, the assumed temperature profile could be correct, and 
the dust profile could be extracted by obtaining a balance between heating and cooling. Note 
that in addition to radiative cooling, energy could also go into condensation. In fact, clouds are 
often observed at these latitudes and seasons. Clouds would also change heating and cooling 
rates by increasing the flux (and heating) up high, and decreasing flux (and heating) in the lower 
10 km. 

While the near equality of heating and cooling rates means that meridional and vertical 
heat transport may not be required to explain the observed winter polar temperatures, it certainly 
does not rule out any meridional or vertical heat transport. Indeed, vertical heat transport could 
also be responsible for the low cooling rates at 20-30 km altitude. Observations do show some 
day to day variability, which could be due to changes in dynamical heat transport, or to pockets 
of high ozone or dust concentrations which would change radiative heating and cooling. 
However, maximum dynamical heating rates are of the order lK/day, with typical winter polar 
heating rates much smaller (Gadian, 1978; Pollack et al., 1981). While this could be significant 
near the surface, it becomes less so at higher altitudes. 

The relative importance of gas and dust is clearly seen in Figure 1. Dust heating is the 
major source of heating at all altitudes, particularly above 20 km. Ozone provides approximately 
10% of the total heating at all altitudes. 57 Jim atmospheres of ozone are used, the maximum 
observed by Mariner 9 (Barth et al., 1973). [1 p.m atmosphere = 2.69 x 10 15 cm' 2 ]. Smaller 
ozone abundances will reduce the importance, but not dramatically. This contradicts earlier 
work by Kuhn et al. (1979) which showed ozone to be a more significant source of heating. 
However, Kuhn et al. (1979) ignored the heating by dust, which is clearly incorrect (Lindner, 
1991b). Indeed, 0.2 vertical optical depths of dust represent the minimum amount of dust 
observed (Leovy et al., 1972; Pollack et al., 1979; Thorpe, 1981; Zurek, 1981). Larger dust 
loading will provide greater heating. Ozone heating is almost the same as that of the NIR bands 
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of CO2, at all altitudes. This occurs because the major CO2 bands are saturated at the large solar 
zenith angles in the winter polar atmosphere. 

Figure 1 can also be used to show what would happen to the thermal structure for the case 
of no dust. CO2 and O3 heating is triple CO2 cooling near the surface, while CO2 cooling is 
triple CO2 and O3 heating at 30 km altitude. This means that a dust-free atmosphere would be 
warmer near the surface and cooler at 30 km to allow for a balance between heating and cooling. 
Higher lapse rates were found in other dust-free studies as well (e.g., Gierasch and Goody, 1968). 

Cooling rates are also dominated by dust. Dust cooling is greater than CO2 15 |im band 
cooling from the surface to 25 km altitude, with CO2 cooling dominant above this altitude. Note 
that cooling is not dependent on latitude, but on temperature. Hence the relative importance of 
dust cooling to CO2 cooling is approximately the same at all winter polar latitudes. 

The dominant ozone heating occurs at 2700 A, with appreciable contributions from 2200 
A to 3100 A (Lindner, 1991b). The Chappuis bands (4000 - 7000 A) provide over 10% of the 
total ozone heating near the surface. CO2, 02, H2O, HO2, and H2O2 produce minor heating of 
the winter polar atmosphere at UV wavelengths, although CO2 and O2 UV heating is appreciable 
above 30 km altitude (Lindner, 1991b). H2O, HO2, and H2O2 number densities are too low for 
any appreciable UV heating. 

CO2 solar heating rates at near-infrared (NIR) wavelengths at 57°N latitude are computed 
for each CO2 band, as shown in Figure 2. Most heating by CO2 in the winter polar atmosphere is 
provided by the 2.7 pm band between 10 km and 30 km altitude, and by the 2.0 pm band below 
10 km. The weak 1.3 pm band provides some significant heating near the surface. The minor 
CO2 bands at 1.3, 1.4, 1.6, 4.8 and 5.2 p.m are all optically thin, and produce negligible heating. 
The 2.0 pm band is more strongly absorbing than these minor bands and becomes optically thick 
at 10 km altitude, which results in a decreasing heating rate with a decrease in altitude below 10 
km. The 2.7 pm band is stronger yet, and becomes optically thick at 20 km altitude. The 4.3 pm 
and 15 |im bands are very efficient, and are optically thick even at 30 km altitude. 15 Jim band 
heating is surprisingly strong, despite the low solar flux at infrared wavelengths. The explana- 
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tion lies in the large bandwidth (from 12 to 19 Jim) and the efficient absorption over the band 
width. 

In addition to being heated from absorption of solar radiation by O 3 and CO 2 , the martian 
atmosphere is also heated by absorption of solar radiation by dust. As dust optical depths are not 
as strongly wavelength dependent as CO 2 and O 3 optical depths, the altitude dependence of dust 
heating is virtually the same for all wavelengths (Fig. 3). The strongest heating occurs in the vis- 
ible where the maximum in solar flux occurs. 

Atmospheric cooling rates due to thermal emission by 0.2 vertical optical depths of dust 
are presented in Figure 4. The maximum cooling occurs at 10 km altitude. Cooling is not as 
efficient in the lowest 5 km for two reasons. The temperature profile is inverted, with a 
maximum near 5-10 km altitude. Thus the near-surface atmosphere absorbs more radiation 
relative to its emission than does the atmosphere at 10 km altitude. Furthermore, the larger 
optical depths near the surface do not allow the emitted thermal radiation to escape the layer as 
easily as at higher altitudes. 

The maximum cooling due to dust below 30 km altitude occurs in the 12-18 pm wave- 
length interval. This is due to the cold temperatures at winter polar latitudes. The overlap of 
strong dust cooling in the 12-18 pm interval with the strong cooling by CO 2 in the 15 Jim band 
(12-19 p.m) is particularly important, as discussed by Lindner (1993b). Hence, to properly 
account for both dust and CO 2 cooling, they must be treated simultaneously (Lindner, 1993b), 
unlike what is usually done. Significant cooling due to dust also occurs at wavelengths longer 
than 18 pm. Cooling by dust at wavelengths shorter than 12 pm is inefficient due to the cold 
winter polar temperatures. (Recall that the peak Planck emission occurs at longer wavelengths 
for colder temperatures.) 

Cooling due to CO 2 near the surface occurs mostly in the wings of the 15 pm band, those 
parts of the 15 pm band where neither absorption nor emission is efficient. Photons emitted near 
the surface in the center of the 15 pm band (the strongly absorbing parts of the band) are rapidly 
re-absorbed. Hence, the cooling is smaller at lower altitudes (see Fig. 1) due to the inability of 
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photons to escape and cool the atmosphere. Photons in the line center can escape to space more 
easily at higher altitudes, explaining the higher cooling rates there. Cooling to the surface also 
occurs, and is included in all calculations. However, in order to cool to the surface, photons must 
pass through large optical depths. Cooling to the surface is only important right near the surface. 

Cooling rates in the 4.3 pm band of CO 2 increase from 10 ' 6 K/Mars day at the surface to 
3X10 -4 K/Mars day at 30 km. Hence, cooling in the 4.3 (im band of CO 2 is about 10*3% of the 
total cooling via CO 2 . Clearly, 4.3 pm band cooling is not an important process in the winter 
polar atmosphere of Mars. The 4.3 pm cooling rate has the same altitude behavior as the 15 pm 
cooling rate. As with the 15 pm band, the line center is optically thick near the surface and all 
emitted photons are quickly reabsorbed, preventing effective line-center cooling in the 4.3 pm 
band. 4.3 pm cooling would be larger at warmer latitudes, as thermal emission would shift to 
shorter wavelengths. However, 4.3 pm band cooling will never be an important cooling source 
in the martian atmosphere. The 4.3 pm band and other NIR bands are important cooling 
processes in the Earth's atmosphere. 

The heating and cooling rates for a late winter (L s = 343°) atmosphere at 51 ° N latitude 
with more dust (t v = 0.5) are shown in Fig. 5. Comparing Fig. 1 and Fig. 5, we see that dust 
heating and cooling rates increase at higher dust opacities, at all altitudes. Clearly, dust heating 
and cooling dominates over that of gas, except possibly for CO 2 cooling above 30 km altitude. 

As Ty = 0.5 was not unusual during the Viking mission, dust heating and cooling would dominate 
for most winter polar latitudes and seasons. Dust heating and cooling at larger dust loadings is 
even more dominant. During global dust storms ( t v ~3 ), heating and cooling by O 3 and CO 2 will 
be negligible compared to that of dust. However, dynamical transport of heat increases during 
dust storms. 

For t v (dust) = 0.5, meridional and vertical heat transport may be even less important than 
for the t v = 0.2 case. Based on observational and modeling evidence, meridional and vertical 
winds do not change much between the 0.2 and 0.5 cases (Haberle et al., 1982). Therefore, 
radiative equilibrium may be a more valid assumption at t v = 0.5 because total radiative heating 
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and cooling rates are twice as large for t v (dust) = 0.5, and are greater than 3 K/Mars day at 
almost all altitudes. 

Heating and cooling rates are well balanced in the lower 20 km of the atmosphere in Fig. 

5. Above 20 km the assumed temperatures are incorrect Temperatures closer to 140 K at 30 km 
altitude would give a better balance between the heating and cooling rates above 20 km. Indeed, 
the inability of dust cooling to keep up with dust heating at higher dust loadings in addition to 
heat transport is the explanation for the higher atmospheric temperatures observed during global 
dust storms. The explanation lies in the negative feedback of larger dust loadings, in that larger 
opacities decrease the ability of photons to escape. CO 2 cooling also becomes less effective at 
higher dust loadings due to dust-gas interaction (Lindner, 1993b). 

Obviously, higher dust loading results in increased dust heating rates. The higher optical 
depth chokes off some light from reaching lower altitudes, which explains why the increases in 
the heating rate for higher dust optical depths are not as large at the surface as at higher altitudes. 
Cooling rates for 0.5 vertical optical depths of dust are twice the cooling rates for x v = 0.2. 

While the optical depth is 2.5 times as large, cooling rates are only 2.2 times larger. This is 
because two negative feedbacks exist in that larger optical depths also hinder the ease of escape 
for emitted photons, and in that larger optical depths increase the thermal flux which increases 
absorption and heating. The thermal heating at altitudes above 25 km also increases. The 
heating above 25 km increases because the upward thermal flux is larger for t v = 0.5, which 
results in larger absorption in the upper atmosphere, while the thermal emission above 25 km 
remains the same. 

The phenomenon of thermal heating is also partly responsible for the low lapse rates in the 
martian atmosphere. Because dust and CO 2 are both radiatively active in the infrared, the atmo- 
sphere near the surface is able to cool very effectively, and keep near-surface temperatures low. 
But the high thermal fluxes are also causing a heating as they are absorbed by the other regions 
of the atmosphere. Any part of the atmosphere that is too cold will be heated by both solar and 
thermal flux. Higher dust loading increases both solar and thermal heating more effectively than 
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thermal cooling. This results in more isothermal conditions as the dust loading increases, as is 
observed. 

Figure 6 shows the heating and cooling rates deeper in the winter polar region at 70°N lati- 
tude for the same season for normal, low-dust conditions. The cooling rates are the same as at 
57°N latitude (Figure 1), because the same temperature profile is used. The same temperatures 
are used at both latitudes to illustrate latitude-dependent changes in the heating rates. The lower 
solar fluxes at 70"N latitude (due to the larger solar zenith angle) result in lower heating rates. 

The atmospheric heating and cooling is approximately equal near the surface and above 20 km. 
However, the assumed temperature inversion is too strong for radiative equilibrium for 70°N 
latitude conditions, as the cooling at 10 km is twice the heating. Slightly lower temperatures in 
the 5-10 km range would provide a better balance between heating and cooling, and would agree 
with observational evidence (Lindal et al., 1979). Dynamics could be relatively more important 
at transporting heat at 70°N latitude because the heating rates are lower. However, observations 
and dynamical modeling indicate that the atmosphere is even more stable against motion at these 
latitudes which would lower heat transport (Haberle et al., 1979). 

The excess cooling could also go into condensation, rather than in changing the tempera- 
ture (Pollack et al., 1990). Indeed, optically thick clouds are frequently observed in the 5-10 km 
altitude range. As the atmospheric temperature near the surface is already at the CO 2 condensa- 
tion temperature, the clouds at this latitude will be at least partly composed of CO 2 ice. Clouds 
would also alter the heating rate profile by shifting the location of solar flux through scattering. 

Comparing Fig. 1 and Fig. 6 shows that the relative importance of O 3 to CO 2 heating is 
virtually the same at 70°N latitude as it is at 57°N latitude. However, both O 3 and CO 2 heating 
are less important relative to dust. CO 2 heating occurs mostly in bands which are optically thick. 
At higher latitudes, the larger solar zenith angles only serve to decrease the transmission in these 
already optically- thick bands and hence decrease their relative importance. O 3 is less important 
to the heat budget at 70°N latitude due to the lower O 3 abundances there (Barth et al., 1973; Lane 
et al., 1973; Lindner, 1991b). The same general altitude and wavelength behavior in heating and 
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cooling is seen at 70°N latitude as at 57°N latitude. Therefore, while dust heating is less at higher 
latitudes due to the decreased solar flux, the importance of dust heating relative to gas heating 
has increased. 

The large solar zenith angle 0 at 70°N latitude increases the effective optical depth (x = 
tv/cos 0) of dust, which increases the absorption of solar flux, and hence increases the heating. 
The larger solar zenith angle also decreases the solar flux, and hence decreases the heating. In an 
optically-thin medium, these effects would cancel and the heating would be the same at 70'N 
latitude as at 57°N latitude. However, the 0.2 vertical optical depths of dust yield an effective 
optical depth of 1.7 at 70°N latitude. Thus, the dust is actually choking off the solar-flux, which 
decreases the heating rate relative to 57°N latitude. The lower heating at the surface is due to the 
large optical-depth of dust. 

CO2 NIR heating at 70° N latitude is half that at 57 °N latitude. The general behavior of 
each band at 70°N latitude is similar to that at 57°N latitude, although the bands saturate at higher 
altitudes. Consequently, the 4.3 and 15 |JLm bands are less important at 70°N latitude, and the 2.0 
and 2.7 pm bands are less important near the surface. The minor bands (1.3, 1.4, 1.6, 4.8 and 5.2 
pm) are optically thin at 70°N latitude, as at 57°N latitude. The heating in all bands is reduced at 
70°N latitude, due to the lower solar flux (via the larger solar zenith angle). Comparing the 
results of 57°N latitude and 70'N latitude, it is apparent that CO2 NIR heating would be 
markedly higher at equatorial latitudes, due to the increased solar flux, and the larger contribu- 
tions by the saturated 2.0, 2.7, 4.3 and 15 |im bands. Hence, the relative importance of dust and 
CO2 will shift in the favor of CO2 at more equatorward latitudes. O3 heating will not be impor- 
tant at equatorial and mid-latitudes because negligible ozone abundances exist there (Barth et al., 
1973). Total heating rates will increase with decreasing latitude due to the decrease in the solar 
zenith angle. The net result is higher atmospheric temperatures with decreasing latitude, which is 
in fact observed. 

SUMMARY AND CONCLUSIONS 

Heating in the winter polar atmosphere of Mars is provided mostly by dust at visible 
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wavelengths, especially for high dust loading, or at high latitudes. CO 2 NIR heating is always 
less than dust heating at 57°N latitude during late winter, but is comparable to dust heating at low 
altitudes for minimal dust loading. CO 2 NIR heating is more important at equatorial latitudes, 
and less important at high latitudes. Most CO 2 heating comes from the 2.7 |im band above 10-20 
km altitude, with most heating by the 2.0 |im band below. 1.3 Jim band heating is appreciable 
near the surface. The importance of minor CO 2 bands requires their inclusion in models of polar 
winter winds and surface energy balance. Ozone heating is only 10% of the total heating at 57°N 
latitude, and is even less at other latitudes due to lower ozone abundances. Ozone heating was 
suggested to be important to the polar heat budget by Kuhn et al. (1979), but the importance of 
dust heating was ignored by Kuhn et al. Heating by CO 2 , O 2 , H 2 O, HO 2 and H 2 O 2 is negligible 
at ultraviolet wavelengths below 30 km altitude. 

CO 2 15 Jim band cooling is the dominant source of cooling at high altitudes for low dust 
abundances, but is ineffective near the surface. CO 2 4.3 Jim band cooling is negligible. Dust 
cooling is the dominant source of cooling at winter polar latitudes under most conditions, with 
the largest dust cooling in the 12-18 Jim wavelength range. Dust cooling increases at higher dust 
loadings; however dust cooling does not increase as fast as dust heating, due to several negative 
feedbacks. As a result, atmospheric temperatures rise with increasing dust opacity, in agreement 
with observations. The warmest winter polar temperatures occur during global dust storms, 
when the largest dust opacity exists. However, whether the high temperatures at polar latitudes 
during global dust storms are due primarily to radiative processes or dynamical heat transport is 
uncertain. Radiative effects of dust have little effect on the overall recession of the polar cap 
(Lindner, 1990). 

Radiative processes are responsible for the low lapse rates in the martian atmosphere. 
Significant dust and gas heating occurs at all altitudes, damping out inhomogeneities in tempera- 
ture. Any region of the atmosphere which is significantly colder than the rest of the atmosphere 
is warmed not only by solar flux but also by the absorption of the large thermal flux. Large lapse 
rates are quickly eliminated by solar and infrared heating. As atmospheric dust loading 
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increases, the atmosphere becomes more isothermal due to the increased solar and thermal heat- 
ing, particularly at higher altitudes. Ignoring dust altogether leads to the opposite situation. CO 2 
cools ineffectively near the surface, but cools readily at high altitudes, which leads to stronger 
lapse rates in the dust-free atmosphere. 
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Figure Captions 


Figure 1. 

Figure 2. 

Figure 3. 

Figure 4. 

Figure 5. 
Figure 6. 


Total atmospheric radiative heating and cooling as a function of altitude on Mars 
at 57°N latitude with x v (dust) = 0.2. Also shown is the contribution by CO2, O3 
and dust. 

Total solar heating rates for CO2 as a function of altitude and the contribution by 
each CO2 band (labeled by wavelength in |im) at 57 °N latitude (a) and 70°N 
latitude (b). 

Total radiative heating rate as a function of altitude due to absorption of solar 
radiation by dust (x v = 0.2) and the contribution from each wavelength bin from 0 
to 5 microns (bins labeled on figure in Jim). 

Total infrared radiative cooling rate by dust as a function of altitude for x v = 0.2, 
and the contribution from each wavelength interval from 5 to 40 p.m (wavelength 
range of each bin is labeled in |im). 

As in Figure 1 except for x v (dust) = 0.5. 

As in Figure 1 except at 70°N latitude for x v (dust) = 0.2 
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ABSTRACT : The efficacy of the UV spectroscopy technique used by Mariner 9 to 
remotely measure ozone abundance at Mars is discussed. Previously-inferred ozone abundances 
could be underestimated by as much as a factor of 3, and much of the observed variability in the 
ozone abundance could be due to temporal and spatial variability in cloud and dust amount 
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INTRODUCTION 

Ozone is a key to understanding atmospheric chemistry on Mars. Over 20 photochemical 
models of the martian atmosphere have been published, and O3 is often used as a benchmark for 
these models (e.g., Lindner, 1988; Shimazaki, 1989; Krasnopolsky, 1993). O3 abundance has 
been inferred from instrumentation on several spacecraft, with the most complete coverage pro- 
vided by Mariner 9 (Barth et al., 1973; Lane et al., 1973; Wehrbein et al., 1979). The Mariner 9 
UV spectrometer scanned from 2100 to 3500 Angstroms with a spectral resolution of 15 
Angstroms and an effective field-of-view of approximately 300 km 2 . The only atmospheric 
absorption in the 2000 to 3000 A wavelength region was previously assumed to come from the 
Hartley band system of ozone (Barth et al., 1973; Lane et al. 1973), which has an opacity of order 
unity at winter polar latitudes (Lindner, 1988). Therefore the amount of ozone was inferred by 
fitting this absorption feature with laboratory data of ozone absorption, as shown in Fig. 1 (Lane et 
al., 1973). Mars O3 shows strong seasonal and latitudinal variation, with column abundances 
ranging from 0.2 |im-atm at equatorial latitudes to 60 jim-atm over northern winter polar latitudes 
(1 itm-atm is a column abundance of 2.689xl0 15 molecules cm -2 ). However, the O3 abundance is 
never great enough to significantly affect atmospheric temperatures (Lindner, 1991; 1993a) or sur- 
face temperatures and frost amounts (Lindner, 1990; 1992; 1993b). Figure 2 (Barth, 1985) shows 
some of the previously-inferred O3 abundances. 

MODELING PROCEDURE 

A radiative transfer computer model is used to re-examine the Mariner 9 UV spectra. The 
discrete ordinate method of Stamnes et al. (1988) is used to treat the scattering and absorption of 
monochromatic radiation through the martian atmosphere. O3 absorption cross-sections from 2000 
A to 3200 A are included (Daumont et al., 1983; Freeman et al, 1984). Ultraviolet absorption 
cross-sections for CO2, O2. H2O, HO2 and H2O2 are negligible (e.g., Lindner, 1988; 

1991; 1993b). The Rayleigh scattering optical depth is computed as in Hansen and Travis (1974), 
using parameters appropriate for Mars. Atmospheric composition is taken as 95% CO2 (Owen et 
al., 1977). The altitude dependence of O3 is based on model results (Lindner, 1988). 
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Atmospheric properties are zonally averaged and assumed azimuthaUy-independent. The 
region from the surface to 40 km altitude is broken into twenty 2-km-thick layers to account for 
vertical inhomogeneity. The Chapman function is used to approximate the slant path in place of the 
secant function (Smith and Smith, 1972), because the winter polar atmosphere always has large 
solar zenith angles, and the secant function is in error for large angles. This work examines 57°N 
latitude in late winter (L$ = 343*), which is when the maximum O3 column abundance was 
observed (Barth et al, 1973; Lane et al., 1973). [Ls, the solar longitude, is a seasonal index; L$ of 
0°, 90°, 180°, 270°, correspond to northern spring equinox, summer solstice, autumnal equinox, 
and winter solstice, respectively]. The surface is covered by ice with an albedo of 0.6 (Warren and 
Wiscombe, 1980; Hapke et al., 1981; James and Lumme, 1982; Warren et al., 1990; Lindner, 
1993a). The altitude profile of temperature rises linearly from 150K at the surface to 160K at 10 
km, and then falls linearly with altitude to 130 K at 40 km, typical for winter polar conditions 
(Lindal et al., 1979; Kieffer, 1979). 

Dust opacities varied from 0.2 to 1.0 for conditions other than global dust storms during 
the Mariner 9 and Viking observations (Briggs and Leovy, 1974; Pollack et al., 1979; Lumme and 
James, 1984), although dust opacities over winter polar latitudes may be slightly less [e.g., 
Lindner, 1990]. Dust opacity varies with wavelength (Zurek, 1978; 1982). The vertical 
distribution of dust is well-mixed to 20 km altitude (Anderson and Leovy, 1978; Zurek, 1982; 
Korablev et al., 1993). The single scattering albedo of airborne dust has a solar average of 0.9 
(Clancy and Lee, 1991), but is only 0.6 in the 2000 A to 3000 A range (Pang et al., 1976; 

Thorpe, 1977). Scattering of radiation by dust is represented by the Henyey-Greenstein phase 
function, with an asymmetry parameter of 0.55 (Clancy and Lee, 1991). Typical cloud and fog 
opacities of 1,0 are taken from theory (Kulikov and Rykhletskii, 1983) and observations (Moroz, 
1976; Pollack et al., 1977; Clancy and Lee, 1991), although large seasonal and latitudinal 
variations exist (Briggs and Leovy, 1974). A single-scattering albedo for cloud particles of 1.0 is 
taken from observations (Clancy and Lee, 1991). The Henyey-Greenstein phase function with 
asymmetry parameter of 0.55 is also used to describe the scattering by clouds (Clancy and Lee, 



1991). Computational difficulties which accompany highly asymmetric phase functions are 
removed with the Delta-L method (Wiscombe, 1977). 

RESULTS AND DISCUSSION 

Using a constant mixing ratio for O3 (e.g., Wehrbein, 1979; Lindner, 1988) and no 
chemical or radiative interaction between O3 and clouds/dust. Fig. 3 shows that when typical 
amounts of dust and cloud are present that significant underestimation of O3 abundance occurs. A 
factor of 3 times as much O3 is needed to generate the same spectrum the spacecraft would measure 
for a cloudy, dusty atmosphere as for a clear atmosphere. If the scattering properties of martian 
clouds and dust were well known, then their appearance would not be a problem, as a model 
would be capable of retrieving the O3 abundance. However, these properties arc not well known, 
which raises doubts about the effectiveness of the UV reflectance spectroscopy technique for 
measuring O3 abundance on Mars. The simulations shown in Fig. 3 arc repeated for a range in 
solar zenith angle (50’-90*), ground albedo (0.3-0.8), altitude distribution of O3, satellite viewing 
geometries, dust scattering properties, and cloud, dust and O3 abundances. A factor of 3 
underestimation is typical, with greater underestimation for high ground albedo or high dust opaci- 
ties. Even if scattering by clouds is properly accounted for (as previously done with Mariner 9 
data reduction in [4]), masking by dust can easily result in factor of 2 underestimation. Results are 
not strongly dependent on solar zenith angle. 

Spatial and temporal variability in temperature and water vapor have been claimed to ac- 
count for the scatter of the data points in Fig. 2 (Barth and Dick, 1974). A decrease in temperature 
results in a decrease in water vapor, assuming the atmosphere is saturated. A decreased water 
vapor abundance decreases the availability of odd hydrogen (H, OH, and HO2), which converts 
CO and O into CO2 catalytically, decreasing the abundance of O needed to form O3. However, 
water vapor is a small source of odd hydrogen in the winter polar atmosphere compared to H2, and 
may not account for most of the variability in Fig. 2 (Lindner, 1988). Masking by clouds and dust 
may also account for some of the observed O3 variability, because the nature and opacity of the 
clouds and dust at winter polar latitudes change significantly spatially and temporally. As the 
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maximum O3 abundance resides near the surface (Lindner, 1988), spacecraft must be able to 
observe through the entire cloud and dust abundance in order to measure the total O3 column 
abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and the airborne 
dust must be traversed twice; first by the incoming solar flux down to the surface, and then once 
again upon reflection from the surface out to the spacecraft In addition, the large solar zenith 
angles at winter polar latitudes mean several times the vertical opacity of cloud and dust must be 
traversed. Indeed, part of the observed latitudinal variation in O3 abundance in Fig. 2 may be due 
to the inability of the spacecraft to observe through the increasing effective optical depths as one 
goes poleward. 

By using a photochemical model which included multiple scattering of solar radiation, 
Lindner (1988) showed that the absorption and scattering of solar radiation by clouds and dust 
should actually increase O3 abundances at winter polar latitudes. Hence, regions with high dust 
and cloud abundance could contain high O3 abundances (heterogeneous chemistry effects have yet 
to be fully understood [Atreya and Blamont, 1990; Krasnopolsky, 1993]). It is quite possible that 
the maximum O3 column abundance observed by Mariner 9 of 60pm-atm is common. In fact, 
larger quantities may exist in some of the colder areas with optically thick clouds and dust As the 
Viking period often had more atmospheric dust loading than did that of Mariner 9, the reflectance 
spectroscopy technique may even have been incapable of detecting the entire O3 column abundance 
during much of the Mars year that Viking observed, particularly at high latitudes. The behavior of 
O3 is virtually unknown during global dust storms, in polar night, and within the polar hood, 
leaving large gaps in our understanding. 

Other possibilities for measuring O3 abundance include solar occultation (Blamont et al., 
1989), IR observations in the 9.6 pm O3 absorption band (Espanek et al., 1991), and observations 
of the O2 dayglow at 1.27 pm, produced by photolysis of O3 (Traub et al., 1979). However, 
further studies of these other techniques are required, especially as regards the effects of clouds 
and dust. 
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FIGURE CAPTIONS 


Figure 1. Ultraviolet spectrum measured by Mariner 9 at 57*N latitude on orbit 144 (taken from 
Lane et al., 1973). To enhance the O 3 absorption feature, this spectrum was divided by one 
obtained at 20*N latitude on orbit 144, where O 3 abundances are minimal. 

Figure 2. Measurements of the O 3 column abundance previously inferred from the Mariner 9 UV 
spectrometer data during the northern winter, L$ = 330-360*, in the northern hemisphere (taken 
from Barth, 1985). 

Figure 3. Synthetic spectra as would be observed by spacecraft for atmospheres with no cloud or 
dust and 30 pm-atm O 3 (solid line), vertical opacities of dust and cloud of 0.3 and 1.0, respec- 
tively, and 30 pm-atm of O 3 (dashed line), and vertical opacities of dust and cloud of 0.3 and 1.0, 
respectively, and 100 pm- atm of O 3 (dotted line). All cases assume a solar a zenith angle of 75* 
(typical for winter polar observations), and a polar cap albedo of 0 . 6 . 
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Rubincam (1992) has splendidly shown that the martian obliquity is dependent on 
the seasonal polar caps. In particular, Rubincam analytically derived this dependence, 
and showed that the change in obliquity is directly proportional to the seasonal polar cap 
mass. Specifically, Rubincam showed 

dT/dt = 3 x 10*10 M(t)/M(0) degrees / Earth year (1) 

where T is the obliquity and M is the mass of the seasonal polar caps, with time t of 0 
being the present. This expression assumes uniformly thick spherical caps with identical 
angular radii of 45°. However even if a very different polar cap mass distribution is used, 
Rubincam estimates the total uncertainty in the constant in Equation (1) to be less than a 
factor of 2. Using the current mass of the seasonal polar cap as typical over geologic 
time, Rubincam calculates that the amount that the obliquity would secularly change is 
only 1.4°. Considering that the current obliquity of Mars is 25°, Rubincam concludes that 
seasonal friction does not appear to have changed Mars' climate significantly. 

Using a computer model for the evolution of the martian atmosphere, Haberle et 
al. (1992a, b) have made a convincing case for the possibility of huge polar caps, about 
ten times the mass of the current polar caps, which exist for a significant fraction of the 
planet's history. Given the large uncertainties in input parameters and in the model itself, 
the results must be regarded as speculative. Also, the Haberle et al. results have been 
unable to favor or rule out a large polar cap scenario versus a small polar cap scenario. 

Nonetheless, even just the possibility of massive polar caps makes the Rubincam 
paper more noteworthy. Since Rubincam showed that the effect of seasonal friction on 
obliquity is directly proportional to polar cap mass, a scenario with a ten-fold increase in 
polar cap mass over a significant fraction of the planet's history would result in a secular 
increase in Mars' obliquity of perhaps 10° (using Equation 1). Hence, the Rubincam 
conclusion of an insignificant contribution to Mars' climate by seasonal friction may be 
incorrect. Furthermore, if seasonal friction is an important consideration in the obliquity 
of Mars, this would significantly alter the predictions of past obliquity as presented by 



Ward (1973; 1974; 1979), Murray et al. (1973), Ward et al. (1979), Rubincam (1990), 
Chao and Rubincam (1990), Bills (1990), Ward and Ruby (1991), Touma and Wisdom 
(1993), and Laskar and Robutel (1993). That in turn would significantly alter the 
predictions of past climate which are based on obliquity predictions (Sagan et al., 1973; 
Ward et al., 1974; Toon et al., 1980; Fanale et al., 1982; Pollack and Toon, 1982; 
Francois et al., 1990). The mechanics of the polar cap system also depend on obliquity 
(Leighton and Murray, 1966; James and North, 1982; Lindner, 1990; 1992; 1993; Wood 
and Paige, 1992.) If obliquities were often much smaller than at present, that could have 
implications for past atmospheric composition (Lindner and Jakosky, 1985). 

Given the enormity of the implications, the work of Rubincam should be given 
more attention and study. Perhaps further modeling of obliquity could be used to rule 
out the possibility of large polar caps for extended times, which would assist modeling of 
atmospheric evolution. Similarly, modeling of atmospheric evolution should be given 
more attention and study because of the implications for obliquity history, and therefore 
climate history. 
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WHY IS THE NORTH POLAR CAP ON MARS DIFFERENT THAN THE SOUTH POLAR CAP? 

Bernhard Lee Lindner . Atmospheric and Environmental Research, Inc. 

840 Memorial Drive, Cambridge, Mass. 02139-3794, USA 




Introduction . One of the most puzzling mysteries about the planet Mars is 
the hemispherical asymmetry in the polar caps. Every spring the seasonal polar 
cap of CO 2 recedes until the end of summer, when only a small part, the residu- 
al polar cap, remains. During the year that Viking observed Mars, the residual 
polar cap was composed of water ice in the northern hemisphere [Kieffer et al.. 
Science . 194 . 1341, 1976] but was primarily carbon dioxide ice in the southern 
hemisphere (Kieffer, J . Geophvs . Res . . 84 . 8263, 1979]. Scientists have sought 
to explain this asymmetry by modeling observations of the latitudinal recession 
of the polar cap and seasonal variations in atmospheric pressure (since the 
seasonal polar caps are primarily frozen atmosphere, they are directly related 
to changes in atmospheric mass) . These models reproduce most aspects of the 
observed annual variation in atmospheric pressure fairly accurately. Further- 
more, the predicted latitudinal recession of the northern polar cap in the 
spring agrees well with observations, including the fact that the CC >2 ice is 
predicted to completely sublime away. However, these models all predict that 
the carbon dioxide ice will also sublime away during the summer in the southern 
hemisphere, unlike what is observed. This paper will show how the radiative 
effects of ozone, clouds, and airborne dust, light penetration into and through 
the polar cap, and the dependence of albedo on solar zenith angle affect CO 2 
ice formation and sublimation, and how they help explain the hemispherical 
asymmetry in the residual polar caps. These effects have not been studied with 
prior polar cap models. 


Ozone. Clouds, and Airborne Dust . Since O 3 is more prevalent in the 
northern hemisphere than in the southern hemisphere, O 3 was suggested as a 
cause for the hemispherical asymmetry in the residual polar caps by Kuhn et al. 
( J . Geophvs , Res . . 84 . 8341, 1979). However, Lindner (submitted to Icarus . 
1991) has shown that O 3 has a minor effect on the atmospheric temperature, and 
hence on the infrared radiation which strikes the polar cap, and Lindner ( J . 
Geophys , Res , . 95, 1367, 1990) has shown that O 3 absorbs less than 1% of the 
total solar radiation absorbed by the polar cap. Thus, O 3 is not an important 
consideration in the polar cap energy budget. 

Lindner (1990) has computed the solar and thermal flux striking the polar 
cap of Mars for various ozone, dust, and cloud abundances and for three solar 
zenith angles. These calculations have been inserted in the polar-cap models 



of Lindner ( Eos Trans. AGU . 67 . 1078, 1986) and Jakosky and Haberle ( J. Geo- 
nhvs . Res . . 95, 1359, 1990). Vertical optical depths of dust and cloud ranging 
from zero to 1 cause little change in the total flux absorbed by the polar cap 
near its edge but increase the absorbed flux significantly as one travels pole- 
ward. Observed hemispherical asymmetries in dust abundance, cloud cover, and 
surface pressure combine to cause a significant hemispherical asymmetry in the 
total flux absorbed by the residual polar caps, which helps to explain the 
dichotomy in the residual polar caps on Mars. 

Light Penetration . Penetration of solar radiation into the cap itself is 
included in my polar cap model, based on the theoretical work of Clow ( Icarus . 
72 . 95, 1987). I find that the inclusion of light penetration slightly 
decreases the albedo needed in the model to keep CC^-ice year-round at the 
south pole by on the order of 1%. The required albedo is decreased because 
some solar radiation is used to heat the subsurface, and not all of this heat 
is transported back to the surface. Overall, I conclude that penetration of 
light into the polar cap has only a small effect ,on the polarcap energy budget. 

Albedo and the Solar Zenith Angle . Warren et al. ( J. Geophvs. Res. . 95 . 
14717, 1990) has computed the dependence of the albedo of the martian polar 
caps on solar zenith angle, and these calculations have been included in my po- 
lar cap model. Since the albedo of ice increases and becomes more forward 
scattering at higher solar zenith angles, and since the solar zenith angle be- 
comes higher as one approaches the pole, the albedo is greatest at the pole. 
This decreases absorption of sunlight, hence increasing survivability of C (>2 
ice. In fact, this Increases the survivability of ice enough to offset the de- 
crease in survivability of ice due to the radiative effects of clouds and dust. 

Discussion. The combination of the effects of solar zenith angle on albe- 
do and the radiative effects of clouds and dust act to extend the lifetime of 
CO 2 ice on the south pole relatively more than on the north pole, explaining 
the hemispherical asymmetry in the residual polar caps without the need of a 
hemispherical asymmetry in polar cap albedo. Another positive aspect this 
solution is that neither the inclusion of solar zenith angle effects on ice 
albedo nor the radiative effects of clouds and dust should appreciably change 
model predictions of the annual cycle of pressure or polar cap recession equa- 
torward of 75” latitude, since approximately 90% of the seasonal CO 2 frost is 
equatorward of 80" latitude. Hence, the good model agreement noted by prior 
researchers to the seasonal cycle in atmospheric pressure and to the recession 
of the polar cap equatorward of 80” latitude is retained. 



MARS SEASONAL C0 2 -ICE LIFETIMES AND THE ANGULAR DEPENDENCE OF ALBEDO 
Bernhard Lee Lindner . AER, 840 Memorial Drive, Cambridge MA 02139 USA 

The albedo of the polar caps on Mars brightens appreciably at high 
solar zenith angle (Warren et al., J , Geonhvs . Res . . 95 . 14717, 1990), 
an effect not included in prior polar-cap energy-balance models. This 
decreases absorption of sunlight by the polar cap, hence decreasing 
sublimation of C0 2 ice. Lindner ( J . Geonhvs . Res . . 95 . 1367, 1990) 
has shown that the radiative effects of clouds and airborne dust will 
increase sublimation of C0 2 ice over that predicted by prior polar-cap 
energy-balance models . Furthermore , . observations hint that more 
clouds may exist in the northern hemisphere, which Lindner (1990) has 
shown would sublime C0 2 ice more quickly in the north than in the 
south. I show here that the effects of the solar zenith angle depen- 
dence of albedo and the radiative effects of clouds and dust offset 
each other, but act to extend the lifetime of C0 2 ice on the south 
pole more than on the north pole, possibly explaining the observed 
hemispherical asymmetry in the residual polar caps without the need of 
a hemispherical asymmetry in polar- cap albedo required by prior 
models. Another positive aspect of this solution is that neither the 
inclusion of the solar zenith angle dependence of albedo nor the radi- 
ative effects of clouds and dust should appreciably change prior model 
agreement with observations of the annual cycle of surface pressure 
and the recession of the polar caps equatorward of 75° latitude. 
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SIMULATIONS OF THE SEASONAL POLAR CAPS ON MARS Bernhard Lee Lindner . 
Atmospheric and Environmental Research, Inc., 840 Memorial Drive, Cambridge, 
Mass. 02139-3794, USA 

Introduct ion . One of the most puzzling mysteries about the planet Mars 
is the hemispherical asymmetry in the polar caps. Every spring the seasonal 
polar cap of CO 2 recedes until the end of summer, when only a small part, the 
residual polar cap, remains. During the year that Viking observed Mars, the 
residual polar cap was composed of water ice in the northern hemisphere [Kief- 
fer et al. , Science . 194 . 1341, 1976] but was primarily carbon dioxide ice in 
the southern hemisphere [Kieffer, J. Geophvs . Res. . 84, 8263, 1979], Scien- 
tists have sought to explain this asymmetry by modeling observations of the 
latitudinal recession of the polar cap and seasonal variations in atmospheric 
pressure (since the seasonal polar caps are primarily frozen atmosphere, they 
are directly related to changes in atmospheric mass). These models reproduce 
most aspects of the observed annual variation in atmospheric pressure fairly 
accurately. Furthermore, the predicted latitudinal recession of the northern 
polar cap in the spring agrees well with observations, including the fact that 
the CO 2 ice is predicted to completely sublime away. However, these models 
all predict that the carbon dioxide ice will also sublime away during the sum- 
mer in the southern hemisphere, unlike what is observed. This paper will show 
how the radiative effects of ozone, clouds, and airborne dust, light penetra- 
tion into and through the polar cap, and the dependence of albedo on solar ze- 
nith angle affect CO 2 ice formation and sublimation, and how they help explain 
the hemispherical asymmetry in the residual polar caps. These effects have 
not been studied with prior polar cap models. 

Ozone. Clouds, and Airborne Dust . Since O 3 is more prevalent in the 
northern hemisphere than in the southern hemisphere, O 3 was suggested as a 
cause for the hemispherical asymmetry in the residual polar caps by Kuhn et 
al. ( J. Geophvs. Res. . 84 . 8341, 1979). However, Lindner (submitted to Icar- 
us . 1991) has shown that O 3 has a minor effect on the atmospheric temperature, 
and hence on the infrared radiation which strikes the polar cap, and Lindner 
( J . Geophvs . Res . . 95, 1367, 1990) has shown that O 3 absorbs less than 1% of 
the total solar radiation absorbed by the polar cap. Thus, O 3 is not an im- 
portant consideration in the polar cap energy budget. 

Lindner (1990) has computed the solar and thermal flux striking the polar 
cap of Mars for various ozone, dust, and cloud abundances and for three solar 
zenith angles. These calculations have been inserted in the polar-cap model 
of Lindner ( Eos Trans. AGU . 67, 1078, 1986). Vertical optical depths of dust 
and cloud ranging from zero to 1 cause little change in the total flux absorb- 
ed by the polar cap near its edge but increase the absorbed flux significantly 
as one travels poleward. Observations hint that hemispherical asymmetries in 
dust abundance and cloud cover exist, and these would combine to cause a sig- 
nificant hemispherical asymmetry in the total flux absorbed by the residual 
polar caps, which helps to explain the dichotomy in the residual polar caps. 

Light Penetration . Penetration of solar radiation into the cap itself is 
included in the polar cap model of Jakosky and Haberle ( J . Geophvs , Res . . 95 . 
1359, 1990), based on the theoretical work of Clow (Icarus, 72 . 95, 1987). 
Lindner and Jakosky ( B . A . A. S , . 22., 1060, 1990) find that the inclusion of 
light penetration slightly decreases the albedo needed in the model to keep 
C 02 *ice year-round at the south pole by on the order of 1%. The required al- 
bedo is decreased because some solar radiation is used to heat the subsurface, 
and not all of this heat is transported back to the surface. Overall, we con- 
clude that penetration of light into the polar cap has only a small effect on 
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Che polarcap energy budget. 

Albedo and the Solar Zenith Angle . Warren et al. ( J- Ceophvs . Res . 95, 
14717, 1990) have computed the dependence of the albedo of the martian polar 
caps on solar zenith angle, and these calculations have been included in ray 
polar cap model. Since the albedo of ice increases and becomes more forward 
scattering at higher solar zenith angles, and since the solar zenith angle be- 
comes higher as one approaches the pole, the albedo is greatest at the pole. 
This decreases absorption of sunlight, hence increasing survivability of CO 2 
ice. In fact, this increases the survivability of ice enough to offset the 
decrease in survivability of ice due to radiative effects of clouds and dust. 

Discussion. The combination of the effects of solar zenith angle on al- 
bedo and the radiative effects of clouds and dust act to extend the lifetime 
of CO 2 ice on the south pole relatively more than on the north pole, possibly 
explaining the hemispherical asymmetry in the residual polar caps without the 
need of a hemispherical asymmetry in polar cap albedo. This does not imply 
that a hemispherical asymmetry in polar cap albedo does not exist, but that 
one is not necessary. 

Observations of the regression of the polar caps and the annual cycle in 
atmospheric pressure are reproduced fairly well by the model, as shown in the 
figures, although further improvement is needed. When CO 2 ice is retained at 
the south pole, the model predictions of the annual cycle in atmospheric pres- 
sure have a phase shift relative to the data, no matter what model input para- 
meters are used. We are investigating other processes not included in prior 
polar cap models. 





Comparison of the model to atmospheric pressure (upper right) , south polar cap 
regression (upper left) , and north polar cap regression (bottom) 
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DOES U V INSTRUMENTATION EFFECTIVELY MEASURE 
OZONE ABUNDANCE? Bernhard Lee Lindner. AER, 840 
Memorial Drive, Cambridge MA 02139, USA. 

Measurements of 0 3 on Mars provide significant information 
about the chemistry and composition of the atmosphere [1], includ- 
ing long-term changes [2], The most extensive and accurate data 
were inferred from the Mariner 9 UV spectrometer experiment, some 
of which are reproduced in Fig. 1 . Mars 0 3 shows strong seasonal and 
latitudinal variation, with column abundances ranging from 0.2 (im- 
atm at equatorial latitudes to 60 jxm-atm over the northern winter 
polar latitudes [ 1 ] ( 1 (im-atm is a column abundance of 2.689 x 1 0 15 
molecules cm -2 ). 

The Mariner9 UV spectrometer scanned from 2100 to 3500 A in 
one of its two spectral channels every 3 s with a spectral resolution 
of 15 A and an effective ficld-of-view of approximately 300 km 2 . 
Measurements weremade for almost half a martian year, with winter 
and spring in the northern hemisphere and summer and fall in the 
southern hemisphere. The detectability limit of the spectrometer was 
approximately 3 (im-atm of ozone. The process used by earlier 
investigators to extract the ozone abundance from the observed 
Mariner 9 spectra is as follows [1|. Each spectrum was filtered to 
remove spurious data points, then compared to the solar flux spec- 
trum and shifted slightly in wavelength in order to compensate for 
any systematic shift in the wavelength calibration of llic spectro- 
meter. Incoming solar radiation was assumed to undergo Rayleigh 
scattering by COj and Mic scattering by the polar hood, and to be 



wavelength (I) 

Fig. 2. Ultraviolet spectrum measured by Manner 9 at 57N latitude on orbit 
144. To enhance the Oj absorption feature, this spectrum was divided by one 
obtained at 20N latitude on orbit 144, where 0 3 abundances are minimal (1 ). 

reflected by a wavelength-independent surface albedo. The only 
atmospheric absorption in the 2000- to 3000-A region was assumed 
to come from the Hartley band system of ozone, and therefore the 
amount of ozone was inferred by fitting this absorption feature with 
laboratory data of ozone absorption, as shown in Fig. 2. 0 3 absorption 
of sunlight is not strong enough to affect atmospheric lemperatureon 
Mars (3), and hence cannot be inferred from temperature measure- 
ments. 

1 use a radiative transfer model based on the discrete ordinate 
method to calculate synthetic radiance spectra. Figure 3 shows that 
when typical amounts of dust and cloud arc present, significant 
underestimation of 0 3 occurs. A factor of 3 times as much 0 3 is 
needed to generate the same S|>cclruni as for a clear atmosphere. If 
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Fig. 3. Synthetic spectra as would be observed by spacecraft for atmospheres 
with (a) no cloud or dust and 30 mm-atm 0 3 , (b) vertical opacities of dust and 
cloud of 0.3 and 1.0 respectively and 30 mm-alm of Oj, and (c) vertical 
opacities of dust and cloud of 0.3 and 1 .0 respectively and 100 mm-atm of Oj. 
All cases assume a solar zenith angle of 75 (typical for winter polar observations) 
and viewing angle of 60. with azimuth angle of 0 (typical for Mariner 9). Polar 
cap albedo of 0.6. 


the scattering properties o f martian clouds and dust were w c l I know n 
Uteri their appearance would not be a problem, as a model would be 
capable of retrieving die 0 3 abundance. However, these properties 
arc not well known, which raises doubts about the effectiveness of U\c 
current UV spectroscopy technique used to measure 0 3 . 

Spatial and temporal variability in temperature and water vapor 
have been claimed to account for the scatter of Utc data points in Fig 
1 I'll. A decrease in temperature would result in a decrease in water 
vajior. if saturated as expected at prevalent temperatures. A de- 
creased water vapor abundance decreases the availability of odd 
hydrogen, which converts CO and O into C0 2 catalytically, decreas- 
ing the abundance of O needed to form 0 3 . However, water vapor is 
a small source of odd hydrogen in the winter polar atmosphere, and 
may not account for most of the variability in Fig. 1 [5]. Masking by 
clouds and dust may also account for some of the observed 0 3 
variability, because the nature and opacity of the clouds and dust in 
the polar hood change dramatically in latitude and even on a day-to- 
day basis. As the maximum 0 3 abundance resides near the surface 
(5J, spacecraft must be able to observe through the entire cloud and 
dust abundance in order to actually see the total 0 3 column abun- 
dance. If reflectance spectroscopy is used, as on Mariner 9. then the 
cloud and the airborne dust must be traversed twice. First by the 
incoming solar flux down to the surface and then once again upon 
reflection from the surface out to the spacecraft. In addition, the large 
solar zenith angles at winter polar latitudes mean several times the 
vertical opacity of cloud and dust must be traversed. Indeed, part of 
the observed latitudinal variation in 0 3 in Fig. 1 may be due to the 
inability of the spacecraft to observe through the increasing effective 
optical depths as one goes poleward. 

The UV spectrometer on Mariner 9 was incapable of penetrating 
the dust during dust storms [1]; the single-scattering albedo and 
phase function of airborne dust and cloud ice are not known to the 
degree required to extract the small UV signal reflected up from near 
the surface. The reflectance spectroscopy technique would also have 
difficulty detecting the total column abundance of 0 3 in cases where 
large dust abundances exist together with the polar hood, especially 
at high latitudes where large solar zenith angles magnify those 
optical depths; yet these eases would contain the maximum 0 3 based 
on theoretical results [5]. It is quite possible that the maximum 0 3 
column abundance observed by Mariner 9 of 60 pm-atm is common. 
In fact, larger quantities may exist in some of the colder areas with 
optically thick clouds and dust. As the Viking period often had more 
atmospheric dust loading than did that of Mariner 9. the reflectance 
spectroscopy technique may even have been incapable of detecting 
the entire 0 3 column abundance during much of the Mars year that 
Viking observed, particularly at high latitudes. 
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How well is martian ozone inferred with reflectance spectroscopy? 

Bernhard Lee Lindner 

AER, 840 Memorial Drive, Cambridge, MA 02139 (617)349-2280 

The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two 
spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field- 
of-view of approximately 300 km^. The only gaseous absorption in the 2000 to 3000 Angstrom 
region was assumed to come from the Hartley band system of ozone, and therefore the amount of 
ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption[l]. 
Mars O 3 as inferred from these spectra shows strong seasonal and latitudinal variation, with col- 
umn abundances ranging from 0.2 jim-atm at equatorial latitudes to 60 pm-atm over the northern 
winter polar latitudes [1]. The detectability limit of the spectrometer was approximately 3 |im-atm. 

I use a radiative transfer model based on the discrete ordinate method to calculate synthetic 
radiance spectra. When typical amounts of dust and cloud are present, significant underestimation 
of O 3 occurs. A factor of 3 times as much O 3 is needed to generate the same spectrum for cloudy, 
dusty atmospheres as for a clear atmosphere. If the scattering properties of martian clouds and 
dust were well known, then their appearance would not be a problem, as a model would be capable 
of retrieving the O 3 abundance. However, these properties are not well known, which raises 
doubts about the effectiveness of the current UV spectroscopy technique used to measure O 3 . 

Spatial and temporal variability in temperature and water vapor account have been claimed 
to account for the scatter of the data points [2], However, water vapor is a small source of odd 
hydrogen in the winter polar atmosphere, and may not account for most of the variability [3]. 
Masking by clouds and dust may also account for some of the observed O 3 variability, because the 
nature and opacity of the clouds and dust in the polar hood change dramatically in latitude and even 
on a day-to-day basis. As the maximum O 3 abundance resides near the surface [3], spacecraft 
must be able to observe through the entire cloud and dust abundance in order to actually see the to- 
tal O 3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and 
the airborne dust must be traversed twice; first by the incoming solar flux down to the surface, and 
then once again upon reflection from the surface out to the spacecraft. In addition, the large solar 
zenith angles at winter polar latitudes mean several times the vertical opacity of cloud and dust must 
be traversed; yet these cases would contain the maximum O 3 , based on theoretical results [3]. 
Indeed, part of the observed latitudinal variation in O 3 may be due to the inability of the spacecraft 
to observe through the increasing effective optical depths as one goes poleward. It is quite possi- 
ble that the maximum O 3 column abundance observed by Mariner 9 of 60jim-atm is common. In 
fact, larger quantities may exist in some of the colder areas with optically thick clouds and dust [3]. 
As the Viking period often had more atmospheric dust loading than did that of Mariner 9, the re- 
flectance spectroscopy technique may even have been incapable of detecting the entire O 3 column 
abundance during much of the Mars year that Viking observed, particularly at high latitudes. 
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IS C0 2 ICE PERMANENT? Bernhard Lee Lindner, At- 
mospheric and Environmental Research, Inc., 840 Memorial 
Drive, Cambridge MA 02139, USA 

Carbon dioxide ice has been inferred to exist at the south 
pole in summertime [1,2], but Earth-based measurements in 
1969 of water vapor in the martian atmosphere suggest that 
all C0 2 ice sublimed from the southern polar cap and 
exposed underlying water ice [3]. This implies that the 
observed summertime C0 2 ice is of recent origin. 

However, Fig. 1 shows that theoretical models of the 
energy budget of the surface that simulate the formation and 
dissipation of C0 2 ice have been unable to preserve seasonal 
C0 2 ice at the south pole and still obtain agreement with 
observations of the polar cap regression and the annual cycle 



Fig. 1. The seasonal recession of the south polar cap as observed 
over the last 20 years [22] and as predicted by models {4,6-9]. (The 
aerocentric longitude of the Sun, L f , is the seasonal index; L t = 0°, 90°, 
180°, and 270° correspond to northern spring equinox, summer 
solstice, autumnal equinox, and winter solstice respectively.) 


in atmospheric pressure [4-10]. This implies that either these 
models improperly treat the energy budget or that C0 2 ice 
from an earlier time is exposed in summer. 

An exact comparison to data is difficult, considering that 
the edge of the polar cap is usually patchy and ill defined 
[18,19], in large part due to terrain that is not included in any 
polar cap model. The edge of the polar cap is also diurnally 
variable since ice frequently forms at night and sublimes 
during the day. There is also some year-to-year variability in 
polar cap regression [20,21]. 

Several processes have been examined that might retain 
the good agreement to observations of the annual cycle in 
atmospheric pressure and to overall polar cap regression, and 
yet allow for better agreement at the south pole, without 
requiring old C0 2 ice. The radiative effects of ozone were 
suggested as important [11], but were shown numerically to 
be unimportant [12,13], However, the radiative effects of 
clouds and dust [12] and the dependence of frost albedo on 
solar zenith angle [14,1 5] do allow for better agreement at the 
pole while maintaining good agreement to overall polar cap 
regression and the atmospheric pressure cycle [16]. Penetra- 
tion of sunlight through the seasonal ice also has a marginal 
positive effect on C0 2 ice stability at the pole itself because it 
allows some solar radiation that would otherwise sublime 
overlying C0 2 seasonal ice to sublime ice within the residual 
polar cap [17]. 

Figures 2 and 3 show my model predictions for polar cap 
regression compared to observations. Before solar longitude 
of 250*, south polar cap regression is predicted to be similar 
to that predicted by earlier models (compare to Fig. 1). 
However, the new model retains C0 2 ice year round at the 
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Fig. 2. Th« regression of the south polar cap, as observed for various 
years (taken from [23]) and as simulated by my model (thin line), as a 
function of the aerocentric longitude of the Sun (L f ). The cap radius is 
that which would be measured on a polar stereographic projection of the 
south polar region; the units of the radius are tactions of the planetary 
radius of Mars. 
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Fig. 3. The regression of the south polar cap, as observed in 1986 
(solid circles), 1971 (crosses), and 1977 (plus signs) (taken from [21]) 
and as simulated by my model (thin line), as a function of the aerocentric 
longitude of the Sun (L $ ). Two-sigma error bars are indicated for the 
1986 data; the errors are smaller for the denser 1971 data and for the 
1977 Viking data. 

south pole. Unfortunately, the model does overpredict polar 
cap extent in early southern summer (see Fig. 2). 

In summary, it appears possible to construct an energy 
balance model that maintains seasonal 0O 2 ice at the south 
pole year round and still reasonably simulates the polar cap 
regression and atmospheric pressure data. This implies that 
the C0 2 ice observed in the summertime south polar cap 
could be seasonal in origin, and that minor changes in climate 
could cause CO z ice to completely vanish, as would appear to 
have happened in 1969 [3]. However, further research 
remains before it is certain whether the C0 2 ice observed in 
the summertime south polar cap is seasonal or is part of a 
permanent reservoir. 
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THE INTERANNUAL VARIABILITY OF POLAR CAP 
RECESSIONS AS A MEASURE OF MARTIAN CLI- 
MATE AND WEATHER: USING EARTH-BASED 
DATA TO AUGMENT THE TIME LINE FOR THE 
MARS OBSERVER MAPPING MISSION. L. J. Martin* 
and P. B. James 2 , 'Lowell Observatory, 1400 West Mars Hill 
Road, Flagstaff AZ 86001, USA, 2 Department of Physics and 
Astronomy, University of Toledo, Toledo OH 43606, USA. 

Seasonal Cycles of Dust, Water, and C0 2 : The reces- 

sions of the polar ice caps are the most visible and most 
studied indication of seasonal change on Mars. Strong, if cir- 
cumstantial, evidence links these recessions to the seasonal 
cycles of C0 2 , water, and dust. These phenomena and their 
interactions will be the subject of an MSATr workshop next 
year titled “Atmospheric Transport on Mars.” Briggs and 
Leovy [1] have shown from Mariner 9 observations that the 
atmospheric polar hoods of the fall and winter seasons are at 
least partially water ice clouds. Around the time of the vernal 
equinox, this water ice may precipitate onto the surface that 
includes C0 2 frosts. The sublimation of the outer edge of the 
seasonal cap begins about this same time, and we begin to 
observe its recession. During the recession of the north cap 
we also observe cirumpolar clouds that are believed to be 
formed by water vapor from the subliming cap [2]. Some 
observations suggest that at least part of the sublimed water 
and/or C0 2 reforms as surface ice toward the cap’s interior. 
This “new” ice is probably the bright component of the polar 
caps that is seen on Earth-based observations. This would 
explain the south cap’s appearance as that of a shrinking do- 
nut during its recession [3]. Near the edge of the shrinking 
cap, dust activity is also evident on the Viking images [4], 
This may result from off-cap winds generated from sublima- 
tion and/or dust that might be released from within or be- 
neath the ices. It has been found that all of Mars’ major dust 
storms that have been observed to date occurred during the 
broad seasons when either the north or south polar cap was 
receding [5], There are short seasonal periods around the be- 
ginning and ending of cap recessions when no major dust 
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The Abundance of Ozone on Mars. BERNHARD LEE LINDNER (AER) 

The most extensive data on martian ozone (O 3 ) come from Mariner 9 obser- 
vations. The Mariner 9 UV spectrometer scanned from 2100 to 3500 
Angstroms with a spectral resolution of 15 Angstroms. O 3 column abun- 
dances inferred from these spectra show strong seasonal and latitudinal 
variability, ranging from 0.2 pm -atm at equatorial latitudes to 60 pm-atm 
over northern winter polar latitudes. Unlike Earth, these O 3 amounts do not 
significantly affect atmospheric temperature ( 1 ). I use a radiative transfer 
model based on the discrete ordinate method to calculate synthetic radiance 
spectra from 2100 to 3500 Angstroms. A factor of three times as much O 3 
is needed to generate the same spectrum for typically observed cloudy, 
dusty atmospheres as for a clear atmosphere. If the scattering properties of 
martian clouds and dust were well known, then their appearance would not 
be a problem, as a model would be capable of retrieving the O 3 abundance. 
However, these properties are not well known, which raises doubts about 
the effectiveness of the UV spectroscopy technique used to measure O 3 
abundances with Mariner 9, and about the accuracy of previously inferred 
O 3 abundances (cloud and dust effects were largely ignored in earlier anal- 
ysis). While scattering and absorption of sunlight by clouds and dust seri- 
ously affect UV spectroscopy, scattering and absorption of sunlight by 
clouds and dust does not seriously affect O 3 abundance on Mars (2). I ac- 
knowledge support under NASA contract NASW-4614. 
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ATMOSPHERIC CHEMISTRY ON MARS 


B.L. Lindner (Atmospheric and Environmental Research Inc . , 840 

Memorial Drive, Cambridge, MA 02139; 617-349-2280; fax:617-661- 

6479) 

The current state of our knowledge of atmospheric chemistry on 
Mars will be reviewed, and differences with the Earth will be 
highlighted. Improvements in modeling work have shown that the 
excessively high atmospheric mixing required by earlier models (eddy 

o 9-1 

diffusion coefficients of 10 cm s ) to explain the atmospheric 
composition observed by spacecraft is no longer necessary. Also, 
recent work has shown that heterogeneous chemistry could be quite 
important on Mars. 

I will focus on the interactions between ozone and clouds and 
airborne dust. The ozone abundance on Mars is sensitive to the 
presence of clouds and airborne dust, in part due to the effects 
clouds and airborne dust have on photodissociative solar radiation. 
Also, the efficacy of the reflectance spectroscopy technique used in 
the past to infer ozone abundance on Mars is questioned due to 
masking by clouds and dust. 

This work is being supported by a grant from the National 
Aeronautics and Space Administration. The views expressed herein are 
those of the author and do not necessarily reflect the views of NASA 
or any of its subagencies . 
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Several fundamental differences in atmospheric chemistry exist between Mars and the 
Earth. The martian atmosphere is primarily CO 2 (95%), with strong vertical mixing, cold 
temperatures (typically 220K), low pressures (6 mb at die surface), high atmospheric dust and 
cloud particle content, and no man-made atmospheric constituents. Earlier difficulties in explaining 
why the atmosphere was not more decomposed into CO and 02 have been placated with models 
using updated reaction rates, ’moderate' eddy mixing of order 10 ? cm 2 $- 1 , and the odd hydrogen 
catalytic cycle (Shimazaki, 1989; Krasnopolsky, 1992). Odd nitrogen and sulfur catalytic cycles 
are of marginal importance, and other catalytic cycles shown to be important in the terrestrial 
atmosphere are unimportant on Mars (Yung et al., 1977; Krasnopolsky, 1992). Currently, much 
work is being undertaken to examine the importance of heterogeneous chemistry (e.g., Atreya and 
Blamont, 1990; Krasnopolsky, 1992), but uncertainties in particle properties make efficiencies 
difficult to evaluate. Also, atmospheric chemistry may significantly alter atmospheric composition 
on climatic timescales, particularly during periods of low obliquity (Lindner and Jakosky, 1985). 

Ozone is a key to understanding atmospheric chemistry on Mars. The O 3 abundance has 
been inferred from UV spectra by several spacecraft, with the most complete coverage provided by 
Mariner 9 (Lane et al., 1973). The Mariner 9 UV spectrometer scanned from 2100 to 3500 
Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 
Angstroms and an effective field-of-view of approximately 300 km^ The only atmospheric 
absorption in the 2000 to 3000 Angstrom region was assumed to come from the Hartley band 
system of ozone, which has an opacity of order unity. Therefore the amount of ozone was inferred 
by fitting this absorption feature with laboratory data of ozone absorption, as shown in Fig. 1. 
Mars O 3 shows strong seasonal and latitudinal variation, with column abundances ranging from 
0.2 |im-atm at equatorial latitudes 
to 60 pm -atm over the northern 
winter polar latitudes [Lane et al., 

1973] (1 pm-atm is a column abun- 
dance of 2.689x10*5 molecules 
cm ~2). However, the O 3 abun- 
dance is never great enough to sig- 
nificantly affect atmospheric tem- 
peratures (Lindner, 1991) or sur- 
face temperatures and frost budgets 
(Lindner, 1990). Figure 2 shows 
some of the inferred O 3 abun- 
dances. 

I use a radiative transfer 
model based on the discrete ordinate 
method to calculate synthetic radi- 
ance spectra. Assuming a constant 
mixing ratio for ozone and no 
chemical or radiative interaction be- 
tween O 3 and clouds/dust, Fig. 3 



Figure 1. Ultraviolet spectrum measured by Mariner 
9 at 57N latitude on orbit 144. To enhance the O 3 
absorption feature, this spectrum was divided by 
one obtained at 20N latitude on orbit 144, where O 3 
abundances are minimal (Lane et al., 1973). 


shows that when typical amounts of dust and cloud arc present that significant underestimation of 
O 3 occurs. A factor of 3 times as much O 3 is needed to generate the same spectrum as for a clear 
atmosphere. If the scattering properties of martian clouds and dust were well known, then their 
appearance would not be a problem, as a model would be capable of retrieving the O 3 abundance. 
However, these properties are not well known, which raises doubts about the effectiveness of the 
current UV reflectance spectroscopy technique used to measure O 3 . 

Spatial and temporal variability in temperature and water vapor account have been claimed 
to account for the scatter of the data points in Fig. 2 (Barth and Dick, 1974). A decrease in 
temperature would result in a decrease in water vapor, if saturated as expected at prevalent 
temperatures. A decreased water vapor abundance decreases the availability of odd hydrogen, 
which converts CO and O into CX >2 catalytically, decreasing the abundance of O needed to form 
O 3 . However, water vapor is a small source of odd hydrogen in the winter polar atmosphere com- 
pared to H 2 , and may not account for most of the variability in Fig. 2 (Lindner, 1988). Masking 
by clouds and dust may also account for some of the observed O 3 variability, because the nature 
and opacity of the clouds and dust in the polar hood change dramatically in latitude and even on a 
day-to-day basis. As the maximum O 3 abundance resides near the surface (Lindner, 1988), 
spacecraft must be able to observe through the entire cloud and dust abundance in order to actually 
see the total O 3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the 
cloud and the airborne dust must be traversed twice; first by the incoming solar flux down to the 
surface, and then once again upon reflection from the surface out to the spacecraft In addition, the 
large solar zenith angles at winter polar latitudes mean several times the vertical opacity of cloud 
and dust must be traversed. Indeed, part of the observed latitudinal variation in 63 in Fig. 2 may 
be due to the inability of the spacecraft to observe through the increasing effective optical depths as 
one goes poleward. 

By using a photochemical 
model which included multiple 
scattering of solar radiation, 

Lindner (1988) showed that the 
absorption and scattering of solar 
radiation by clouds and dust ac- 
tually increased O 3 abundances at 
winter polar latitudes. Hence, 
regions with high dust and cloud 
abundance could contain high O 3 
abundances (heterogeneous chem- 
istry effects have yet to be worked 
out). It is quite possible that the 
maximum O 3 column abundance 
observed by Mariner 9 of 60fim- 
atm is common. In fact, larger 
quantities may exist in some of 
the colder areas with optically 
thick clouds and dust. As the 
Viking period often had more 
atmospheric dust loading than did 
that of Mariner 9, the reflectance 
spectroscopy technique may even 
have been incapable of detecting 
the entire O 3 column abundance 
during much of the Mars year that 
Viking observed, particularly at 
high latitudes. 
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Figure 2. Measurements of the O 3 column 
abundance inferred from the Mariner 9 UV 
spectrometer data during the northern winter, Ls = 
330-360, in the northern hemisphere (Barth, 1985). 



However, observing the UV 
spectrum from the surface would 
greatly reduce the effects of clouds and 
dust, and hence significantly increase 
the accuracy of ozone abundance re- 
trievals. A lander with a simple pho- 
tometer having only 2 channels, one in 
the O 3 absorption band at 2500 
Angstroms and one out of the band at 
3500 Angstroms, would achieve this. 
Other possibilities for measuring ozone 
include solar occultation (Blamont et 
al., 1989), IR observations in the 9.6 
pm O 3 band (e.g., Espenak et al., 
1990), and observations of the O 2 
dayglow at 1.27 Jim, produced by 
photolysis of O 3 (Traub et al., 1979). 
However, further studies of these other 
techniques are required, especially as 
regards the effects of clouds and dust 
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Figure 3. Synthetic spectra as would be 
observed by spacecraft for atmospheres with 
(A) no cloud or dust and 30 jim-atm O 3 , (B) ' 
vertical opacities of dust and cloud of 0.3 and ; 
1.0, respectively, and 30 pm-atm of O 3 , and 1 
(C) vertical opacities of dust and cloud of 0.3 1 
and 1 . 0 , respectively, and 100 pm-atm of; 
O 3 . All cases assume a solar a zenith angle of 
75 (typical for winter polar observations), 
and viewing angle of 60, with azimuth angle 
of 0 (typical for Mariner 9). Polar-cap albedo 
of 0.8 (new ice). 
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Wavelength (in Angstroms) 

Fig. 3. Synthetic spectra as would beobsctved by spacecraft for atmospheres 
with no cloud or dust and 30 pm-atm O, (solid line), vertical opacities of dust 
and cloud of 0.3 and 1 ,0. respectively, and 30 pm-atm of 0 3 (dashed line), and 
vertical opacitiesofdustandcioudof0.3 and I.O.iespeclively.and lOOpm-atm 
ofOj(doUedline).AIIcascsassumeasolarzeniihangleof75°(typicai for winter 
polar observations), and a polar cap albedo of 0.6. 


to observe through the increasing effective optical depths as one 
goes poleward. 

By using a photochemical model that included multiple scatter- 
ing of solar radiation. Lindner (3) showed that the absorption and 
scattering of solar radiation by clouds and dust should actually 
increase O, abundances at winter polar latitudes. Hence, regions 
with high dust and cloud abundance could contain high 0 3 abun- 
dances (heterogeneous chemistry effects have yet to be fully under- 
stood |2.91). It is quite possible that the maximum 0 3 column 
abundance observed by Mariner 9 of 60 pm-atm is common. In fact, 
larger quantities may exist in some of the colder areas with optically 
thick clouds and dust. As the Viking period often had more atmos- 
pheric dust loading than did that of Mariner 9, the reflectance 
spectroscopy technique may even have been incapable of detecting 
the entire 0 3 column abundance during much of the Mars year that 
Viking observed, particularly at high latitudes. The behavior of O, 
is virtually unknown during global dust storms, in polar night, and 
within the polar hood, leaving large gaps in our understanding. 
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THE EFFECT OF POLAR CAPS ON OBLIQUITY. B. L. 
Lindner, Atmospheric and Environmental Research, Inc., 840 
Memorial Drive, Cambridge MA 02139. USA. 

Rubincam ( I ] has shown that the martian obliquity is dependent 
on the seasonal polar caps. In particular, Rubincam analytically 
derived this dependence and showed that the change in obliquity is 
directly proportional to the seasonal polar cap mass. Specifically. 
Rubincam showed 

d ip/dt » 3 x IO- 10 M(t)/M(0) degrees/Earth year ( I ) 

whereqr is the obliquity and M is the mass of the seasonal polar caps, 
with time t of 0 being the present. This expression assumes uniformly 
thick spherical caps with identical angular radii of 45°. However, 
even if a very different polar cap mass distribution is used, Rubincam 
estimates the total uncertainty in the constant in equation (I) to be 
less than a factor of 2. Using the current mass of the seasonal polar 
cap as typical over geologic time. Rubincam calculates that the 
amount that the obliquity would secularly change is only 1.4°. 
Considering that the current obliquity of Mars is 25°. Rubincam 
concludes that seasonal friction docs not appear to have changed 
Mars' climate significantly. 

Using a computer model for the evolution of the martian atmos- 
phere. Habcrlc et al. (2.3) have made a convincing case for the 
possibilityofhugepolarcaps.about I Ox the mass ofthe current polar 
caps, that exist fora significant fraction of the planet's history. Given 
the large uncertainties in input parameters and in the model itself, the 
results must be regarded as speculative. Also, the Habcrle et al. 
results have been unable to favor or rule out a large polar cap 
scenario vs. a small polar cap scenario. ; 

Nonetheless, since Rubincam showed that the effect of seasonal 
friction on obliquity is directly proportional to polar cap mass, a ; 
scenario with a ten-fold increase in polar cap mass overa significant 
fraction of the planet's history would result in a secular increase in ;; 
Mars' obliquity of perhaps 10° (using equation (I)). Hence, the '■ 
Rubincam conclusion of an insignificant contribution to Mars' cli- 
mate by seasonal friction may be incorrect. Furthermore, if seasonal { 
friction is an important consideration in the obliquity of Mars, this ■ 
would significantly alter the predictions of past obliquity as prr- : 
sented by Ward [4-6], Murray ct al. [7J. Ward ct al. (8). Rubincam 
(91. Chao and Rubincam (10). Bills (111. Ward and Ruby (12]. J 
Touma and Wisdom (13J. and Laskar and Robutel (14], That in turn 
would significantly alter the predictions of past climate, which are }.<) 
based on obliquity predictions (15-20). The mechanics of the polar 3 
cap system also depend on obliquity (21-26). If obliquities were | 
often much smaller than at present, that could have implications for > 
past atmospheric composition (27). % 

Given the enormity of the implications, the effect of the polar caps 4 
on the obliquity of Mars should be given more attention and study. 
Perhaps further modeling of obliquity could be used to rule out the .!>; 
possibility of large polar caps for extended times, which would S 
assist modeling of atmospheric evolution. Similarly, modeling of $5 
atmospheric evolution should be given more attention and study # 
because of the implications for obliquity history, and therefore 
climate history. t‘ 
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ESCAPE OF MARS ATMOSPHERIC CARBON THROUGH 
TIME BY PHOTOCHEMICAL MEANS. J. G. Luhmann'. J. 
Kim 2 , and A. F. Nagy-', 'Institute of Geophysics and Planetary 
Physics. University of California. Los Angeles CA 90024-1567. 
USA. 2 KARl,Scoul, Korca.'Spacc Research Laboratory, University 
of Michigan, Ann Arbor MI 48109, USA. 

Luhmann et al. ( 1 ] recently suggested that sputtering of the mar- 
tian atmosphere by reentering O* pickup ions could have provided a 
significant route of escape for C0 2 and its products throughout Mars' 
history. They estimated that the equivalent of C in a - 1 40-mbar CO : 
atmosphere should have been lost this way if the Sun and solar wind 
evolved according to available models. Another source of escaping 
C (and O) that is potentially important is the dissociative recombi- 
nation of ionospheric CO 4 near the exobasc (2). We have evaluated 
the loss rates due to this process for "ancient" solar EUV radiation 
fluxes of 1.3, and 6x the present flux in order to calculate the 
possible cumulative loss over the last 3.5 Gyr. (Earlier estimates of 
loss by McElroy [2] used the present-day rates and thus represent 
underestimates.) The inputs and assumptions for this calculation 
arc the same as used by Zhang el al. [3] for an evaluation ofhistorical 
O escape by dissociative recombination of ionospheric 0,\ We find 
loss rates of C that arc at least comparable to the sputtering loss 
rates, thereby potentially accounting for another 100 mbar or more 
of Mars' original atmosphere. 

References: [1 ] Luhmann J. G. et al. (1992) GRL, 19. 2151- 

2154. (2] McElroy M. B. (1972) Science, 443. [3] Zhang M. H. G. 
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MIGHT IT BE POSSIBLE TO PREDICT THE ONSET OF 
MAJOR MARTIAN DUST STORMS? L. J. Martin', P. B. 
James 2 , and R. W. Zurek', 'Lowell Observatory, Flagstaff AZ 
86001, USA, 'University of Toledo, Toledo OH 43606, USA, 'Jet 
Propulsion Laboratory and California Institute of Technology. 
Pasadena CA 91 109, USA. 

This was done very successfully by the late "Chick" Capen in 
1971, but we now believe that the chance of having a planet- 
encircling storm in any given Mars year is less than 50% [ I ]. Capen 
suggested that these storms begin around the time of perihelion. 
More recent storms have extended this season to nearly one-third of 
a martian year, during the same interval that the south polar cap is 
receding (2). There is no observational evidence that storms of this 
size have occurred outside of that period, although smaller dust 
storms have been observed throughout most of the martian year. The 
circumstances that allow a limited storm to become a runaway or 
encircling storm are not well understood. Seasonal effects are appar- 
ently just one aspect of these circumstances, but apparently a critical 
one. Dust activity seen by Viking near the edges of the receding cap 
and data showing that the cap may be receding at a faster rate prior 
to these storms suggest that the seasonal south cap may be influenc- 
ing dust activity. 

We have also determined that the north polar hood recedes during 
major dust storms, but it is not clear whether impending storms 
might have an effect upon this atmospheric phenomenon. Viking 
images do show local storm clouds near the hood prior to the first 
1977 planet-encircling dust storm, but the hood is such a dynamic 
feature that minor changes may not be meaningful. Wc are, how- 
ever, continuing to analyze these data. 

Several datasets indicate that Mars’ atmosphere was less clear 
before the first 1977 encircling storm, although we cannot discount 
the possibility that this was merely a seasonal change. Data from 
other Mars years are less detailed and comprehensive, but the 1977 
Viking data from both imaging [3] and infrared [4] suggest that dust 
in the atmosphere was increasing prior to the storm. Peter Boyce 
found that, prior to the 1971 planet-encircling storm, there was 
“violet haze" present on Mars. He attributed this to the impending 
storm, which may have been correct, but this condition, which could 
be due in part to atmospheric dust on Mars, is not uncommon at times 
when no storm is on the way. This may also be true for other 
indicators of increasing atmospheric dust mentioned above. 

Capen also believed that smaller, precursor storms occurred be- 
fore a planet-encircling storm. This generally seems to be the case, 
although the data arc not conclusive. These earlier storms certainly 
provide a good vehicle for raising dust into the the atmosphere and 
regional dust storms may be a sign of an impending larger storm. 
However, many of these storms occur without any subsequent dust 
activity, even during the dust storm "season." 

Investigations of dust-storm observations show that that the Hcl las 
Basin is the most active area on Mars for all sizes of storms (2). This 
area is probably their primary dust source. 

Earth-based observations suggest that, during the expansion 
phase of planet-encircling storms, diurnal cycles often begin at 
Hellas, presumably with a new load of dust, as mountain climbers 
return to a base camp for more supplies to be cached along their 
route. Each day the storms carry an increasing supply of dust farther 
to the west, until Hellas is reached from the east, completing the 
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HOW WELL WAS TOTAL OZONE ABUNDANCE INFERRED WITH MARINER 9? Bernhard 
Lee Lindner, AER, Inc., 840 Memorial Drive, Cambridge, Massachusetts 02139 USA 

Ozone is a key to understanding atmospheric chemistry on Mars. Over 20 photochemical 
models of the martian atmosphere have been published, and O3 is often used as a benchmark for 
these models[l-3]. O3 abundance has been inferred from instrumentation on several spacecraft, 
with the most complete coverage provided by Mariner 9[3,4]. The Mariner 9 UV spectrometer 
scanned from 2100 to 3500 Angstroms with a spectral resolution of 15 Angstroms and an effective 
field-of-view of approximately 300 km 2 [4]. The only atmospheric absorption in the 2000 to 3000 
A wavelength region was previously assumed to come from the Hartley band system of ozone[4], 
which has an opacity of order unity at winter polar latitudes[3]. Therefore the amount of ozone 
was inferred by fitting this absorption feature with laboratory data of ozone absorption, as shown 
in Fig. 1[4]. Mars O3 shows strong seasonal and latitudinal variation, with column abundances 
ranging from 0.2 pm-atm at equatorial latitudes to 60 jim-atm over northern winter polar 
latitudes[4] (1 pm -atm is a column abundance of 2.689xl0 15 molecules cm -2 ). However, the O3 
abundance is never great enough to significantly affect atmospheric temperatures[5] or surface tem- 
peratures and frost amounts[6]. Figure 2 shows some of the previously-inferred O3 abun- 
dances[7]. 

A radiative transfer computer model is used to re-examine the Mariner 9 U V spectra. 
Assuming a constant mixing ratio for O3 and no chemical or radiative interaction between O3 and 
clouds/dust. Fig. 3 shows that when typical amounts of dust and cloud are present that significant 
underestimation of O3 abundance occurs. A factor of 3 times as much O3 is needed to generate the 
same spectrum the spacecraft would measure for a cloudy, dusty atmosphere as for a clear atmo- 
sphere. If the scattering properties of martian clouds and dust were well known, then their appear- 
ance would not be a problem, as a model would be capable of retrieving the O3 abundance. 
However, these properties are not well known, which raises doubts about the effectiveness of the 
UV reflectance spectroscopy technique for measuring O3 abundance on Mars. The simulations 
shown in Fig. 3 are repeated for a range in solar zenith angle (50*-90*), ground albedo (0.3-0.8), 




altitude distribution of O3, satellite viewing geometries, and cloud, dust and O3 abundances. A 
factor of 3 underestimation is typical, with greater underestimation for high ground albedo or high 
dust opacities. Even if scattering by clouds is properly accounted for (as previously done with 
Mariner 9 data reduction in [4]), masking by dust can easily result in factor of 2 underestimation. 
Results are not strongly dependent on solar zenith angle. 

Spatial and temporal variability in temperature and water vapor have been claimed to ac- 
count for the scatter of the data points in Fig. 2[8]. A decrease in temperature results in a decrease 
in water vapor, if saturated as expected at prevalent temperatures. A decreased water vapor abun- 
dance decreases the availability of odd hydrogen (H, OH, and HO2), which converts CO and O 
into CO2 catalytically, decreasing the abundance of O needed to form O3. However, water vapor 
is a small source of odd hydrogen in the winter polar atmosphere compared to H2, and may not 
account for most of the variability in Fig. 2[3]. Masking by clouds and dust may also account for 
some of the observed O3 variability, because the nature and opacity of the clouds and dust at winter 
polar latitudes change significantly spatially and temporally. As the maximum O3 abundance 
resides near the surface[3], spacecraft must be able to observe through the entire cloud and dust 
abundance in order to measure the total O3 column abundance. If reflectance spectroscopy is used, 
as on Mariner 9, then the cloud and the airborne dust must be traversed twice; first by the incoming 
solar flux down to the surface, and then once again upon reflection from the surface out to the 
spacecraft. In addition, the large solar zenith angles at winter polar latitudes mean several times the 
vertical opacity of cloud and dust must be traversed. Indeed, part of the observed latitudinal varia- 
tion in O3 abundance in Fig. 2 may be due to the inability of the spacecraft to observe through the 
increasing effective optical depths as one goes poleward 

By using a photochemical model which included multiple scattering of solar radiation, 
Lindner(3] showed that the absorption and scattering of solar radiation by clouds and dust should 
actually increase O3 abundances at winter polar latitudes. Hence, regions with high dust and cloud 
abundance could contain high O3 abundances (heterogeneous chemistry effects have yet to be fully 
understood[2,9]). It is quite possible that the maximum O3 column abundance observed by 



Mariner 9 of 60p.m-atm is common. In fact, larger quantities may exist in some of the colder areas 
with optically thick clouds and dust. As the Viking period often had more atmospheric dust 
loading than did that of Mariner 9, the reflectance spectroscopy technique may even have been 
incapable of detecting the entire O3 column abundance during much of the Mars year that Viking 
observed, particularly at high latitudes. The behavior of O3 is virtually unknown during global 
dust storms, in polar night, and within the polar hood, leaving large gaps in our understanding. 
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FIGURE CAPTIONS 


Figure 1. Ultraviolet spectrum measured by Mariner 9 at 57°N latitude on orbit 144(7]. To 
enhance the O 3 absorption feature, this spectrum was divided by one obtained at 20'N latitude on 
orbit 144, where O 3 abundances are minimal[7]. 

Figure 2. Measurements of the O 3 column abundance previously inferred from the Mariner 9 UV 
spectrometer data during the northern winter, L* = 330-360*, in the northern hemisphere[4]. 

Figure 3. Synthetic spectra as would be observed by spacecraft for atmospheres with no cloud or 
dust and 30 (tm-atm O 3 (solid line), vertical opacities of dust and cloud of 0.3 and 1.0, respec- 
tively, and 30 (im-atm of O 3 (dashed line), and vertical opacities of dust and cloud of 0.3 and 1.0, 
respectively, and 100 (im-atm of O 3 (dotted line). All cases assume a solar a zenith angle of 75° 
(typical for winter polar observations), and a polar cap albedo of 0 . 6 . 



RELATIVE INTENSITY 


ORBIT 144 



I I 1 I I I I I I I I l I I I 


2100 2500 3000 3500 

WAVELENGTH <A) 


70 


60 H 



Relative Intensity 



2100 2400 2700 3000 

Wavelength (in Angstroms) 





N94- 27892 


In Situ Mars Ozone Detector 

Dr. B. Lee Lindner 

Atmospheric and Environmental Research, Inc. 

840 Memorial Drive, Cambridge, MA 02139 
(617) 547-6207 

Dr. Elliot M. Weinstock 
Dept, of Chemistry and Center for Earth and Planetary Physics 
Harvard University, Cambridge, MA 02139 

We propose sending a balloon-borne UV photometer sensor package to measure atmospheric 
ozone on Mars, and this package could be a Discovery Program sensor candidate. Past measure- 
ments of ozone on Mars are highly uncertain, perhaps a factor of 3 or so uncertain (Lindner, 1993), 
due primarily to interference and masking by cloud and dust In-situ balloon measurements would 
avoid these problems, and would provide ‘ground truth’ which remote sensing techniques cannot 
We have explored this approach to measure ozone abundance in the terrestrial stratosphere with a 
balloon-borne UV absorption photometer (described in detail in Weinstock et al., 1986). 
Atmospheric pressures and temperatures and ozone concentrations near the surface of Mars are simi- 
lar to those in the terrestrial stratosphere. 

In brief, the instrument uses 245-nm radiation from a low-pressure mercury discharge lamp 
in a 40-pass white cell (1200-cm path length). The instrument has measured terrestrial stratospheric 
ozone with better than 3% precision and 5% accuracy. The instrument is small (approx. 30cm 
cubed), lightweight (approx. 20 kg), self-scrubbing (i.e., self-calibrating), and has in-flight diag- 
nostics and a low data rate. This balloon-borne photometer is expected to be able to measure ozone 
throughout the polar winter region (poleward of 40* latitude) within the lower 20 km of the atmo- 
sphere but possibly not at equatorial latitudes where O3 concentrations are extremely low (10 9 cm' 3 ), 
although slight improvements may allow for measurement of these low concentrations as well. 

The objectives for this instrument are to observe O3 variability with cloud and dust amount, 
to measure the altitude dependence of O3, and to measure diurnal fluctuations in O3. There are no 
measurements to date of any of these variations. The abundance of ozone during the great dust 
storms on Mars is totally unknown due to the inability to see through the dust. The possibility of 
heterogeneous chemistry of ozone with airborne dust and cloud on Mars is the most active topic in 
Mars aeronomy today, and the question can only be settled through an in situ study such as we pro- 
pose. The dependence of ozone on altitude has never been measured, only the column abundance. 
An improved understanding of the altitude dependence of ozone would also determine the degree of 
vertical atmospheric transport. Because the instrument proposed here contains a lamp, ozone could 
be measured throughout the day/night cycle. Ozone in polar night has never been measured, and 
measurements here would greatly increase our understanding of heterogeneous chemistry and atmo- 
spheric transport. 

Hence, in-situ ozone measurements would significantly improve our understanding of Mars 
atmospheric chemistry and composition and the interaction of atmospheric chemistry with atmo- 
spheric radiation and dynamics. Much interest has been generated over the past 5 years for a Mars 
Aeronomy Mission. This is a low cost alternative which in many respects will learn more informa- 
tion. This mission would complement a proposed Japanese mission to study the chemistry of the 
martian upper atmosphere (Tsuruda, 1991). 
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